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Abstract 

Recent developments of the covariant Bethe-Salpeter (BS) approach with the use of the sepa- 
rable interaction for the deuteron are reviewed. It is shown that the BS formalism allows a covari- 
ant description of various electromagnetic reactions like the lepton-deuteron scattering, deuteron 
electro-disintegration, deep inelastic scattering (DIS) of leptons on light nuclei. The procedure of 
the construction of the separable nucleon-nucleon (NN) interaction is discussed. The BS formalism 
facilitates analysis of the role of the P- waves (negative energy components) in the electromagnetic 
properties of the deuteron and its comparison with the nonrelativistic results. Furthermore the 
covariant BS approach makes it possible to analyze DIS of leptons from the deuteron in a model 
independent way and to extend the formalism to DIS reactions on the light nuclei. 
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1 Introduction 



The study of electromagnetic properties of light nuclei, A < 4, facilitates the construction of the 
theory of strong interactions and, in particular, the nucleon-nucleon interaction (see, for example, |T]]). 
Theoretical research in this field is of topical interest which is reflected in recent review articles 0-0. 
A large amount of available experimental data stimulate a further development of theoretical methods 
which are often restricted to qualitative predictions. The forthcoming experiments are expected to 
provide high precision data which will allow us to explore the region of large momentum transfer in 
elastic, inelastic and deep inelastic (DIS) electron-nucleus reactions. 

These data will be able to furnish qualitatively new information about the fine nuclear structure 
at small distance. This is the reason why the role of the non-nucleon degrees of freedom as mesons, 
A-isobars, quarks etc. on intermediate and high energy phenomena is widely discussed (see, for ex- 
ample. Clear understanding and consistent interpretation of the experimental information is 



not possible without the consideration of the relativistic kinematics of reactions and the dynamics of 
the interaction. For this reason, the construction of a covariant approach and a detailed analysis of 
relativistic effects in electromagnetic reactions with light nuclei become the tasks of the highest priority. 

One can identify three steps for the construction of a theoretical framework of the lepton-nucleus 
interaction. The first step is to introduce the dynamical degrees of freedom parameterizing hadronic 
and nuclear structure. The second one is to construct from these degrees of freedom bound states 
which are hadrons and nuclei. The third step is to find a formalism of the interaction of these degrees 
of freedom with incident particles which are photons 7 and leptons I in the present study. 

In some sense this division is artificial, since the three steps are just different parts of the same 
problem of dynamics of interacting fields to be solved in one consistent approach. The problems are 
interrelated by the underlying physics in such a way that by fixing one of them one can find solutions 
for the others. Nevertheless, such a division helps us to separate the task within a consistent set of 
approximations and to use a phenomenological technique to set constrains on different parts of the 
approach. 

The fact that nuclei consist of bound nucleons introduces a major problem for theoretical description 
of relativistic / — A interactions. The deuteron is naturally the first object in the row of many other 
nuclei, and has received a vast number of treatments within many different approaches. One finds also 
that non-relativistic schemes of calculations are widely employed in the analysis, which can be justified 
for a few particular cases. On the other hand, the consistent consideration of the relativistic bound 



states is offered within the Bethe-Salpeter (BS) formalism p2|, |23[, which makes it most promising for 
the class of the tasks considered in the present review. What is even more important is that the BS 
formalism allows qualitatively a new interpretation of the physics of the relativistic bound state and 
should not be regarded as an alternative scheme only. 

The first practical applications of the BS formalism were based on the three dimensional reduction. 
The first reduction of the BS formalism, the so-called quasi-potential equation, has been made in 



several works The main idea behind these approaches is to fix the relative time of bound 

constituents to a certain value. Since the relative time (or relative energy in momentum space) is 
considered as an unphysical feature of the BS formalism, the different ways of fixing the relative time 
should produce equivalent quasipotential approaches. The comparison of the different quasipotential 
equations was done in review articles ^. This line of realization of the BS formalism allows one to 
take into account some of the relativistic effects, but it looses the general relativistic covariance pSI (see. 



papers [^, ^ also). Moreover, a number of difficulties arises which does not allow one to establish a 
direct link to the non-relativistic calculations, for example, the absence of a non-relativistic reduction of 
arbitrary kernels, the problem of the interpretation of the abnormal parity states. This fact motivates 
us to use the original BS formalism which offers a consistent covariant description of the interacting 
particles and their bound states. The qualitative connection with the nonrelativistic results can be 
made on the level of observables. Such exotic features of the BS formalism as relative time of the bound 
constituents and abnormal parity states receive their nonrelativistic interpretation through the relation 
to dynamical degrees of freedom. Analysis of the relative time dependence gives qualitatively a new 



This review is devoted to the analysis of the three steps considered above within the Bethe-Salpeter 
formahsm and its apphcation to the study of the electromagnetic properties of light nuclei. We empha- 
size the covariant description of the BS formalism by taking the separable interaction, which is still a 
stage of infancy. In particular, the role of the abnormal parity states, in not yet confronted with exper- 
imental data, though the necessity is demonstrated in this paper. The review is divided into 6 sections. 
In section 2, we consider basic properties and definitions of the Bethe-Salpeter formalism. By using an 
example of the NN system we discuss the basic properties of the four dimensional (4D) bound states. 
In section 3, we investigate the simplest process for such a study — elastic scattering of electrons off the 
deuteron. In section 4, we discuss the application of the BS formalism to the problem of the deuteron 
electro-disintegration, which is an inelastic process with finite momentum transfer. The analysis of the 
deep inelastic scattering of leptons on the lightest nuclei — A = 2 ... 4 is presented in section 5. The 
consideration of the DIS reaction in the infinite momentum frame leads to many simplifications in the 
analysis and allows us to understand clearly some of the important physical consequences of the 4D 
consideration and to draw model independent conclusions. Section 6 is devoted to the summary of the 
results of the AD BS formalism. 



2 Bethe-Salpeter Approach 
2.1 Formalism 

In study of processes, involving bound states, in local field theory we have to consider particles which 
are not asymptotically free. The standard field theory reduction technique gives a way to extract all 
the information about physical states contained in the matrix elements to a product of field operators 
(see for example However, this technique is based on the assumption that the physical states 



in the matrix elements can be treated as being asymptotically free, and the interaction is switched 
on adiabatically. Therefore, bound states are excluded from consideration. In order to include bound 
states, the reduction formalism must be supplemented by a procedure allowing expectation values in 
such states to be expressed in terms of vacuum expectation values. 

This problem is solved in the nonrelativistic field theory by introducing an external classical field 
which allows the bound state dynamics to be described as a particle motion in a potential well. As 
a result, the calculation of the expectation value in a bound state reduces to the calculation of the 
expectation value in a one-particle state, and binding effects are taken into account by introducing the 
momentum distribution of this particle. This simple approach can be used to obtain amplitudes of the 
lepton scattering off bound states in the form of convolution. However, in this case, the role played by 
relativistic corrections and off-shell effects remains unclear. In quasipotential approaches, the solution 
of the bound state problem reduces to deriving an analog of the Shrodinger equation with a covariant 
three-dimensional potential |2^-[^. The calculation of expectation values in bound states should be 



essentially the same as in the nonrelativistic case. However, in contrast to nonrelativistic approaches, 
the quasipotential method allows a qualitative study of the role of the relativistic and the off-shell 
effects. 

A method of calculating the expectation values of T-products of local operators in bound states 



was suggested in |^ . The essence of the method is that the expectation values in bound states are 
expressed in terms of the vacuum expectation values of a T-product of local operators and the matrix 
elements for the transition between the vacuum and the bound state. We shall consider the application 
of this method to processes involving a bound state of n nucleons. 



2.2 The Bethe-Salpeter Amplitude 

A consideration of a scattering process is based on the analysis of the matrix element of a T-product 
of local operators rii{yi) . . . rjkiyk) in the initial A and final A' bound states. : 



P and P' are the initial and final momenta of the bound state while the a and a' are full sets of 
their discrete quantum numbers. The local field operators, rjiiiji), are current operators determining the 
nucleon interactions with external fields. To express this matrix element in terms of vacuum expectation 
values we shall use the fact that in a definite kinematical region the joint propagation of n interacting 
nucleons occurs via the formation of a bound state. In this case the first term of the series expansion 
of the n-nucleon Green's function in intermediate state has the form: 

(0|T(^(a;i) . . . ^{xn)^{x\)) . . .^«)|0) = 
d^P 



5^(0|T(vI/(xO...^(x„))|A(a,P)) 



X 



(27r)3 

x{A{a, P)\T{^{x[) . . .^«))|0)e(mm{|xo,|} - max{\x',^\}) + .... (2) 

The functions and are nucleon field operators in the Heinsenberg representation. The 6 function 
is arising due to the expansion of the T-product in the matrix elements. This function ensures that the 
causality condition is satisfied: 

min{xoi} > max{x'Q-}. (3) 

The extreme values (minimum and maximum) of the nucleon coordinates can be determined by intro- 
ducing the average coordinate of n fields conjugate to the total momentum P, which in the case of n 
particles of the same mass has the form, 

n 

X = ^ . (4) 

n 

The nucleon coordinates relative to this point — X. As a result, the maximum and minimum 

coordinates can be defined as 

maxjxoi} = Xq + max{\xoi\}, 

min{xoi} = Xq — max{\xoi\}. (5) 
Using the integral representation of the 6 function 

+00 

0{xo) = ^ / ^-,dPo. (6) 

— oo 

we obtain the necessary expression relating the vacuum expectation value of the nucleon operators and 
the matrix elements for the transition from the vacuum to the bound state: 

(0|T(V|>(X1) . . . ^{Xr.mx\) . ..^{X'M = E / ^^^"^^'^"''^^''""''''^ ^ 

..,i(Po-i.)(n^a.modHma.mui) (0|T(^(xi) . . {xn))\A{a, P)) {A{a, P)\i:{^ {x[) . . . vl'K))|0) 

Pq-E + i5 

This expression is inconvenient because of using two sets of coordinates simultaneously, {xj} and {X, Xi). 
To avoid this, we shall change over to using the second set of coordinates everywhere in the calculations. 
Owing to the translational invariance, the following relation holds for the fields 

^(x + x) = e*^"*(x)e-'^". (8) 

Replacing Xj by X + Xj in the operators \E' and ^ and shifting by X, using transformation (^, we 
obtain the matrix elements for the transition from the vacuum to the bound state in the space of 
relative nucleon locations: 

(0|T(e*^^vl>(a;i)e-^^^ . . . e'^^vl>(a;„)e-*^^)|A(«, P)) = . . . f „)e-^(^^°-P-^), 



Although the functions ^^^"''^^ formally depend on n variables, only n — 1 of them are independent, 

n 

owing to the equation = 0. As a result of these transformations, Eq. (|^) takes the form: 

i 

Since the integral with respect to Pq is determined by the behavior of the integrand near the pole 
Po = E, we have omitted the exponential factor exp{i{Po — E){max{\xoi\} + max{|xoj|})}. 
The unknown functions <|)^("'^) and <|)'^^"''^'' 

. . . x„) = (0|T(M/(xi) . . . <il{x^))\A{a, P)), 

...<) = {A{a, P)\Tmx[) . . .^K))|0), (11) 

introduced above and entering into the matrix elements for the transition from the vacuum to the bound 
state in Eq. (|n|) describe the nuclear state in terms of the degrees of freedom of the virtual nucleons. 
These functions are called the Bethe-Salpeter amplitudes. They give the solution of the fundamental 
problem — the expression of the expectation values in bound states in terms of the vacuum expectation 
values. 



2.3 Analysis of the Matrix Elements in the BS Formalism 

In order to explain how the matrix element (||) is related to the nucleon Green's functions and the 
Bethe-Salpeter amplitudes, we consider the matrix element 

(o|T(M/(xi) . . . ^{xMyi) ■ ■ ■ Vk{yk)Wi) ■ • •^«))|o) (12) 

near the singularity of the n- nucleon bound state at P^ = M^. We expand the time-ordered product 
in the matrix element (|l|) in a product of matrix elements of \E', r], \E'. In order to do this, we choose 
the maximum and minimum zeroth components from the set {xi} in accordance with (j^) and keeping 
only terms which correspond to the first term in @ we write the T-product in the form: 

T(*(xi) . . . ^ixn)viiyi) ■ ..VkivkWix',) . -MO) =_ 

= T(vi/(xi) . . . ^{xn))T{r]i{yi) . . . Vk{yk)m^{x[) . . . ^{x'J) x 

xe{Xo -Yq- max{\xoi\} - max{\yoi\})e{YQ - X'q - max{|£oJ} - max{\yoi\}) + 

+T(^(x;) . . .^K))T(r/i(yi) . . .fer/(|/,))T(vI/(xi) . . . ^(x^)) x 

x6'(Xo - Yq- max{\xQi\} - max{\yoi\})9(Yo - Xq - max{\xoi\} - max{\yoi\})- 

Inserting a complete set between T-products, we rewrite Eq. (|l]) as a sum over states from the 
complete set: 

(o|T(vi>(xi) . . . ^{xn)viiyi) ■ .-VkiykM^'i) ■ • •^«))|o) = 

= 5^(0|T(M/(a;i) . . . v|/(x„))|i?)(i?|T(r7i(yi) . . . Vk{yk))\R') x 

R,R' 

x{R'\T{^{x[) . ..^{x'J)\0)e{Xo -Yo- max{\xoi\} - max{\yo,\}) x 

x9{Yo - Xq - max{\xQi\} - max{\yoi\}) + 

+ J](0|T(^K) . .Mx'^))\R){R\Tiv,iy,) . . .Vkiyk))\R') x 

B,R' 

X {R'\T{^{xi) . . . ^{xn))\0)e{X'^ -Yo- max{\x'oi\} - max{\yo,\}) x 



N/ CiiV^ 'ry-trtnr'S I t'^ . I \ 'ry~trtnr'S 



where denotes a summation over discrete quantum numbers and an integration over continuous 

R,R' 

variables. 

Since we are interested in the behavior of Eq. (|I|) near the pole at = M^, it is sufficient to study 
the contribution to Eq. ([13[) from the lowest bound state corresponding to this pole. The bound state 
corresponds to the ffist term of the series ( ]T3| ) with R = A{a, P) and R' = A'{a', P'). Using the integral 
representation for the 9 function, we obtain the following expression for the matrix element (p!2|): 



(0|T(vl/(xi) . . . ^(a;„)r/i(yi) . . . . . . ^K))|0) =J]y^^(0|T(*(xi) . . . ^{x^))\A{a, P)) 

(A(a,P)|T(r^i(yO . . .r,u{yu)\^{a\P'))){A{a\P')\T(^{x\) . . .^«))|0) x (14) 

dPQ (IPq e^-iiPo-E){Xo~Ya)^-~i{Pl^-E'){Ya~X[;)^-~i{Po~E){rnax{\x^^^^ 



{271) (27r) (Po - ^ + i5) {Po -E' + i6) 

Using Eq. (j^) to get over to the relative variables, we rewrite this expression as 

(o|T(vi/(xi) . . . ^{xn)vi{yi) . . . Vk{yk)Wi) ■ -MOm = (15) 

Ef d P d P ^j(Po„g)(maa:||xo,||-max||go,||) -i(P/,-gO(max||x;,J|-max||yo,|}) 
/ (27r)4(27r)4 



a,a' 



{Po-E + t6){P^-E' + t6) 

On the other hand, owing to the unitarity of the ra-nucleon Green's function related to the amplitude 
(I) as 

G2n{Xl ...Xn, x[, . . . x'J = (0|T(^(xi) . . . . . . ^«)) |0), (16) 

the matrix element (0) can be rewritten identically as 

(o|T(vi/(xi) . ..^{xMyi) ■ ..VkivkMO ■ • •^«))|o) = 

/ 74 74 j4 / ^4 / 74 // 74 // 74 /// 74 HI r~< / \ 

a Z\ . . . a Zfitt z-^ . . . a z^a z^ . . . ct z^a z^ . . .a z^ '-J2n\X\ . . . Xn^ zi, . . . Zji) 

xG,^{z, . . . ^1, . . . 4)(0|T(vI/(^0 • • • ^{^'nWiVi) . . . rik{yk)^{z'i) . ..^{zlW) 

j d'^Zi . . . d'^Znd'^z'C . . . d'^z'^G2n{Xl ■ ■ ■ Xn, Zi, . . . Zn) 

xG2n+k{zi ...Zn,yi... yk, z'l', . . . z'")G2niz'l' , . . . x[, . . . ), (17) 
where (j2n+fc is the truncated Green's function defined as: 

G2n+kiXi . . . Xn, Z/l • • • Z/fc, x[, . . . x'„) = 

= J d'^Zi . . . d'^Znd'^z[ . . . d'^z'^G2n{Xl ■ ■ ■ ^n, ^fcj • • • ^n) X 

X (0|T(vl/(^0 . . . ^izn)viiyi) ■ ■ ■ Vkiyk)^iz[) . . .^iz'J)\0)G,M. • • • <, ^'i, •••<), (18) 

where index n is related to the nuclear operators and index k is related to the external field operators 
77. Taking into account the behavior of ra-nucleon Green's function (^) near the pole at P^ = M^, we 
compare this expression with ([151) : 



5^ / / d^^---dzndz[...dz'^e 



Pa-E + i6 



x'^in^kyZi . . . Zn.,yi . . .yk, ■ ■ ■ z^) 



E 



Pq-E' + i5 

d^P d'^P' -i(^p^-E){max{\xo,\}-max{\yo,\}) -i{P^-E'){max{\x'^^^^^ 



(27r)4 (27r) 



Multiplying both sides of the integrand in (|1^) by {Pq — E){Pq — E') and passing to the limit Pq ^ E 
and Pq ^ E', we obtain an expression for calculating the expectation value in the bound state of the 
T-product of a set of local operators: 

{A{a, P)\T{viyi) ■ ■ ■ viyk))\A\a', P')) = [ dz, . . . dzndz[ . . . dz'^'^^^"'''\z,, . . . z„) 



G^n+ki^i ...zn,y,...y,,z[... z'J<l>^'^^'''"\z[, . . . <). (20) 

This expression relates the scattering amplitude of n-nucleon bound state with irreducible Green's 
function G2n+k, describing scattering on n virtual nucleons, to the BS amplitudes and 
describing the nuclear states in terms of the nucleon degrees of freedom, the equation of which is to be 
found. 



2.4 The Bethe-Salpeter Equation 

The relation between the BS amplitude $^("'^) and the n-nucleon Green's function is established 
in Eq. (p!0|). Thus, having an equation for G2n, we obtain an equation satisfied by the BS amplitude 
<j)-4(a,P) this task, however, it is insufficient to know G2n only perturbatively, since the analysis of the 
behavior of the Green's function near a bound-state pole requires summation of the entire perturbation 
series. Let us therefore examine what general equations for G2n can be obtained without invoking to 
any perturbation theory. 

The propagation of a free nucleon from a point xi to a point X2 is described by the free Green's 
function S'(i) satisfying the equation of the form ||55|] : 

(i4i - m)S{i){xi,X2) = S{xi - X2), (21) 

where m is the nucleon mass. In the case of a nucleon interacting with its own field, a term taking 
into account the self-interaction appears on the right-hand side of the equation for the exact Green's 
function G2, 

{id^-, - m)G2{xi,X2) = S{xi - X2) + j dx[G2{xi,x[)G2{x[,X2), (22) 

Comparing (PT|) and (P^, we see that the function G2{xi,x[) satisfies the Dyson equation, 

G2ixi, X2) = S~^^{xi, X2) - G^^ixi, X2), (23) 

i.e., it coincides with the one-nucleon irreducible self-energy part. The Green's function S^2) describing 
the joint propagation of two physical nucleons which do not interact with each other satisfies an equation 
of the form: 

(i4i -m*) (id^^ -m*)S'(2)(xi,X2,i/i,?/2) = 5{xi - ?/i)(5(x2 - ^2), (24) 

where m* = m + G2{xi,X2)- Inclusion of the interaction between the nucleons leads to the appearance 
of an additional term on the right-hand side: 

{idxi - m*) (g) {id^2 - m*)G4,{xi,X2, yi, 1/2) = 
= 5(xi-yi)5(x2-y2) + f dy\dy'oGA(xi,X2,y\,y'o)GA(y',,yLyi,y2). (25) 



Comparing (plf) and (p5|), we obtain the two-particle analog of the Dyson equation, namely, the inho- 
mogeneous Bethe-Salpeter equation: 



G4(xi,X2,?/i,?/2) = 5'(2j(xi,a;2,?/i,?/2) - G^^{xi,X2,yi,y2), (26) 

In analogy with the one-nucleon case, the function describing the interaction between the nucleons 
coincides with the irreducible self-energy part of the two-nucleon system. Generalizing to the case of n 
nucleons, we obtain the equation: 

G2n{Xl ■ ■ ■ Xn, X^ . . . X^) = ^(^n) ("^1 " ' ' ^1 • • • ^n) ^2n (-^l • • • -^n, Xi . . . X^) , (27) 

where the function S(^n) is the direct product of ra-nucleon propagators: 

^(„)(xi . ..Xn,x[ . . .<) = (0|T(^(xi)^(x;))|0) ® . . . (0|T(*(x„)^(x'J)|0). (28) 
Using Eq. (|27| ) for G2n, we obtain the integral equation with G2n as the kernel: 

G2n{Xl ■ ■ ■ Xn, X^ . . . Xn) = ^{n) {Xl ■ ■ ■ Xn, X^ . . . X^) ~\~ 

J dZi . . . dZndz'i . . . dz'^S(n)iXi . . . Xn, Z\ . . . Zn)G2n{Z\ . . . Zn, z[ . . . z'^) 

xG2n{z[...z'n,x[...x'n). (29) 

Thus, the exact n-nucleon Green's function is the solution of the integral equation which relates the 
two unknown Green's functions G2n and G2n to each other. 
There are several ways of solving this problem: 

• dispersion method with Nakanishi integral representation of perturbation theory; 

• separable ansatz for G2n', 

• perturbative method. 

In analogy with the Dyson equation, the Bethe-Salpeter equation can be studied by using the 
technique of dispersion relations. This can be realized by introducing a generalization of the spectral 
representation for the exact one-particle Green's function for the case of n particles — the Nakanishi 
integral representation of perturbation theory |^ . This approach has been used successfully for solving 



the Bethe-Salpeter equation in the case of two scalar particles |?^, 33]. 

On the other hand, Eq. ( ^9]) offers an excellent possibility to model the n-nucleon Green's function 
if G2n is introduced exphcitly. Both the separable ansatz and perturbative methods are related to this 
strategy. 

In the case of a separable form for the kernel of Eq. (^), we write: 

N 

G2n{zi . . .Zn,z[. . .Z'J = ^ KjQi^Zi . . . Zn)gj{z[ ...<). (30) 

i,3 

In this case the problem of solving the integral equation is replaced by the problem of solving a system 



of linear equations. This approach has been used successfully to describe a two-nucleon system 
pO[] . Recently, the combination of the approaches based on the use of a separable potential and the 
spectral representation taking into account the analytic properties of the two-nucleon Green's function 
is demonstrated to serve as the foundation for the construction of a relativistic separable ansatz for the 
function G2n • 



Most commonly used form of G2n can be obtained by perturbative methods. Let us consider the 
iterative solution of Eq. (p9|). We take the zeroth iteration as 



Substituting this expression into Eq. (^), we obtain the first iteration: 

+ J dZi . . . dZndz'i . . . dz'^Sl^n) {^1 ■ ■ ■ Xni Zi . . . Zn)G2n{Zl . . . Zn, z[ . . . z'^) X 
xS(n){z[...z'^,x[...x'J. (31) 

In the course of successive iterations we obtain the expansion of the exact n-nucleon Green's function 
in an infinite series in powers of G2n'- 

G2ni,Xl ■ ■ ■ Xji, X-^ . . . X.^'j ^(n) (-^1 • • • -^nj X^ . . . ~l~ 
oo „ 

+ / dZi . . . dZndz[ . . . dz'^ . . . dz^^'^ . . . dz[^^ S{n) {Xi . . . Zi . . . Zn) X 

m 

XG2n(^l • • • Zfij Zi . . . Z^)S{n){^Z-y . . . Z^, Xi . . . X„) ... X 

xG^nizi""-'^ . . . zlr-'KzT^ . . . zt^^)Sazi-^ . . . zt\ x[... x'J. (32) 

On the other hand, the function G2n can be expanded in a perturbation series in a specific nucleon- 
nucleon interaction model (for example, the meson-exchange model): 

G2n(^Xi . . . Xni X^ . . . X^) ^ ^ ^2n('^l ' ' ' "^"5 "^1 ' ' ■ "^n)) 

i 

^2n ("^l • • • ■^m X^ . . . X^) = 5'(„)(xi . . . X^, X^ . . . X^. (33) 

Comparing (^) and (|33D , we obtain: 

G2n{Xl ■ ■ ■ Xjii X^ . . . X^) = 

°° 1 /" 1 

= ^ ^ ^ _^ ^ / dZi... dZndz[ . . . dz'^G^"''\^ {Xi . . . Xn, Zi...Zn) X 

m mi+m2=m 

XG("^+^)2„(^1 ...Zn,z[... 4)G("^^)2n (^1 . . . z'^, x[ . . . x'J. (34) 



In the meson-exchange model of A^A^ interaction, the first term of the series (m = 1) corresponds to 
the one-boson exchange approximation in the kernel of (p9D (the ladder approximation). 

We substitute Eq. (p!0| ) into Eq. (p9|), multiply both sides of the resulting expression by (Pq — E), 
and take Pq ^ E. We obtain: 

<l>^("'^)(xi,...x„)= (35) 

dZi... dZndz[ . . . dz'^Sl^ri){Xi ...Xn,Zi... Zn)G2n{Zl . . . Z^, z[ . . . z'^)^'^^'^'^\z[, . . . z'^). 

Thus, the matrix element for the transition between the vacuum and the n-nucleon bound state satisfies a 
homogeneous integral equation with kernel G2ni which is related to the exact ra-nucleon and one-nucleon 



Green's functions by (]27|) . 

We shall need Eq. (|35| ) for the rest of the calculations. By means of the Fourier transform, the 
Bethe-Salpeter amplitude can be rewritten in momentum space as: 



$^(")(pi ...pn) = j d\i . . . rf^x„e'^>'"»$^("'^)(xi 



x„), (36) 



where Pi is momentum of z-th nucleon, Pi ~ ^- more convenient to use a set of the variables 
which includes the total momentum P explicitly. We use the set of the momenta {P, k{\^^^ where 
ki = Pi — P/ n is the nucleon relative momenta. The momentum fc„ is not included in the set because 
it is not independent, kn = — K. In terms of this set the expression ( |36D can be written as: 

^^^''\P,k,...kn-l) = ! d^X^...d^Xne'^^^^^^^^^''^^\x^...Xn). (37) 



The BS equation in the momentum space takes the form: 



$^(-)(P,A;i,...A;„_i)= (38) 
5(„)(P, h,... K_,) / ^ ■ ■ ■ ^^^2n(P, h... K.u k[... K^,)^^^-\P, k[, . . . k'^_,). 

Since (|38D and (|35|) are homogeneous integral equations, the BS amphtude is defined up to some 
constant. In order to determine this constant we consider the expectation value of the baryon current 
at zero momentum transfer, 

(A(a,P)|J^(0)|A(a',P)) = znP^5"'"'. (39) 
Using Eq. (|20|) we obtain the normalization condition for $^("): 

d%^d^ d^^^i-\p^ . X (40) 

J (27r)4 (27r)4 {271^ (27r)4 K , i, n ij k J 

xG2„+i^(g = 0,P,k^... A;„_i, k[... k'^_,)<l>^^'''\P, k[, . . . k'^_,) = mP^d'^^''' . 

Using the fact that at zero momentum transfer the exact truncated photon-n-nucleon vertex G2n+i^ is 
related to the derivative of the n-nucleon Green's function with respect to the total momentum, 

— d 

G2n+i^{q = 0, P, fci . . . A;„_i, k[... k'^_^) = --^G^^XP, ki... kn-i, k[... A;;_i), (41) 
and expressing with the help of Eq. (^), we obtain the normalization condition: 

J (27r)4 (27r)4(27r)4 (27r)4 ^ ' ^' " 

$^('^')(P, A;;, . . . A;;_i) = -znP^(5"'"'. (42) 



X 



d — 

+ -gJ^G2n{P, ki . . . kn-l, ki - . . k^_i) 

We conclude this section by noting the important features of the BS amplitude: 

• the BS amplitude depends on the zeroth component of the relative coordinate (relative time) of 
the nucleons, which, according to (0), is refiected in the dynamical observables of the n-nucleon 
bound state. In momentum space this leads to a dependence on the zeroth component of the 
nucleon relative momentum (relative energy). The relative time dependence is manifested as 
observable effects in DIS of leptons which will be discussed in section ^, 

• the analytic properties of in (p!OD are related with singularities of Green's function G2n- 
This connection can be used to derive nonperturbatively the kernel of the BS equation |ET|. In 
section p.8.2| we will consider this in details. There are poles associated with the external nucleon 
propagators, cuts in the relative momenta, and poles associated with the various bound states 
formed either by several or by all nucleons. The latter ones are isolated in ([T0|) and therefore do 
not contribute to <l>"^("). The poles associated with a bound state of n-nucleons (n < A) can be 
isolated by special procedure [0, |5^ , discussed in section |5l^ , where it is applied to the study of 
DIS of leptons off light nuclei. The consideration of these poles gives the nuclear amplitude in 
terms of nucleons and lighter nuclei amplitudes. The first type of singularities can be explicitly 
isolated by introducing the BS vertex function: 

S'(„)(P, fci . . . kn-i)T^^''\P, k,... kn-l) = $^(")(^, k,,... kn-l), (43) 



which is widelv used in this review. 



2.5 Basic Properties of Two-Nucleon BS equation 

We consider now the two-particle case of the BS equation, which allows one to understand its basic 
properties in detail. Starting from the formula (p9|) with n = 2, 

G^4a/3,75(a;i, a;2, x\, X2) = 5'(2)^^ X2, x[,X2) (44) 

f ^ - 



where we have introduced explicitly the spinor indices noted by Greek letters. The repeated spinor 
indices are assumed to be summed up. The functions G4 and G4 are the exact and the truncated two- 
nucleon Green's functions, respectively, and 5(2) is the Green's function of two noninteracting nucleons, 
and equals to the direct product of full one-nucleon propagators. It is a widely used assumption to omit 
self-energy part in one-nucleon propagators. In this case 

^map,'y5^^^^ ^2, X[, 4) = S^^^iXi - x[)S^^^{x2 - 4). (45) 

To write the BS equation in momentum space, we take Fourier transform and introduce 

V{P,k',k) = J d'^xid'^X2d^yid'^y2G4:{xi, X2,yi,y2) X (46) 

exp + tk ik , 

where P is the total momentum, k and k' are the relative 4-momenta of the two nucleons before and 
after the interaction. They are connected with 4-momenta of first (pi) and second (^2) particles: 

P = Pi+P2, k = {p,-p2)/2, (47) 
Pi = P/2 + k, p2 = P/2 - k. 

We introduce the function V{P] k' , k) in momentum space for the kernel (j4(a;i, X2, x'^^, Xg) of the BS 
equation (pOf) . Similar formula could be written for functions G4{xi,X2,x[,x'2) and S(^2){xi,X2,x[,x'2). 
Thus, the BS equation for full Green's function of the two-nucleon system can be written as, 

G^p.4P, k', k) = So.,{^ + k)Sps{j - k)6^'\k' - k) (48) 

+^Sc..{^ + k')Spxi^-k') J 0^^V,,,,^iP,k',k")G,^,,siP,k",k), 

where one-nucleon propagator Sapip) with assumption (^) has the form, 

S^p{p) = [l/{p.^-m + ie)]^p. (49) 

Here m is the mass of the nucleon, p ■ 7 denotes p^'y^, and 7^ are Dirac matrices. The Greek letters af3 
denote the component in the 4x4 matrix. 

Introducing two-nucleon T-matrix by equation 

P P f d^k" 

Saa{^ + k')Sf3p{— - k)T^p^^s{P, k', k)= j j^^Ga,p,e\{P, k' , k")Vex,js{P, k", k), (50) 

we can write the BS equation for the T-matrix as 

T,^,^,(P, k\ k) = V^f^^^siP, k', k) (51) 
+^ j ^K./3,.A(P,A;',A;")5,,(P/2 + A;")5Ap(P/2-A;")T,,,^5(P,/e",fc). 

To solve the BS equation, we should assume some form for the interaction kernel. Considering a 
model with exchange particle (for instance, meson with mass /i) interacting with nucleons we could 



1. 



if the two-nucleon system forms a bound state, T-matrix has a simple pole in the total momentum 
squared (s = P^) at the point corresponding to the mass of the bound state s = M^; 



2. in the region s > 4m^, T-matrix has so-called unitarity cut which corresponds to the elastic 
nucleon-nucleon scattering (A^A^ NN); 

3. in the region s > {2m + nfi)'^ , T-matrix has cuts which correspond to the inelastic nucleon-nucleon 
scattering resulting in the production of n mesons with mass fi [NN — > NN{nfi))] 

4. in the region s < (4m^ — [nfiY), T-matrix has cuts which correspond to the inelastic nucleon- 
antinucleon scattering (in a cross-reaction) with production of n mesons with mass fi [NN —>■ 
NN{n^)) (so-called left-hand cuts). 

Another choice of the kernel (for instance a widely used separable form) leads to analytic properties 
different from the ones considered here. It will be discussed in section |2.8| . 

The equation for the BS amplitude in momentum space can be written by using Eq. (|55D. We use 
A{a) = JAi because the total momentum J defines the bound state in the two-nucleon case. 

^i^{P, k) = zS^,{^ + k)Sp,{^ ^K,p,a(P, K k")^i^{P, k"), (52) 

and the normalization condition (^) takes the form: 

74 U ^4 z 



5^^\k - k')S-\P/2 + k)S'\P/2 - k) 



^ {S{P/2 + k)S{P/2 - k)} S-\P/2 + k)S-\P/2 -k) + 



dP^ 



$^^'(P, k) = -2tP,,6^^'. (53) 



+ ^np,k',k) 

If the interaction kernel does not depend on the total momentum P, then Eq. (|53| ) becomes 

Wf'""^^' {SiP/^ + k)S{P/^ - k)} T'-^{P, k) = -2^P,6^^\ (54) 

where r"''^(P, k) is the two-nucleon vertex function defined as 

Sa,{^ + k)Spsi^ - kW.siP, k) = <I>,^(P, A;). (55) 



2.6 Partial-Wave Decomposition of the BS Amplitude 

In order to solve the BS equation and to calculate the cross sections of the electromagnetic reactions 
with a two-nucleon system, we use the partial wave decomposition of the BS amplitude separating 
the radial and the spin-angular parts. The two representations for the partial-wave decomposition are 
considered. 



2.6.1 Direct Product Representation 

In the direct product representation, we determine the two particle spinor basis in laboratory frame as 
UP^ik) ® UP^{—k), where /i is the spin projection, pi 2 is the so called p-spin [^], which distinguishes 
the positive and negative energy states. Both the positive and negative energy states are necessary in 
order to prepare a complete set for the two-particle bound state. The spinors U^{k) are connected with 
the Dirac free spinors, Us{k) and fs(fc), as 



UP(k) = I , ^ (56) 



The Dirac spinors are determined as |5 

Us{k) = L{k)usiO), Vsik) = L{k)vsiO 
and the boost operator for a particle with spin 1/2 and mass m is [|126 

m + k ■ 77o 



(57) 



(5J 



Here k = {Ej^, k) is the 4- momentum of a particle on mass shell, Ej^ 
particle. In laboratory frame the spinors can be written as: 



+ is the energy of the 



uJO) 



Xs 





vJO) 





X-s 



where Xs are two-component Pauli spinors. The normalization conditions are: 



m 



Us{k)us'{k) = —S, 



m 



E 



k 



Vs{k)vsi{k) = -—6ss'- 



E 



k 



(59) 



The BS amplitude $'^'^(P, k) of the two-particle system with total angular momenta J and projec- 
tion A4 in laboratory frame can be written as 



^■'^ip,k)= J2 <PjLsp,,,ih,\k\)y''^"'^'^''^ik), 

LSpip2 



(60) 



where P = (M, 0) is the total momentum, k is the relative momentum of the two-particle system 
{k = {ko, k), ko 7^ Ef^). Here JLSpip2 is combination of quantum numbers of total angular momentum 
J, orbital momentum L, spin S and p-spins. We define the spin ang ular function as yJLSpip2,M(^ky_ 

yjLSp^p^M^f^^ ^ J2 {Lm,Sms\JM)CrSi'r^2\Sms)YL^,ik){U^^^ (61) 

where (.|.) are Clebsch-Gordan coefficients and k = k/\k\. The spin-angular function y-^^''^PiP'^^^(^k) is 
a matrix, (1 x 4) ® (1 x 4) = (1 x 16) in spinor space. The spinor indices ajS specify the component of 
this matrix. 

The orthogonalization condition for the spin-angular functions is 



/ rffifc3^i^^^^''^'^\fc)3^i^'^'^'^''-^'(fc) = S,L,6, 



'MM' ^SS' Spip[Sp2p'2- 



(62) 



Where: dfl^ = dcf)^ dcosO^ and the conjugated spin-angular function can be obtained by substitution, 
Us^{k)^Ut{k). 

We can write the inverse propagators [S{P/2 + A;)] [S{P/2 — k)] ^ in laboratory frame as: 



[S{P/2 + k)]-^ = P-^/2 + k 



■ '-f — m 



(63) 



(pi ■ 7 - m)S''y + (p2 ■ 7 + rn)S. 



-1 



[S{P/2 - k)]-' 



P 



(1)" 



1 



2E 



k 



2 — A; ■ 7 — m = 

{p2 ■ 'y - 171)3^^^ + (pi ■ 7 + m)S'+^ 



where pi = {Ef^, k),p2 = {Ef^, —k), s = P^, and functions Sp, with p = ± are 
q(i) _ 1 //', /7/0 I i.„ _ ^Z7. ^ c(2) _ 1 /7/0 _ z,„ _ „ J7. ^ 



Using expressions (|56D , (|57D , (pSf ) and the Dirac equation, we can write 

[S{P/2 + k)]-' U^^ik) = p5,« f/;(-fc), 
[5(P/2-A:)]-it/;(-fc)=p5/^)f/;(fc) 

Here, we can write the expansion of the BS vertex function r^^(P, k) as 

r;:5"(^, ^) = E (-1)^ ^jL5p,p2(fco, |fc|) y^'^^'^'^i-k). (65) 

LS'pip2 

Here the radial part gLSpip2 of the vertex function is connected with the radial part 4>LSpiP2 of ths BS 
amplitude (|6^) through 



'rJLSpip2 

It is further convenient to introduce symmetrical notation for the positive and negative energy states. 



We define states with total p-spin g: \g >= J^p^pii 


\pi 


\p2\Qli'e)\\pi > ® 


\p2 


>, as 




(+) = |11>, (- 




|1-1>, 






(67) 


(e) = 110 >= ^(1-- > 


1 

'2 


1 ,11 


11 
'22 


>), 


(68) 


= 00 >= ^ — > 
^ ' ' y2^'22 


1 

'2 


1 ,11 

-2>-l2-2> 


11 
'22 


>)■ 


(69) 


Eq. (p^) can be written in the following way 












(pJLSgiko, \k\) = S^iko, 




s)gjLSsi.h, |fc|), 






(70) 


where Sg{ko, \k\; s) is 














-^(- 


^-^o-%)-\ 






(71) 


S- = {^ + ko + E}^) 


-\- 


^-fco + i?fc)-\ 









Se = {\- kl - El)iC- - kl - Elf - Ak:iEl)-\ 
So = {-2koEf.){C- - - Elf - AklEl)-\ 

2.6.2 Matrix Representation 

In the matrix representation, we replace the spinor of a second particle by the transposed spinor and 
then calculate the direct product, 

U'slik) ® U:^{-k) f/;;(fc) ® U^fi-k). (72) 

We use matrices (4 x 1) (g) (1 x 4) = (4 x 4) instead of matrices (1 x 4) (g) (1 x 4) = (1 x 16) in spinor 
space. Then matrices (4 x 4) can be expanded by Dirac 7-matrices. 
The BS amplitude in the rest frame can be written as: 

$^^(P, k) = x'^iP, k)Uc = J2 ^JLSoih, \k\) T'^'^'^ik) Uc, (73) 

LSg 

where Uc is the charge-conjugate matrix 

Uc = «727o, (74) 
and the spin-angular part Y'^^^^'^{k) has a structure similar to 3^'^^'^^'-^(fc), but with the replacement 



We show here, as an example, the function r^'^^'-^(fc), with notation ^'^+^Lj for partial states 
{JLSg): 



SlS2 

V / 

1 + 7o / -cr ■ ^ -ia2 



(m + pi • 7) iM-l{m-p2-l)Uc . 



Here we make use of the relations for Pauli spinors and Clebsch-Gordan coefficients: 

V2j2iW2^2\lM)xsal = (o- ■ C^) (ta^). (75) 

£^ is a 3-vector of the polarization of a particle with spin one and the components in the rest frame, 

e+i = (-l,-^,0)/x/2, ^_, = {l,-t,0)/V2, |o = (0,0,1). (76) 

The polarization 4- vector = (0,^^) is determined in the rest frame. 

In the general case, we can separate p-dependence and rewrite the spin-angular functions r'^^^^''^(k) = 
pjL5pip2,^(fc) as: 

.jLS—,Mtu\ _ P2-l-m -1 +7of^jL5,^/ X pi-'y + m 



^2E^{m + Ef^) 2 V2^fc(m + £;fc)' 

^2Ef^{m + Ef^) 2 V2^fc(^ + ^fe)' 



'2Ej^{m + E^) 2 ' ^2Ej^{m + Ej^y 

The conjugate functions can be written in the following form: 

f''^^^'^(fc) = 70 [r^^^^'^(p)]So, (78) 
Then the orthogonalization condition can be presented as 

j dQf^ TT{r'^^^'^\k)r''^'''^''^'{k)} = Smm'SluSss'S,,'. (79) 

Using the Pauli principle for identical particles 

$-^^(P, A;) = -Pi2$-^-^(P, k), (80) 

where P12 is the permutation operator of two particles, we can write the BS amplitude X'^'^ {P-i k) in the 
rest frame as: 

X^^(P,A;) = {-iy+'Uc[x'^{P,-k)fUc , (81) 

where / is the isospin of the system. These relations give us the symmetry properties of the radial 

functions (pjiSeiko, \k\) for the transformations, when we replace ko —ko. Radial functions for 
different LSg will be odd or even under this transformation. In order to have a radial function with 



-""So— Channel. The BS amplitude of the two-nucleon system in the ^^o-channel has four states: ^Sq , 
^Sq , ^Pq, ^Pq (or '^-Pq^ ' "^-^0^^)- For A^A^-scattering we take ^/s = \/P^. The corresponding spin-angular 
parts are shown in table where pi = {Ef,, k) and p2 = {Ej^, —k). 

Table 1: Spin-angular parts f"'^'^^'"^ for ^^o-states 



i^rHpi ■ 7 -P2 ■ 7)75 



Using the Lorentz invariant expressions for ko and Ej^: 



we can rewrite the expression for the BS amplitude in a covariant form. It is, of course, more convenient 
to use the direct covariant form of the BS amplitude in the ^5*0 channel, written in the 4x4 matrix 
form. For this purpose, we introduce four Lorentz covariant functions, hi{P ■ k, k"^): 

v/4^X°°(^,^) = /ii75 + /i2 - (Pi ■ 775 + 75P2 ■ 7) + (83) 

m 

, ,Pi - 1 - p2- 'y + m. , pi-7-m p2- 'f + m 

hs ( 75 - 75 ) + 75 , 

m m mm 

where 4-momentum of the two particles, pi and p2, are determined by (|47|) , and functions hi{P ■ k, k"^) 
can be expressed via the radial functions 4>jLSe{ko, |k|) defined in Eq.(|7^): 



h = 72/^1(0, ++0. J + 72/4(0^ ,-0. J- /iA;o|ferV3,.-8m|fcriDo03po, 
h2 = im|fc|-V3,e, (84) 
hs = 8aom2(0 +0 )-|/iA;o|/^^V3pe-8ao"^|fc^^e("^-%)..p 

,3|K|-1 



-^0 ^0 



hi = -4ao\/2m^(0^ ^ + 0^ _) + 8aom^|fc 



Here 



ao = l/{l6MEf^), 6 = 2m + Ef^, fx = m/M 

Do = ao(4A;2 + 16m2 -4E| -M^), (85) 
Di = ao(-MV4 + A;2-3E| + 4m2). 

We note that the functions /i2 and 03pe are odd under ko —ko, while other functions are even. 

^S'l — ^ Z?i— Channel (Deuteron). The BS amplitude for the deuteron^ has eight states: ^S^, 
^S^,^Dt, ^D^, 3pe^ 3po^ ipe^ ipo^ 3 p+- ^ 3p-+^ 1 p+- ^ ^Pf^), which are numbered as 1,...,8. 
The corresponding spin-ang ular parts r^^^'^{k) are tabulated in table ||. 
The BS amplitude has the following covariant matrix form, 

V4^x''^{P,k) = h,iM-l + h2^+ (86) 

m 

\ m m J 



Table 2: Spin-angular parts f ^^'^'•^ for the deuteron 



ILS 






3 C 




-7 

7+i(Pi-7-P2-7)(^ex)|fc|"'] 


'Pi 
'Pi 


Y 2 


• 7(Pi • 7 - P2 • 7) - (k^u)] \k\~' 
VS{k^M)\k\-' 



\ m J m 

/^ Pi-7-m p2-7 + m \ 

\ m m J 

^ / pi ■ 7 - P2 ■ 7 - 2m \ k^M ^ 
^ \ m J m 

m \ m J m 

Here we have introduced eight covariant functions hi{P ■ k, k"^), which are connected with the radial 
functions (pjLSgiko, \k\) in the rest frame via relations: 

hi = D+ (03^^ - V2</.3^J + L»r (03^_ - V2</.3^_) + lV6i^ko\k\-' 03pe + V6mDo\k\-' 03^„, 

h = V2{D^ <f>^^^ + Dt cl>^^_) - D+ <f>^^^ - <l>^^_ (87) 



^3 = -jV6m\k\ Vspe, 

/i4 = 8aiV2mko (^3^+ - ^3^, ) + 8a2£ m/co(03^+ - 03^- ) 

-16aoA/3m^|fc|"^ (/^o^ipe - %0ipJ> 
/i5 = 16aom2 [03^_ + 03^^ - V2 (03^^^ + 03^_ )] + 8aoV6m\k\-^ %E^. cj>,^^ + 

(2m2-4)03^J, ' 

16aoV^m2|fc|-i (m03 - J, 

1 1 

/i7 = 4aom2 [V2(03^^^ + 03^_) - (^3^^ + (t>^^_)] - Aao^^^w?\k\-^ cj>,^^, 
hs = 4aom^|fe|-2[v^(m-£;fc)(03^^+03^_)-(2£;fc + m)(03^^+03^_)] + 



The coefficients in these expressions are: 



oi = aQm/{m + Ej^), = ai/ {m - Ej^), 

Df = ao(4A;^ + 16m2-M^-4£;|±4M£;fc), 

= ai{16m'^ + 16mEf^+4:El, + -4:kl±4:Me), 

Df = a2[-12mEl^ + 2M^Ef^-8k^Ef^ + 16m' + mM^ -Amkl + SEl, 



4- /'I fi ^ A/f _l_ /I 



Df = ao(-2E| +4m2 + 4m% ±eM), 
Di = ao(2£±M), 

and ao, e, /i, are given in Eq. (^51). Here functions h^, /i4 and 03 ^, 0^ ^ are odd and others are even 
under ko —ko. 



2.7 Construction of the Light-Front Wave Function from the BS Am- 
plitude 



Here, we consider the relation between the BS approach [22| and the hght front dynamics (LFD) 
approach for a two-nucleon system 0, which is one of possible 3D relativistic dynamics proposed by 
Dirac |5^. In this approach the state vector describing the system is expanded in Fock components 



defined on a hypersphere in the four-dimensional space-time. The LFD approach is intuitively appealing 
since it is formally close to the nonrelativistic description in terms of the Hamiltonian, and state vectors 
can be directly interpreted as wave functions. 

The equivalence between these two approaches has been a subject of recent discussions presented in 
ref. 0, |5^ and references therein. The application of the two approaches for a two-nucleon system 
(deuteron) serves to clarify the structure of the different components of the amplitude. It is also useful 
in the context of the three-particle dynamics, where the proper covariant and/or light-front construction 
of the nucleon amplitude in terms of three valence quarks (including spin dependence and configuration 



mixing) is presently discussed [|59|, Although the relation between the light-front and the Bethe- 
Salpeter amplitudes for a two-nucleon system has been spelled out to some extend in a report 0, we 
provide here some useful details. 

To proceed, we present different ways of construction of a complete (covariant) Bethe-Salpeter 
amplitude, see, e.g. [0, Beside the spin structure of the wave function (amplitude), we also present 



a comparison of the radial pa.Tt of the amplitude on the basis of the Nakanishi integral representation [43 



The representation is well known and elaborated for the scalar case, see, e.g. [H, 53]. However, this 
formulation was not followed up later in the actual applications. We employ this integral representation 
in order to establish a connection between the two different approaches. 

In this context, the direct product representation used in the rest frame of the two-nucleon system 
using the p-spin notation is close to the nonrelativistic coupling scheme and provides states of definite 
angular momentum. To construct the covariant basis, this form is transformed into a matrix repre- 
sentation which will then be expressed in terms of the Dirac matrices. A generalization to arbitrary 
deuteron momenta finally leads to the covariant representation of the Bethe-Salpeter amplitude. This 
was explained in the previous section along with an explicit construction of the deuteron (J = 1) and 
the J = two-nucleon state. 

We now compare the BS amplitude to the covariant LFD form. The state vector defining the light- 
front plane is denoted by u, where u = (1, 0, 0, —1) leads to the standard light-front formulation defined 
in the frame t + z = 0. We start from the integral that restricts the variation of the arguments of the 
Bethe-Salpeter function to those of the light-front plane. 



/ = y d^xi d^X2 5{u ■ xi) 5{ul! ■ X2) $^(a;i, X2) exp(iZi ■ xi + ih ■ X2) , (88) 
where h, I2 are the on-shell momenta: if = I2 = 1 and 

$^(xi,a;2) = <l>^='(""-'*''^)(a;i,a;2) =< Q\T{m{xr)^!{x2))\P > (89) 

is the Bethe-Salpeter amplitude (|rT]). 

To go further we represent first the S -functions in (|8^) by the integral form, 

5(uj ■ Xi) = — [ exp(—iuj ■ Xiaijdai, 5(uj ■ X2) = — f expi—icu ■ X2a2)da2 , (90) 



We introduce the Fourier transform of the Bethe-Salpeter function $(P, k) as the Fourier transform of 
<l>^ defined as, 



X2) = {2n)-^/^ exp [-iP ■ (xi + X2)/2] l'^(x) , x = ^^—^ . (91) 

The BS function is 

*,P,.,^/*-Wexp(ft..M*., (92, 

where k = (pi — p2)/2, P = pi + P2, Pi and p2 are usual off-mass shell four- vectors. Making the change 
of variables ai + 0:2 = t, (^2 — cn)/2 = /?, we obtain: 

Pi = /i — ljt/2 + 
P2 = h- ujt/2 - uj(3, 

and the integral (^) takes following form: 

/+00 j-+oo / J 7 \ 

5^^^ (/i + /2 - P - t^r)rfr y $ U + /2 - cc^r, + cj/?J d(] . (93) 

On the other hand, the integral (^8|) can be expressed in terms of the two-body light-front wave 
function. We assume that the light-front plane is the limit of a space-like plane, therefore the operators 
\l/(xi) and \l/(x2) commute, and, hence, the symbol of the T-product in (|89|) can be omitted. In the 
considered representation, the Heisenberg operators \l/(x) in (^) are identical on the light front u-x = 
to the Schrodinger operators (just as in the ordinary formulation of field theory the Heisenberg and 
Schrodinger operators are identical for t = 0). The Schrodinger operator \l/(x) (for the spinless case for 
simplicity), which for u ■ x = is the free field operator, is given in p[: 

^(x) = - — \t7k [ ^(k)exp(ik ■ x)d'^k 

{27iy/^ J 

1 f d^k 
= ^2n)^/^ J ^^''^'^ exp(-ifc ■ x) + a^{k) exp{ik ■ x)] -j= . (94) 



We represent the state vector \P) in (|89D as the Fock components of the state vector. Since the vacuum 
state on the hght front is always "bare", the creation operator, applied to the vacuum state (0| gives 
zero, and in the operators \E'(x) the part containing the annihilation operators only survives. This cuts 
out the two-body Fock component in the state vector. Thus the integral I can be obtained in the 
following form B: 



(27r)3/2(cj . P) 
2{uj ■ h){uj ■ I2) 7_ 



/+00 
ij{h,l2.P,uoT)6^^\h + l2-P -uot) dr , (95) 
-00 



where ip{li,l2, P,ujt) is two body wave function defined on the light front specified by u. We make a 
comparison (|93|) and (|95|) and find the formal connection between the LFD wave function and the BS 
amplitude: 

^(/i,/.,P,^r)= r<t>(p,^-l^ + Pu)df3, (96) 

7r{uj-P) V 2 J 

Since li and I2 are on the mass shell, it is possible to use the Dirac equation after making the 
replacement of the arguments indicated in Eq. (^). This will be done explicitly for the J = and the 
deuteron channel in the following. 

2.7.1 ^Sq - Channel 

Using the Dirac equations 'u(/i)(7 ■ /i — m) = 0, and (7 ■ I2 + m)Cu{l2Y = O5 o'^^ obtains the light-front 
wave function from the Bethe-Salpeter amplitude using Eq. (|73D , (pBf ) 



The functions Hi{s, x) and H2{s, x), depending now on x = uj ■ li/u ■ P and s = (/i + ^2)^ = 4(g^ + m^), 
are obtained from the functions hi{P ■ k, fc^) through the remaining integrals over (3 imphed in Eq. (p^), 

Hf'is,^) - iv/M(l-2.„.-MV/5..P.-V4.™^.(2.-l)^)..P.^ 

= N j hi{s,x,P')dp', 

Ht\s,x) = N [k{s,x,P'){PTdP' (98) 



where the variable /?' = /3a; • P has been introduced, and N = x{l — x), 1 — x = cu ■ I2/U! ■ P. We would 
like to stress that in this case the functions and do not contribute. Instead of the four structures 
appearing in the Bethe-Salpeter wave function, the light front function consists of only two. Note, that 
the second term in the parenthesis is defined by the pure relativistic component of the Bethe-Salpeter 
amplitude. 

2.7.2 ^Si-^Di - Channel (Deuteron) 

In the deuteron case, starting from the formula Eq. (P^D, replacing the momenta Pi, and applying the 
Dirac equation we, arrive at 

m muj ■ P 

I Off (1) ^ ' ' ^ + 2^(1) 7 • eM7 ■ uj - 7 ■ ujj ■ eM 



+ |2/ff + 2fff']^:,^. (99) 



m'^uj ■ P CO ■ P 

r(2) , n u(2)i ^ • gM7 • ^ 

where if^^'^'' are defined in Eq. (pSl) . In this case the functions /i4 and hs do not contribute. The 
expression (7 ■ eM7 • a; — 7 ■ co'7 ■ eA/ ) at the term given in Eq. (^) can be transformed to a different 
one to compare it directly to the light-front form. Using in addition the on-shellness of the momenta h 
and I2, the resulting form is 

■ eul ■ uj - ■ u-f ■ eM)Cu2 = -^Uil-i'^^e^^p^ie^hyhpUJ^ 

+k ■ eM ^ ■ P — ■ eM ^ ■ +P 
_l(s_M2)(x-iv-eM 

+ \m{s - M^) :L^^^^]CuI (100) 

The final form of the light-front wave function is then 

^, k ■ eu , jj, ^ ■ Cm , jj, k - eul - ^ 

m 



XiM H[-i-eM + H'^ - + ^3 ^ + ^4 2 w 

— muo ■ P m^LJ ■ P 



+H^i'j5e^^p^e^lij2pUJa + Hq^- — , (101) 



with the functions 







H'2 


= iff + l2i7«. 


Hs 










= -2Hi'\ 




ms 


m 





[s - M^] 



{2x - l)2Hi'\ 



2t:iJ^ni^2HP . ri02^ 



Provided the invariant functions hi are given from a solution of the BS equation, the above relations 
allow us to directly calculate the corresponding light-front components of the wave functions. 

Thus, the projection of the Bethe-Salpeter amplitude to the light front amplitude reduces the number 
of independent functions from eight to six for the — Di channel, and from four to two for the ^Sq 
channel. The reduction takes place because the nucleon momenta pi and p2 are on mass shell in the 
LF formalism. The result is based on the application of the Dirac equation and the use of the covariant 
form. Any other representations (e.g. spin-angular momentum basis) also lead to a reduction of the 
number of amplitudes for the two-nucleon wave function that is however less transparent. For an early 



consideration compare, e.g. ref. [33 



2.7.3 Integral Representation Method 

A deeper insight into the relation between the Bethe-Salpeter amplitude and the light front wave 



function will be provided within the integral representation proposed by Nakanishi 0, This 
method has recently been fruitfully applied to solve the Bethe-Salpeter equation both in the ladder 
approximation and beyond within the scalar theories [HI, In this framework the following ansatz 



for radial Bethe-Salpeter amplitudes (^) of orbital momentum L has been proposed, 

<PjLS,{ko, \k\) ^ <t^j,sAP ■ k, e) = da dz ^^^,.l'"jiTp.k_^^y (103) 

where gjLSe{c(, z) are the densities or weight functions, k = m? — M2/4, and the integer n > 2. The 
weight functions gjLSgioi, z) that are continuous in a vanish at the boundary points z = ±1. The form 
of Eq. ( p.03|) opens the possibility to find the Bethe-Salpeter amplitude in the whole Minkowski space 
while commonly used solutions are restricted to the Euclidean space only. In fact, the densities could 
be considered as the main object of the Bethe-Salpeter theory, because knowing them allows one to 
calculate all relevant amplitudes. 

For the realistic deuteron we need to expand the Nakanishi form to the spinor case, which has not 
been done so far. The key point in this procedure is a proper choice of the spin-angular momentum 
functions as well as the integration over angles in the Bethe-Salpeter equation. The choice of the 
covariant form of the amplitude allows us to establish a system of equations for the densities gij{oi, z), 
suggesting the following general form for the radial functions hi{P ■ k, k"^) (even in P ■ k) 



r"00 + 

hi(P-k,k'^) = I da dz 







P- A;)" 

(104) 



1 {{a + k"^ - k"^ - zP ■ k)"" 

9i2ia,z) P 

(a + /s:2-A;2-zP. A;)"+i 

gisja, z) {P ■ kf 
{a + - k^ - zP ■ k)^+^ 

For the functions that are odd in P ■ k the whole integrand is multiplied by factor P ■ k. Although now 
the number of densities is larger, the total number of independent functions is still eight. The form 
given in Eq. ( |105|) is valid only for the deuteron case. The continuum amplitudes of the ^5*0 state, e.g., 
require a different form. 

It is a major advantage that we can perform the integration over (3' in the expressions of Eq. ( pH]) 
by using the integral representation. Substituting the arguments of the functions hi into the integral 
representation Eq. (|105D leads to a denominator of the form 



V''{a, z; X, s, /?') = (« + |(1 + (2a; + 1)^) - (3'{2x - 1 + ^) - ie)'' (105) 
Using the identity for an analytic function F{z) 

f^^ , Z""" F(z) m Fil -2x) , , 



we can express the radial amplitudes in terms of Nakanishi densities. Thus, H^{s, x) reads 



x{l-x) f r^5i(a,l-2x) 052(0, l-2x)sx(l 

H5{S,X) = / da < H -. -; r-TT > (107) 

^ ' s J {{a + sx{l-x)) {a + sx{l-x)y J ^ ^ 

Note, that the dependence of the amplitude on the light front argument x is fully determined by the 
dependence of the density on the variable z = 1 — 2x, which has also been noted in ref. |0 for the 
Wick-Cutkosky model. 

This fully completes the connection between the Bethe-Salpeter amplitude and the light front form. 
The evaluation of the Nakanishi integrals does not lead to cancellations of functions. Although some 
functions are cancelled for reasons given above, all spin-angular momentum functions (or all densities) 
in principle contribute to the light-front wave functions. 

Once the Bethe-Salpeter amplitudes are given (or the Nakanishi densities) the light front wave 
function can be calculated explicitly. The Nakanishi spectral densities of the Bethe-Salpeter amplitudes 
lead directly to the light-front wave function. 

We would like to stress that the two relativistic approaches have shown qualitatively similar results 
in the description of the electro-disintegration near the threshold (see chapter (D). The functions /s 
and g2 (notation of ref. [0]) may be related to the pair current in the light-front approach whereas the 
functions and h2 play this role in the Bethe-Salpeter approach. The results presented here allow us 
to specify this relation on a more fundamental level. 

2.8 Solution of the BS Equation 

In the previous section we have considered the general properties of the BS amplitude. The different 
representations for spinor and angular parts were proposed. The most important part defining the 
dynamical structure of the bound state is the radial part which can be found by solving the BS equation. 

In this section, we consider the solution of the BS equation with a separable interaction (separable 
ansatz). The separable interaction is well known in nonrelativistic approaches. In particular, this 
suggestion allows one to reduce the system of integral equations, in which the Lippmann-Schwinger 
equation is transformed, after partial expansion, to a system of linear equations ^ The Bethe- 



Salpeter equation with separable kernel can be transformed to a system of linear equations as well. 
2.8.1 Separable Interaction 

Let us consider the BS equation for T- matrix after partial decomposition (see subsection p.6| ): 

'^jL's',',jLsM^ l^'l' ^0, \k\;s) = V^,^,5,,,,,^s,(A;o, \k'\,ko, \k\; s) + (108) 

i 



/dkQ / k d\k j^//g//„// (^Q) |fe I7 ^Q) s) Sg'i{kQ, \k"\; s) Tj^„^„^„ jj^^^ik^, \k"\^ kg, \k\; s) 

r ft c// _// 



L"S"g" 



where Sg{ko, \k\; s) is given by Eq. ([7T]). 

The separable ansatz for the interaction is introduced in the following manner [^, |49 

N 

V.ws',,j.sSK, \klko, \k\;s) = J] X,,gf''<^\ki |fc'|) ^^(^o, |fc|). A,, = A,„ (109) 

where is a rank of separability, Xij are parameters of the interaction kernel, and gf^^^ko, \k\) are 
functions which define the interaction. Then according to Eq. ( |108|) the T-matrix can be expressed in a 
separable form too. We assume the following separable form for it 

N 

Tj,,s',,jLs^ko,\k'\,ko,\k\;s) = ^ n,{s)g^'^''''\ki\k'\)g';'^''{ko,\k\). (110) 



Substituting Eqs. (|1Q9D and {^) in Eq. ([lOSD we can obtain an expression for r(s), 



{r-\s)),, = {X-%-H,,{s), (111) 
where Hij{s) is determined by equation: 

The solution for the radial part of the BS amplitude can be presented in the following form, 

N 

(pjLSs(^o,\k\) = Sf,{ko,\k\; s) Xij g-^^%ko,\k\) Cj{s), (113) 
where the coefficients Cj{s) satisfy the following system of linear homogeneous equations: 

AT 

Ci{s) - J2 Hik{s)XkjCj{s) = 0. (114) 
k,j=i 

2.8.2 Dispersion Analysis of Nucleon-Nucleon T -matrix 

To analyze the analytic properties of the solution for the T-matrix with separable interaction, let us 
consider a rank I case. We take only an S'-state in ^5*0- and ^Si — ^ Di-channels (namely ^Sq- and 
'^S'j^'-waves). Omitting all partial waves and separability indices one could write a solution for T-matrix, 

t{k',, \k'lk,, \k\; s) = T{s)g{k',, \k\)g{ko, \k'\) (115) 

with the function r(s), 

T{s) = l/{X-' + h{s)), (116) 

and function h{s), 



[gjko, \k\ 



Here we have introduced small letters for the functions T and H in the simple indexless case. 
As a result, the T-matrix can be rewritten in the following form: 

«(,;jfe'|,,,j*|,,) = £M4M^, („8) 



X-^ + h{s) 



and the on-mass-shell expression is: 



t(s) = ^= f^^Q'^'^]' (119) 

d{s) x-^ + h{sy ^ 



It should be noted that the Eq. ( |119D has the so-called N/D-ioim widely used in the nonrelativistic 



T-matrix theory ||l[ and some methods of relativization of the theory. 

We use the following representation for the on-mass-shell T-matrix valid in the region of unitarity: 

t{s)^t{0,p,0,p;s) = -^E^e^^(^) sin5(.), (120) 



Using Eq. ([TT9|) , it is easy to relate the T-matrix and the phase shift S{s). To achieve this, we 
assume the imaginary part of the function n{s) satisfies the following condition: 



Imn(s) = 0. 



:i2ii 



The condition is related with the specific choice of (^-functions for the A^A^- vertex which will be discussed 
later. Taking into account Eqs. (|119|) and (|121|) , the phase shift 6{s) can be given as 



cot 6{s) 



Ret{s) 
Imt(s) 



X-^ + Reh{s) 
Im h{s) 



(122) 



To express the low-energy parameters in term of the T-matrix solution, it is suitable to expand the 
function h{s) in a series of p terms: 



h{s) = ho + iphi + p^h2 + ip^h + C'(p^), 
Reh{s) =ho+p^h2 + 0{p^), 
Im h{s) = p{hi + p^h + 0{p^)). 



(123) 
(124) 
(125) 



The low-energy parameters of A^A^-scattering are introduced by expanding the T-matrix into series of 

(126) 



p-terms following the expression suggested in ref. ||66|| : 

pcot 6{s) = 



-1 



and taking into consideration only the first two terms of the decomposition ( |126| ). 

Using now the definition ( |126| ) and Eqs. (|122|) - (|125|) one can find the parameters Oq and tq: 



ao 



A-i + ho 
2 



hi 



+ ho)-^ - h2 
til 



(127) 
(128) 



At the mass of the bound state squared M^, T-matrix has a simple pole at the total momentum 
squared s, and the bound state condition can be written in the following form: 



^(^0; 1^ I' ^0; 1^1) ■s) 



B{k'„\k'\,ko,\k\;s = M^ 



R{kQ, \k I, /cq, \k\] s). 



(129) 



where functions B and R are regular at the point s = M^. The bound state energy Eb is connected to 
Mf, as: Mfe = 2m — E^. The bound state condition ( |129|) , with the help of Eq. ( [118| ), can be presented 
in a form: 



A 



-1 



Ms = Ml). 



(130) 



Let us consider now the analytic properties of the solution (namely, function h). The simplest choice 
of the function (7(A;o, \ k\) is the Yamaguchi type function |46| , ^ : 

(131) 



g{koAk\) = {kl-k''-l3^ + ie)-\ 
In this case, the function h can be rewritten as follows: 

1 



h{s,(3) 



47^3 



dkn 



k'd\k\ 



(v/i/2 -E^ + lef - kl kl -EI + 



le 



(132) 



where 5^2 = djdfp' and E^ = ^/^~+~^. We introduced the second argument (3 to underline the 



ovrvl 1 r>i + r\ o-rv on r\ on r-o r\-f Vi o f n n r>+ ir^n Pi r^n ■fViic! t^ot*o tyi o^f ot* 



Analyzing Eq. (p2D one can identify four poles in the complex plane ko, namely: 

ki^\s) = ^-Ef^ + ^e k^^\s) = -^ + E^-te (133) 

kf^ {s) = -Ep + le fcj') {s) = Ep-ie 

The variation of s results in the move of the poles fcg^'* and /cq^"*, and one confronts the situation 
where two poles "pinch" the real fco axis. It means the function h{s) has in this s-point the leap and 
imaginary part. First points in which this condition is satisfied (branch points) can be found from the 
following equations: 

fc«(s) = kf\s) So = 4m2, (134) 

A;«(s) = A;i^)(s) si=4(m + /?)l (135) 

Summarizing the situation, one could say that the function h{s) has two cuts starting in points sq 
and Si respectively, and therefore can be written in a dispersion form: 

(136) 

J s' — s — le 

p{s', P) = 9{t - Am^)p,i{s', P) + e{t - 4(m + PY)pUs', (3) 

with two spectral functions Pei;in (e/ stands for elastic and in — for inelastic) which are connected with 
the imaginary parts as follows: 

p{s\P) = -\mh{s\P) = ^.{h-h*). (137) 

Let us note some analytic properties of the obtained solution for T-matrix (see also subsection |2.5| ): 

1. A bound state (deuteron) appears in the T-matrix as a simple pole in total momentum squared 
s at the mass of the bound state squared, s = M^; 

2. in the region s > 4m^, T-matrix has the cut corresponding to the elastic NN scattering (Eq. 
O); 

3. in the region s > 4(m + /3)^, T-matrix has the cut corresponding to the inelastic A^A^ scattering 
(Eq. dm)); 

4. T-matrix has no left-hand cuts but there is a pole of second order at the point s = 4(m^ — jS"^). 
To find spectral functions one should perform fco-integration in Eq. ( |132| ) which results in: 

1 

27^2 -I'' J ^I'-i s/A - y/^E}^ + m^-(3^ + ie ~ ~ ~ ^ ~ 
Taking into account the following symbolic equation: 



v/i - 2% + ie 2Ef3 



138) 



1 V 

T i7i6{x - xq) (139) 



X — xo±ie X — Xq 
it is easy to find spectral functions: 

s' > 4m^, 

Pei{s', (3) = V^Vs' -4:my{7c\s' - 4m2 + Af3^f), (140) 
s' > 4(m + /3)\ 

Pin{s', P) = -{64(3^ + 16/3^s' - 192/5^m2 - 20ph'^ 
+192/3^772^ - 32p^s'm^ + 16m^s' + 3s'^ - 64m^ - 12m^s'^)/ 



To perform integration in Eq. ( |132| ) it is suitable to introduce new variables: 



m 



P 



1 



4m^ 



s t 

2 m + (3 1 + Afl 

V = = -, 

m- j3 1 - A/3 

If the following conditions are valid: 



Am? ' 



w 



a 



m?-(3^ 1-A/' 
s-A{m + (3f _ t-{l + \pf 



s-4(m-/3)2 t-{l-\pf' 



(142) 
(143) 

(144) 



m > /3 > =^ 
4(m + I3Y > s> Am? 



1 > A« > 0, 



(145) 
(146) 



the parameters are real and positive: 



1 > p2 > 0, 
i;^ > 1 > 0, 



> 0, 

> 0. 



(147) 



The condition ( p,46| ) means that the second (inelastic) imaginary part does not contribute to the function 
lm.h{s) when phase shifts are calculated in the region 4(m + /3)^ > s > Am? and, therefore: 



2\2 



Im/i(s,/5) = \v!ihf,i{s^ [3) 



Am'^n{p'^ + w'^y 



if 4(m + PY > s > Am^ 



(148) 



Performing integration (|136|) one can obtain the real part of the function h{s) (imaginary is given by 

Eq. (ina)): 



Re h{s, (3) = hei{s, (3) + /ii„(s, (3), 
hei{s,/3) = . „ „ ^In 



1-p 



l + P 



(149) 

(1 - p^)(l + W^)\W^ - p^) 1 (l-p2)(l+«;2) 

arctan h 



2\2 



32m%2(y4^^2) 
(l + a2)(l + i;2)2 



In 



+ 1 



16m'^TT'^{v'^ — l)(f^ — cr^) 



16m^7r^(f^ — l)2f (i;^ — o"^)^ 
(l + a2)(l + t;2)3 



16m%V(t;4 + (j2)(i;2 _ ^2)2^2(^2 _ 1)2 



v"" + v^cr^ - 4f ** + Av^a^ - v- a^) arctan - 

V 



[-v^ + 2t;V^ - 3t;V2 + 3(t%2 - 2(t^ + a^) arctan 



0" 



To find also the expressions for the low-energy parameters we should return to Eqs. ( |125D -( P!^ ) and 
expand function h{s) in a series of p terms: 



ho{l3) 



(1 + ^2) 2N . In 

- — „ „ „ (w + [l + w ) arctan — ) 
Am'^Ti'^w^ w 

5 + 20t;2 + 5t;8 + I4t;4 + 20t;6 l3v^ + 1 , /3t;2 + 1 

H : — ^ ^, r, , ^ — \ I — ^ arctan 1 / 



(150) 



4m27r2(i;2 _ 1)3(3^,2 + 1) y y"^ + 3 
(v^ + 1)^ 11 

arctan 



h2m 



Am'^7r^v{v'^ — 1)^ 



In 



- 1 



4m27r2(ti2 - ly 



, ^ ^-(3w + (3 + w'^) arctan — ) 



(151) 



-(304t;2 + 1212i;^2 + 3790w^ 



4m%2(^2 _ 1)5(^2 + 3) (3^,2 + 1)2 
_LQn/1^,14 I 070/1-,, 10 _1_ 0Q„,16 I 1 01 0.1,4 , 07n/1.,,6 I nQ\ 



3^2 + 1 l3v^ + 1 

arctan 1 / 1 



^2 + 3 



f2 + 3 



(v^ + lOv^ + l)(v^ + 1)^ 1 

H — - — — arctan - 

Am^Ti^viv^ — 1) f 



+ 



In 



1 



1 



+ 



llv^ + 52t;6 + 66f^ + 52f2 + 11) (t.^ + 1) 



/^3(/3) 



2\2 



(l + W^) 

(1 + ^2)2(4 + 3^2) 



(152) 
(153) 



At this moment, we can find internal parameters of tlie separable interaction (A, (3) to reproduce 



experimental values for low-energy parameters a| 



exp 
Os 



-23.748 ± 0.010 fm, rgf 



2.75 ± 0.05 fm 



for singlet channel (^5'o), and Oq^^ = 5.424 ± 0.004 fm and bound state (deuteron) energy 
2.224644 ± 0.000046 MeV for triplet channel (^^i). Experimental data are taken from ||67|| . 
In the case of -"^So-channel we use Eq. (|128|) with = a^^f to find A: 



exp 



Inserting the above expression into Eq. (|128|) with tq = Tq^^ we find: 

2 



exp 



(154) 



(155) 



Solving nonlinear Eq. ( |155| ) we find the value of (3, and then using Eq. ( |154| ) — the value of A. 



In the case of Si-channel we obtain A from the bound state condition ( |130|) with Ei, = E\ 

\~' = -h{s = {MT')\(3l 



7\exp_ 



where Md = 2m — Ed- Inserting the above expression into Eq. (128) with qq = Ciot^ we find: 



exp 



hoi/3) -his = iMry,/3y 



(156) 



(157) 



Solving the nonlinear Eq. ( [157| ) we can find (3, and then using Eq. (|156|) — the value of A. As a result 
we find: 



for ^So channel: A 
for ^Si channel: A 



-0.29425404 GeV"^ 
-0.79271213 GeV-^ 



/3 = 0.22412880 GeV, 
P = 0.27160579 GeV, 



;i58) 



The phase shifts 5s (s) and 6tis) calculated with these parameters are displayed in Fig. |1|. Experi- 
mental data are taken from ref. ||68|. As it is seen from the figure, the simplest choice of the separable 
interaction — rank I with just two parameters A and (3 — is able to provide the low-energy parameters 
of elastic NN scattering a<j and in singlet channel, and at and bound state (deuteron) energy Ed in 
triplet channel, with required accuracy and to reproduce phase shifts up to Tiab — 100 MeV. 

We note in conclusion that the dispersion form of the T-matrix for the elastic A^A^ scattering obtained 
for the separable interaction allows to perform analytical calculations and explicitly connect parameters 
of the kernel and the observables [^. It is interesting to note that the rank I separable interaction is 
able to provide the low energy data even up to Tiab — 100 MeV and the deuteron properties. 



2.8.3 Covariant Graz-II Interaction 



In the actual calculations for various electromagnetic observables, we use a more complex separable 




T,ab(MeV) 



Figure 1: The triplet 5t (top panel) and singlet 5s (bottom panel) phase shifts calculated with the 
separable interaction in Yamaguchi form with rank I. 



are taken into account ^ ^Df). In this case the functions Qi have the following form 



1*1) = ^^-#^. (159) 

ifci)-- <'-°-'''> 



(/Cg k I3I2Y 



(/Cq fc /^2l)(^0 ^ /^22)^ 

r^i'^M^o, |fc|) = ^^a'^M^o, |fc|) = 9?hko, \k\) ^ 0. 

The parameters of these functions are given in table ^. 

The solution of the BS equation for the vertex function can be explicitly written as 

3 £) + 

g-ig+{ko,\k\) = (ciAii + C2A12 + C3Ai3)5fi ^ {ko,\k\) + (160) 



35+ 

(CiAi2 + C2A22 + C3A23)fl'2 ' (^0, 

g3jj+{ko,\k\) = (ciAi3 + C2A23 + C3A33)5(3 ' (/co, |fc|), 

where we take into account that the matrix A is symmetric. Vertex functions (73^+ {k^, \k\) and g3jj+ {k^, \k\) 
in Euclidean space {k^ = iko) are shown in Fig. ^ and Fig. ^, respectively. 

To calculate the phase shifts, we use the following parametrization for on-mass-shell T-matrix: 

, . 8 / cos2e e^^^s _ 1 ^sin2e c^^^^+^d) \ 

~ ~^/^s - 4m2) V ^sin2e e*(^s+5D) cos2e 6^*^° - 1 ^ ' 

where Ss {Sd ) are phase shifts of ^Si- {^D^-) waves and e is the mixing parameter. 

The calculated results are given in table ^ and in Fig. ^ The Bethe-Salpeter approach with the 
separable interaction provides the deuteron properties and also phase shifts of the nucleon-nucleon 



1 




Figure 2: Vertex function gas+i^^, \ k\). 



Table 3: Parameters of the covariant Graz-II separable interaction 



7i 


28.69550 








GeV-2 


All 


2.718930 


X 


10- 


-4 


GeV^ 


72 


64.9803 








GeV-2 


A12 


-7.16735 


X 


10- 


-2 


GeV^ 


Ai 


2.31384 


X 


10- 


-1 


GeV 


Al3 


-1.51744 


X 


10- 


-3 


GeV*' 


/3l2 


5.21705 


X 


10- 


-1 


GeV 


A22 


16.52393 








GeV^ 




7.94907 


X 


10- 


-1 


GeV 


A23 


0.28606 








GeV^ 


P22 


1.57512 


X 


10- 


-1 


GeV 


A33 


3.48589 


X 


10- 


-3 


GeV^ 



2.9 Separable and One-Meson Exchange Interaction 

In this section, we show the meson exchange interaction (meson-nucleon) in ladder approximation and 
its connection with the separable kernel with the Yamaguchi-type (7-functions. To achieve this, we 
introduce simple approximation for ladder kernel: 

V-(fc[„|fc^|;0,0)V-(0,0;fco,|fc|) 
y(0,0;0,0) 

Using the expressions of the kernel for scalar meson exchange (sc) [|^, |6l 

1 



V{ki\k'\;ko,\k\)^V{ki\k'\;ko,\k\) 



(162) 



V.cik'o, \k'\; ko, \k\) 



9sc 1 



Air 7r'^A\k'\\k\Ef^'Ef^ 



{E^,E}^ + m')Qo{z)-\k'\\km 



(163) 



where z 



(k 



12 



fi'^) / {2\k' \\k\) , /i is the meson mass and Qi{z) are the Legendre 



functions of second kind, Qo{z) = 2 -"^ In (2 + l)/{z — 1), Qi{z) = zQq{z) — 1. It can be shown that 

Vsc{k'o,\k'\;0,0)= lim Vscik'„\k'\;ko,\k\) = au>g{k'o,\k'\), (164) 

ko^Q,\k\^0 



where 




Figure 3: Vertex function gs^+ik^, \k\). 



Table 4: Properties of the deuteron and the low energy A^A^-scattering parameters in ^5*1 chanel with 
rank III Graz-II kernel 



Pd(%) 


eo 


Qd 




Pd/s 


tq (Fm) 


a (Fm) 




(MeV) 




(e/2m) 








NR Graz 11 4 


2.225 


0.2484 


0.8279 


0.02408 


1.7861 


5.4188 


Gov. Graz 11 4.82 


2.225 


0.2812 


0.8522 


0.0274 


1.78 


5.42 


Exp. 


2.2246 


0.286 


0.8574 


0.0263 


1.759 


5.424 



~g{k',,\k'\) = l/{k','-k''-f,'). (165) 

Expressions for Vsc{0,0; kQ,\k\) can be obtained from Eqs. (|164[ - p.65| ) by the following substitutions: 
k'o — >■ kf) and — > |A;|. To make the connection between parameters, we perform |fc'|/m-decomposition 
in the function a^^/ up to 0{k' /m^) term: 

9l 1 



Kc(0,0;0,0) = -|^i^i^. (166) 
4tt tc^ fi'^ 



Using the expression (|162|) , we find: 



2 

Vsc{k'„ \k'\- ko, \k\) = -^{^^fg{k'„ \k'\rg{ko, \k\). (167) 
By comparing this expression with the separable form of kernel introduced as 




Figure 4: (a) ^Si- and (b) '^Di-phase shifts calculated by using rank III Graz II kernel. Experimental 
data are taken from ref. 



we find the following relation between parameters: 



/3 = /i, 



(169) 



Equations ( |167|jl69| ) are valid also for vector-meson-exchange kernel {vc) with substitution gsc Qvc 
and /i is the vector-meson mass. The separable form and can be derived from the following expression, 



ale 1 



1 



47r vrMlfe'llfclE, 



-2{2Ef^,Ef^+m^)Qoiz) 



(170) 



For pseudoscalar-meson-exchange, the kernel of interaction has the form: 



Vps{kQ, \k \', kf), \k\) 



9sc 1 



1 



{Ef^'E}^-m^)Qo{z) + \k'\\k\Qi{z) 



in 7i'^4\k\\k'\Ef^>Ei^ 
and for |fc| ^ we can write: 

V,sik',, \k'\; ko, \k\) ~ k''e~gik',, \k'\rg{ko, \k\). 



(171) 



(172) 



In the latter case Vps tends to zero for |fc| — ^ 0, and it is impossible to find a relation between 
parameters similar to Eq. ( |169| ). 

To illustrate the relations between parameters we used Eq. (|169|) and calculated values corre- 
sponding to parameters fi and g'^/4:7c for scalar- and vector-exchange- mesons from ref. Results are 
given in table ^. 



3 Elastic Electron Deuteron Scattering 

The previous section was devoted to the analysis of the basic objects and methods of the BS formalism. 
It was shown how different dynamical processes involving bound states of particles can be included 
into the field-theoretical consideration. In this section, we consider the application of these methods 
to elastic electron scattering off the simplest bound system — the deuteron. We will highlight most 



Table 5: Connection between parameters of two kernels. 



li (GeV) g'/AT^ A (GeV^") 

0+ NN5 0.983 0.64 -.06265954891 

0+ NNa' 0.550 7.07 -.2166930823 

1- NNp 0.769 0.43 -.02576448048 

I- NNuj 0.7826 10.6 -.6577877886 



Our interest in the electron-deuteron scattering is connected first of all with recent experimental 
studies of the deuteron electromagnetic forms factors in the region of high transfer momenta (see, for 
example, |7^-[^), where the relativistic effects a priori play essential role, and with recent tensor 



polarization data ||73[|-||75[| 



Traditional nonrelativistic methods are based on the impulse approximation with allowance for 
relativistic corrections such as the meson-exchange currents (MEC) and retardation effects. In a number 
of investigations (see for example [|I^,|]TB[-|7^ and references therein) it has been shown that the correct 
account of these effects is necessary to explain the experimental results. Generally, the deuteron elastic 
form factors are known to be sensitive to the choice of the strong nucleon form factor and to the MEC 
models. On the other hand, the recent relativistic investigations |j79|-[0 show that some of the meson- 
exchange currents (in particular, the pair current) are automatically included in the relativistic impulse 
approximation. 

Beside the relativistic effects, it is important to study the contribution of the nucleon form factor to 
the deuteron elastic form factors and to its polarization properties. There exists a number of theoretical 
and phenomenological models of the nucleon form factors. One finds largest differences in the results 
of evaluation of the electric neutron form factor G|(g2) as well as of the ratio Ge/Gm for the proton at 

= —q^ > 1 (GeV/c)2. A further progress is related with the appropriate choice of the polarization 
observables, which would allow the consistent analysis of the structure of the bound nucleon. In this 
section we apply the BS approach to the analysis of elastic electron deuteron scattering including such 
topics as elastic form factors and polarization tensor of the deuteron. 

3.1 Relativistic Kinematics 

The differential cross section for unpolarized elastic electron-deuteron scattering in the one-photon- 
exchange approximation (Fig. ^) is expressed in terms of the Mott cross section and deuteron structure 
functions Alq"^) and B{q^) (the electron mass is neglected): 



da\ _ a^cos'^ej2 ^^^^^ 



dQ'^ J Mott 4^2 (1 + 2EJM sin^ 9^/2) ' 

where 6e is the electron scattering angle, M is the deuteron mass, is the incident electron energy, 
and 



8 0^2/ ''^ 2 

9 ' -Q^ 



A(g^)=W) + ;;r^^F^(g^) + ;r^F^(g^ 
4 



B{q') = -r]{l + r])Fl{q'), (175) 

where rj = —q^/AM"^ = Q^/AM^. The electric Fc(g^), the quadrupole Fqi^q"^) and the magnetic -Fm(q'^) 
form factors are normalized as 

Fn(Q) = l. Fr^(Q) = Or.. F^.(Q) = ij.T.— (176) 



where m is the nucleon mass, Qd and /^d are quadrupole and magnetic moments of the deuteron, re- 
spectively. The tensor polarization components of the final deuteron are expressed through the deuteron 
form factors as follows: 



T20 [A + 5 tan 



2 "Cl 



i=[^ryFcFQ + Ir^'F^ + iry(l + 2(1 + r^) tan^ 



T21 [A + 5tan2 |] = + ^2 gij^2 ^_^y/2p^p^ 1^ 



(177) 



T22 [A + Btan 



2 '^e 



2v^ 



Equation (|173|) can be obtained by using the standard technique |^2| from the following amplitude 
of the process 



(178) 



where Um{l) denotes the free electron spinor with 4-momentum / and spin projection m, and q = 
I — I' = P' — P is the 4-momentum transfer, P{P') is the initial (final) deuteron momentum; \DA4) is 
the deuteron state with total angular momenta projection Ai, and is the electromagnetic current 
operator. 




Figure 5: Electron-deuteron elastic scat- 
tering in the one-photon approximation. 



The deuteron current matrix element is usually parameterized in the following way (due to P- and 
T-parity conservation and gauge invariance): 



{D'M'\J,\DM) = - eCM'{P')^PM{P) 



{P' + P),{9-'F,{q^)~'^^F,{q^) 



(179) 



where ^m{P) and ^X^/(-P') are the polarization 4- vectors of the initial and final deuteron, respectively. 
Form factors Fi^2{q'^), Gi{q^) are related to Fc(g^), Fq^{q^) and Fyi[q^) by the equations 



Fc = Fi + -r/[Fi + (1 + 77)^2-^1], Fq = Fi + (l + r/)F2-Gi, Fm = Gi. (180) 



The normalization condition for the deuteron current matrix element has the form (in contrast to (BSI) 
i is included into the matrix element): 

lim {D'M'\J^\DM) = 2eP^ 5 mm' ■ 

To calculate the deuteron form factors one should use the particular system of reference. In the 
laboratory frame the 4- vectors have the following form (the Z-axis is along the photon momentum): 



P = (M, 0), P' = (M(l + 2t]), 0, 0, 2M^^/1 + r]), 

r, — (OA/T^ n n OA/T /I _1_ nn\ 



Cm=+i{P) = ^M=+i{n = --^(0,1,^,0), 
^M=-iiP) = U=-i(P') = ^(0, 1, -t, 0), 



U=o(P) = (0,0,0,1), U=o(P') = (2v^V'l+^>0,0,l + 2r/). (182) 
Using expressions ( |181| , p.82| ) and the parameterization in the form of Eq. (|179|) one obtains: 



{M'\Jo\M) = 2Me (1 +r/) {f^5mm' + 2r/[Fi + (1 + 7^)^2 - Gi]5m'o5mo}, 
{M'\J,\M) = ^ 0y /TT^ Gi {Sm'm+i - Sm'm-i]- (183) 

To calculate the deuteron form factors, one should know three matrix elements with different total 
angular momentum projections and current components. 

3.2 Gauge Invariance and Gauge Independence in the Bethe-Salpeter 
Approach 

It is well known that the principle of gauge invariance imposes stringent constraints on the amplitudes 
of electromagnetic interactions with bound systems. In the first order of perturbation theory in the 
charge e, this principle leads to the continuity equation for the electromagnetic-current-density operator 
and to the Ward-Takahashi (WT) identity for the five-point Green's function. We will consider the WT 
identity for the five-point Green's function and its implications for the Mandelstam current, which 
determines the amplitude of electron scattering on deuterons in the Bethe Salpeter formalism. 

The continuity equations for the bare (Noether) current and effective currents (for example, the 
Mandelstam current or conserved currents in nonrelativistic quantum mechanics) are not sufficient for 



ensuring the gauge invariance of the amplitudes of electromagnetic transitions ||82|-||86|. Moreover, 
both initial and final states must correspond to the current used in the analysis. 

As a rule, the conserved deuteron electromagnetic current includes two-particle contributions as- 
sociated with meson-exchange currents or interaction currents. At the same time, it was shown in 



^, that the amplitude of elastic eD scattering can be gauge invariant in the relativistic im- 
pulse approximation, which is based on the concept of the one-particle scattering mechanism. This 
result was obtained in the BS formalism with one boson exchange (OBE) potentials and with separable 
interactions. At first glance, it is at odds with the common point of view on the problem of gauge 
invariance 



Our objective here is to study in detail the conditions under which the gauge-invariant description 
of elastic electron scattering can be achieved in other models for constituent interaction. In addition, 
we discuss a certain extension of the WT identity 38 1 for an arbitrary system of charged particles. 



For a two-fermion system, the Mandelstam current can be represented as 
where J^^'' = Jfj;^^ and are, respectively, the one-particle, and two-particle contributions that 



satisfy the relations ||82| 



< Jf'^\k\ k; P', P) = eiS{k' -k- q/2) [S^^\P/2 + k)'^ - S^^\P' /2 + k'Y^] S^^\P/2 - k)'^ + 
+ e26{k' -k + q/2) [5(2)(P/2 - k)-^ - S^^\P' 12 - k')-^] S^^\P/2 + ky\ (185) 

Jl^\k', k; P', P) = eiV{P, k' - q/2, k) - V{P', k', k + q/2- )e, + 
+ e2V{P,k' + q/2,k)-V{P',k',k-q/2)e2. (186) 



Here, k and P {k' and P') are the relative and total 4-momenta in the initial (final) states, respectively; 

P' — P _L .Q(l,2)/'t.N ;„ ^-^,„ J 



V{P, k', k) is the kernel of the BS equation. This kernel depends on the relative momenta k and k' and 
on the total momentum P. In the isospin formalism, we have 



e. = (187) 

where e is an electron charge, and Tz{i) are the Pauli matrices (i = 1, 2). 

It is worth noting that equations ( p.85|) and ( |186|) do not define the current completely. They only 
impose certain constraints on the longitudinal component of the current. 

The amplitude of elastic eD scattering can be represented in the form (see (|178| )) 

M'm,m = e^l 0^,^,^''''{P',k')Uk',k,P',P)^'''{P,k) ^ e'^iMUi),, (188) 

where $-^^(P, fc) ($-'^'(P', A;')) (^ is the BS amplitude, which describes the initial (final) state, and 
e'^ is the virtual-photon polarization vector. The gauge-independence condition 

I'iMUi), = q^mi,^)^'^ + {MU,Y^^] = (189) 

is met if the current satisfies identities ( p.85| ) and (186) and if the amplitudes $'^*^(P, k) and $'^^^'(P', k') 
satisfy the BS equation with the same kernel (|52D. 

In the relativistic impulse approximation illustrated in Fig. the deuteron matrix element Ai^^^, 
defined by the Eq.( |189| ), for on-shell 7A^A^ vertex satisfies the gauge-independence condition ||91| 



9"(-^M'A/)f ^ = 0. (190) 

It implies that 

riMU,)^ = J ^^$^^'^'(P', A:')jf (fc', A;,P',P)$^*^(P, A;) = 0. (191) 

It turns out that, for certain models of interaction, this relation holds. Let us isolate isospin structure 
in V. We have 

ViP, k', k)= ^tVAP, k\ k), (192) 
r=o,i 

where 11^ is the projection operator onto the state with total isospin T, and Vt is the corresponding 
component of interaction. In the ladder approximation, 

VT{P,k\k) = VT{k' -k), (193) 

we can verify that the condition ( [L91| ) is satisfied for any function Vrik' — k) . This result is independent 
of the isospin value in the initial or final state, although only one Vo{k' — k) component contributes to 
(|191|) in the case of eD-scattering. We note that all OBE interactions can be represented in the form 
(|193|) and that the corresponding amplitudes are gauge independent, in accord with the result reported 
in 



B5, 86 



As the next example, we consider the interaction described by a separable potential ( |108|) . The 
authors of |^ Q proved that the amplitude in the impulse approximation is a gauge-independent 
quantity in this case. Their proof is based on the transformation properties of the quantities in Eq.( |191|) 
under boosts. 

However, this result can be obtained in a simpler way by calculating the contraction in ( |189| ) and 
by considering that, for interaction of the form ( |108| ), the BS amplitude in ( |191| ) is independent of the 
total momentum P (see [pl| ). 

Thus, we can conclude that, in the ladder and separable approximations for the kernel of the BS 
equation, the amplitude of elastic scattering in the impulse approximation is gauge independent. The 
common feature of the kernels in this section is that they are independent of the total momentum of 



These conclusions do not mean that other gauge-independent contribution can always be disregarded 



(for example, the contribution of the exchange currents). Indeed, the results reported in |9^, p3 
where the qq system with separable interaction was studied, revealed that two-particle currents exert 
a noticeable effect on the pion charge form factor, especially at high momentum transfers. We can 
state with confidence that, for a certain model of A^A^ interaction, two-particle and more complicated 
electromagnetic currents must be taken into account in a consistent manner even in calculating the 
elastic nuclear form factors. However, there is no universal recipe for constructing these currents for a 
given interaction of particles in a bound system. 

3.3 Relativistic Impulse Approximation 

In the relativistic impulse approximation illustrated in Fig. |^, the deuteron current matrix element can 
be written as 



D'M'\j;^'^\DM) = te J ^Tr|x^-'(P',fc')rf (g)x^^(P,fc)(P-7/2-fc-7 + m)}, (194) 



where (P) is the BS amplitude of the deuteron, P' = P + q and k' = k + q/2. The vertex of 'jNN 
interaction, 

rf (,) = 7„ff ) - ''''''';J'''''' F<\/), (195) 

is chosen to be the form factor on mass shell. The isoscalar form factors of the nucleon Pi2^(g'^) = 

(P^2^(g^) -|- P]^2^(g^))/2, which appeared due to the summation of two nucleons, are normalized as 

F^^\o) = 1/2 and F2^\o) = {xp + >Cn)l'i' with Xp = /ip — 1 and x„ = /i„ being anomalous parts of the 
proton (/Xp) and neutron (/i„) magnetic moments, respectively. 




Figure 6: Relativistic impulse approx- 
imation for the elastic electron deuteron 
scattering. 



It should be noted that the choice of the 'jNN vertex in a form (|195|) and disregard the current 
interaction (CI) (two-nucleon currents), in general, breaks the gauge independence of the reaction. 
Nevertheless, as it was shown in there exist two cases when RIA and CI are gauge independent 
separately. The first case corresponds to the separable interaction with no dependence on the total 
momentum in the radial functions, and the second one to the one-meson exchange kernel. 

First, the trace was taken in Eq. ( |194| ). The covariant form for the BS amplitude ( p6[ ) was used. 
After taking the trace, the scalar products of 4-momenta (P, k, q) and the deuteron polarization 4- 
vectors {iM^iXi') with definite spin projections were inserted. Then, using Eqs. (0), functions hi were 



expressed in terms of functions (i)jLSe (^^e section |2.6.2|) . All scalar products were evaluated in the 
laboratory frame. 

The resulting expressions for the deuteron current matrix element can be written as 
{D'M'\Jr\DM) = I^:^{q') Fi'\q')+I^'^{q') {q% (196) 



^d'^'(?') = / dko \k\' d\k\ dicose)J2K,s'^K^ |fc'l)0..s.(^o, \k\) I^'/m'm, i^o, |fc|,cos^,g^ 



L'S'e',LSg 



where the function Ii2^ij^ uiko, \k\, cos 6*, g^) is a result of the trace calculation and the substitution of 



^J■ 



In Eq. ( [L96| ), the radial part of the BS amplitude for the final deuteron (pii'S'g'iK^ W\) depends on 
the components of the 4- vector k' calculated in the rest frame {rf). The vectors in the rest frame and 
in the laboratory frame are related via Lorentz transformation: 



where the Lorentz transformation matrix C is of the form: 



(197) 



C 



2^y^TTv \ 

1 
1 

\ 1 + 27] ) 



l + 2r] 





(198) 



To simplify the notations, the components of the 4- vector k'^^ are denoted by k' = k'^j- = (/cq, k'^, k'y, k'^), 
and = = ./k^TK^TW- Using relations (|l8l]), and (|9|) one finds 



k' 

k ^ 



1 + 27])ko - 2^a/1 + r]k^ - Mr], 



1 + 2r])k^ - 2^^Jl + r]kQ + M^^Jl + 77 



(199) 



where fco, k^, ky, k^ are the components of the 4- vector k in the deuteron rest frame. 



3.4 Numerical Results 

The calculations were performed with a covariant separable kernel of A^A^-interaction Graz II (for details, 
see section p.8.3| ). Considering only ^S^ and ^D^ states one writes (see Eq. ([70|)) 



1LS+ 



{ko, \k\) = S+{ko, \k\; s)g^^g^{ko, \k\), 



where g^^j^g^ is the radial part of the vertex function. Thus, the Bethe-Salpeter amplitude involves 
singularities in the fco plane, which are infinitesimally close to the real axis. Some of the singularities 
arise from the propagator, while the others come from the radial part of the vertex function — in other 
words, from the functions g]^^^ defined in Eqs. ( |159| ) . 

For the initial deuteron, the singularities do not depend on (or if) and always remain in the same 
quadrant: 

for the propagator k^ = ±M/2 =1= V k"^ + rn^ ± ze, 

ILS+ u _ _L. /i,2 , ^2 



for the functions g^^ ko = ±y k + f3f + ie. 

The situation changes for the final deuteron. Due to the boost of the arguments ( |199D of the 
amplitude, the singularities depend on (or rj) and can go across the imaginary axis and appear in 
another quadrant (mobile singularities). The positions of the singularities are the following: 

for the propagator 

A;o = -(1 + 4r])M ± \lk^ + ni^ + 4A/r/(l + rf)M\k\ cosO + 47^2(1 + r;)^^ ± ie, 
for the functions g^^^ 



ko = -riM ± ^ k^ + + 2^/^i(lTv)M\k\ cose + ri^{l + ri)^M^ T ie. 

The mobility of the singularities does not affect the calculations if the Cauchy theorem is applied. 
But for the Wick rotation procedure, this means that the additional contributions (the residues at 
these mobile singularities) should be taken into account. The minimal value of — for which the 
imaginary axis is traversed is: for propagator — = M(2m — M) ^ 4.17 x 10~^ (GeV/c)^, for g]^^^ 



negligible (about 1%) in the region —q^ < 2 (GeV/c)^, and are getting larger with growing momentum 
transfer. But the contribution of the residue from the propagator is very large and can modify the 
curves significantly even in the region —q^ < 2 (GeV/c)^. 

The contribution of the residue from the propagator is shown in Fig. |^ for the functions A{q^) and 
B{q^). The contribution is substantially large both for the function A{q^) and the function B{q^) (for 
the function B{q^), this contribution fills the minimum, which does not exist in the experimental data). 
This result can be considered as a specific relativistic effect caused by the Lorentz transformation in 
the arguments of the Bethe-Salpeter amplitude (vertex functions and propagator). The comparison 
of the relativistic (RIA) and the nonrelativistic (NRIA) calculations with the separable kernel of the 
A^A^-interaction Graz II is shown in Fig. ^ also. One can see that the difference between the RIA and 
the NRIA calculations rather small. However we must keep in mind that the RIA calculations does not 
include the contribution of the P-waves to the deuteron (further development of the RIA calculations 
must take into account this contribution) and the NRIA calculations does not take into account the 
mesonic exchange currents. Only after that we can draw the conclusion about the full contribution of 
the relativistic effects to the structure functions of the deuteron. 

Yet another interesting result of the investigations is the dependence of the deuteron form factors 
on the nucleon form factors — in particular, on the neutron electric form factor G'^{q^). The electric 
and the magnetic form factors of nucleons (G'e('?^) and Gm('?^), respectively) are related to the Dirac 
and Pauli form factors (Fi(g^) and ^2(5^), respectively) by the equations 



G^{q') = F,{q') + ^F,{q 
Guiq') = F,{q') + F,{q'). 



(200) 



Three sets of the nucleon form factors are used in the calculations. The first set is so-called dipole 



fit 



Gl,{q') = (1 + «:p)G|(g^), = K^Gl{q'^ 



GP(g2) = l/(l-g70.71(GeV/c)^) 



2x2 



Gl{q') = 0, 



(201) 



where = 1.7928, k,^ = —1.9130 are the anomalous magnetic moments of the nucleons. The second 
set was suggested by the vector meson dominance model (VMDM) 

,2 



Ff (t) 



Ff (t) 



— 



7a 



m 



ml - g2 



l^i^lu, + + Kp 



A? 



+ Ai + g2^ 



^2L 



A? 



(202) 



2L, 



A? + g2 VAi + g2 



-q 



,\n {Xl-q')/Xl 
InAi/Ai ■ 



where Ai = 0.795 (GeV/c), A2 = 2.27 (GeV/c), A3 = 0.29 (GeV/c), = 0.163, 7^ = 0.411, = 0.784 
(GeV/c). The third set is that from relativistic harmonic oscillator model (RHOM) 



(203) 



E 



2m' 



;I^'\q'), 



GW) 

1 + K,p 

I^'\q') 



\ r\ 1 



O /l-> Q\0 



exp 



aa = 0.42(GeV/c)l 



The first model assumes that the neutron electric form factor is equal to zero. Two other models lead 
to a nonzero G^. 

Figure || shows the charge and quadrupole form factors (Fc(g^) and Fq^q"^), respectively). The zero 
of the form factor Fc{q'^) is in the range of —q^ = 1.2 — 1.35 (GeV/c)^, but experimental data yield 
— = 0.69 — 0.83 (GeV/c)^ This dip comes from the specific choice of the separable Graz II 

kernel (in the calculations with non-relativistic Graz II potential, the zero of Fc(g^) is shifted too). The 
nucleon form factors do not shift the zero in the form factor Fc{q^). The nucleon form factors with the 
nonzero electric from factor for the neutron (VMDM and RHOM) are more suitable for the description 
of the experimental data on the quadrupole form factor FQ(g^). We stress that the contriburion of the 
P-waves (the negative energy states of the BS amplitude) for the deuteron must shift the dip of the 
Fc{Q^) and must give better agreement with the experimental data analogously the NRIA calculations 
with the taking into account of the mesonic exchange currents |Jl^. The structure functions Alq"^) and 
B(q^) are shown in Fig. y. The RIA calculations of the structure functions of the deuteron show the 
strong dependence on choice the form factors of the nucleon. After taking into account the contribution 
of the P-waves for the deuteron we can make more exact conclusion about the choice of the nucleon 
form factors. Figures 10 and p!T| a show the tensor polarization components T2o{q^),T2i{q'^) and T22{q^^ 



for the final deuteron calculated with 6^ = 70°. It is seen that T2o(g^) and T22{q^) have very weak 
dependence on the nucleon form factors, but for T2i{q^) it is more pronounced. The T2i(g^) calculated 
with three different electron scattering angles 9e = 19.8°, 70°, 80.9° is shown in the Fig. [TTlb. The 
change in the electron scattering angle affects sizably the component T2i(g^). Figure |ll|c shows the 
result of calculation of the component T2i(g^) with different nucleon form factors at electron scattering 
angles fixed in experiment. This result can be used to choose between the models for the nucleon form 
factors. Unfortunately, large uncertainties in experimental data prevent from choosing between the sets, 
and future measurements of the component T2i(g^) can be very useful for the analysis. 

Note that the calculated T2o{q'^) differs from the experimental data in the region —q^ > 0.6 (GeV/c)^. 
This fact, apparently, could be explained by several reasons. It is necessary to improve the description 
of the zero of the charge form factor Fc{q^) by changing the separable kernel of A^iV- interaction and 
taking into account the negative energy states of the Bethe-Salpeter amphtude for the deuteron. It is 
also important to investigate the contribution of the two-body electromagnetic current. Information 
about the effect of these factors provide a powerful tool in the study of the on- and off-shell behavior 
of the nucleon form factors in elastic electron deuteron scattering. 

3.5 Electromagnetic Moments of the Deuteron 

In the previous section, the electromagnetic form factors of the deuteron were discussed. The relativistic 
effects in the form factors were analyzed in detail. The sensitivity of the polarization observables to the 
models for the nucleon form factors was discussed too. 

Now we consider static electromagnetic characteristics of the deuteron — magnetic and quadrupole 



moments — in details. As it was shown in ref. ||109||, the relativistic effects in BS approach, such 



as effects of relativistic kinematics, retardation effects or negative energy states, play essential role in 
evaluation of static moments. Our main task is to obtain full expressions for the moments in RIA and 
compare results with nonrelativistic ones. 

3.5.1 Definitions 

We define the magnetic moment, /iD, and the quadrupole, Qd, moment of the deuteron from the 
normalization condition ( |176| ) of the for form factors 

Tfl 1 

/iD = — limPM(g'), QD = ^limPQ(g2). (204) 

Analytical expressions for the moments can be obtained from the deuteron form factors. We use the 
Breit system which is defined as 

P^=P' = Er.= \/m^+P^. P = --. P' = -. (2Qh) 
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Figure 7: (a) Structure function A(g^). Long dashes represent the calculation without the contribution 
of mobile singularities (no mobile singularities - NMS). The solid curve shows the full relativistic im- 
pulse approximation calculation. Short dashes correspond to the nonrelativistic impulse approximation 
calculation (nonrelativistic Graz II potential). Experimental data are taken from ||96[-||9^,[^0|, (b) 
Structure function B{q^). Notations for the curves are identical to those of Fig. Experimental data 



are taken from 96, gl noq, 101, 72 . 
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Figure 8: (a) Charge form factor Fc{q^). Long and short dashes represent calculations with the VMDM 
and RHOM nucleon form factors, respectively. The solid curve corresponds to the dipole fit. The 
experimental data are taken from ref. [^. (b) Quadrupole form factor FQ(g^). Notations for the 
curves are identical to those of Fig. |^a, and the experimental data originate from the same source. 
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Figure 9: (a) Deuteron structure function A(g^). Notations for the curves are identical to those of 
Fig. Experimental data originate from the same source as in Fig. |^. (b) Deuteron structure 
function B{q^). Notations for the curves are identical to those of Fig. [pD. The experimental data 
originate from the same source as in Fig. |^. 
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Figure 10: (a) Tensor polarization component T2o{q^) calculated at 9e = 70°. Notations for the curves 



are identical to those of Fig. Experimental data are taken from ref. |7g] (values from |73, 74|,| 102 |- 
108|| were recalculated at 9e = 70°). (b) Tensor polarization components T22(g^) calculated at 9e = 70°. 



Notations for the curves are identical to those of Fig. |a. Experimental data are taken from ref. ^ . 
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Figure 11: (a) Tensor polarization component T2i(g^) calculated at 9e = 70°. Notations for the curves 
are identical to those of Fig. |^. Experimental data originate from the same source as in Fig. ^b. (b) 
Tensor polarization component T2i(g^) calculated at different angles 9e = 19.8°, 70°, 80.9° (solid curve, 
long and short dashes, respectively) with dipole nucleon form factor. The data points are from the 
same source as in Fig. |ll|b. (c) Tensor polarization component T2i(g^) calculated with dipole fit (V), 
VMDM (o) and RHOM (o) at experimental electron scattering angle values. The data points originate 
from the same source as in Fig. |l( 



Taking the vector q in the Z-direction we can write the 4-vectors as 



p = M( v^TT^, 0, 0, -07), p' = M( v/TT^, 0, 0, 

g = (0,0,0,2MV^), 

1 



iM=+i{.p) — iM=+i{.p') 



V2 



(0,1, z,0), 



(206) 
(207) 



iM=^i{P) = iM=-i{n = ^(0, 1, -I, 0), 

iM=o{P) = (-v^, 0, 0, yiT^), iM=o{P') = (v^, 0, 0, yiT^). 

Taking into account the relations between the form factors Fi, F2, Gi and Fc, -Fm, Fq 
the formulas (|206[ - p07|) and the deuteron current parametrization (|179|) , we obtain 



e {P'M' = +l\J,^\PM = Q) 

e (P'M' = 0\Jo\PM = 0) - jP'M' = +l\Jo\PM = +1) 
2M 2^i^/TTv ' 



and using 

(208) 
(209) 



With the help of the definitions ( p04| ) and the last two equations, the moments /id and Qd can be 
written as 



Qb 



m 



^ lim 

1 



{P'M' = +1\J^\PM = 0) 



lim 



{P'M' = 0\Jo\PM = 0) - {P'M' = +l\Jo\PM = +1) 



(210) 
(211) 



2M3 v-o 2ri^TTri 

These formulas are essential in the calculation of the deuteron moments. To find the analytical 
expressions we refer to the RIA deuteron current matrix element (|194|) . It should be noted that in 
contrast to the BS amplitude in the rest frame (r/), the moments are calculated in the Breit system. 
Therefore, the BS amplitudes and the propagators should be transformed (boosted) from the rest frame 
to the Breit frame. These quantities have the form (a ■ 7 = 0/^7^) 



X-\P.k) = K{C)x''^{PrfX-'k)K-\C), 



-JM'lTDl U'\ 



X—\P',k') = A~\C)x''^\P,f,Ck')A{C), 
A-\C){P ■ 7/2 - A; ■ 7 + m)A(£) = (P,/ ■ 7/2 - C-^k ■ 7 + m), 

where A is the Lorenz transformation of the BS amplitude, 

M + P ■ 770 



A(£) 



v/2M(Fd + M) 



(212) 
(213) 
(214) 



(215) 



CPrf, P' = C ^Prf, P = '^Prf, p' = C ^p'rf- The matrix of the Lorenz transformation is 



and P 

written as (compare with Eq. ( |198D ) 



C 



f v^TT^ \ 

10 

1 

V -v/^ v^rq^/ 



(216) 



Applying expressions ( pl2[ - pl4| ) to Eq. ( |194| ) one finds 

d^k 



{P'M'\J^'^\PM) 



I e 



Tr\x''^'{P^f,k')r,{q)x''^{Prf,k) 



where k' = k'^r is connected with k = krf and q as (components are given by ( |199| )) 



k' = C^k + ^Cq. (218) 

The boosted 'jNN vertex has the form: 

f^{q)=A{C)T^{q)A-\C). (219) 

To find the expressions for the moments we should expand Eq. ( pl7| ) in rj up to the order y/rj for the 
magnetic moment and the order rj for the quadrupole moment, and then take the hmit rj —>■ following 
Eqs. ( piCI| , 211 ). As it follows from Eq. (pi7|) , there are objects which give term of order t] {^), namely: 

(i) the BS amplitude of the final deuteron, including vertex function itself and propag ator x^^'{Prf,k'); 

(ii) boosted 7A^A^-vertex T^{q) ; (iii) transformation operator of the BS amplitudes [A^^(£)]^. 
Analyzing these factors we could note the following peculiarities of the current matrix element, which 

come from the covariant and the relativistic nature of the formalism: (i) effects connected with the 
dependence of the BS amplitude on ko ^ a/A;^ + m^; (ii) effects connected with Lorenz transformations 
of the BS amplitude arguments (|216|) ; (iii) effects connected with Lorenz transformations of the BS 
amplitude itself ( ^15| ); (iv) effects connected with negative energy states and transitions between ^S"^-, 
and ip(^)'(°)-, 3p(°)'(-).waves. 
We note here that A(£) is written in the limiting case, 77 ^ as 

A(/:) = 1-^7370 + | + O(r/^/'). (220) 

3.5.2 Magnetic Moment of the Deuteron 

From Eqs. ( |210[ ), ( pi7| ) and ( |220| ) we obtain the following result for the deuteron magnetic moment: 

/iz) = /i+ + yUi_ + yU2- + /is-, (221) 

where the matrix elements corresponding to the transitions between states with positive energies are 
denoted by subscript +, while the subscript — means that at least one negative energy state is included 
in the matrix element. These terms /I's have the form 

3 1 

/i+ = (/ip + /i„)(P35+ + P3D+) - ^if^p + yUn - 2'>^^D+ + (222) 

= i(/Xp + /i„)(P3pe + P3po) + i(P3pe + P3p.) + ^(Plpe + Plpo) + , (223) 

/i2- = -ifip + l^n)Pzs- + ^35- + ^(/ip + lin)PzD- " ^^3^- + ^2-, (224) 

/i3-= C^i^^^'i^'^'^', (225) 

LSq,L'S'q' 

where Puse are the pseudo-probabilities of the corresponding states defined as 
Pilsq = 1, 

LSg 

^^^^c, \S,{ko, \k\) g,,s,{ko, \k\)\', (226) 



UJ+ = 2Ej^ - M, = -2Ej^ - M, uJe = uJo = -M, Ef^ = Vk^ + m^. (227) 



In Eqs. ( |222[ - |225D , the diagonal terms are given explicitly while nondiagonal ones are included into /i3_ 
and P's, which are defined as 



1, 2m, 3o+ m 3c.+ m 3c+ 2m, 
B.^ = — u„ + u„ - 1 ^ ff, 1 UJ"^ UJ"^ - 1 P,<.+ 228 



+ -^l/ip + /in - 1 - ^j^^l 

1 2m 1 1 2m 

^1- = -^il-j^)if^p + f^n + -)iP^Pf + P^P;)--il-j^)iRpS+Rps) (229) 

2 1 pe 1 DO 4 1 pe 1 DO ; — AiTTl^ 3 pe 1 po ^ — 3 po 1 pe 

--H/^ ' + -i^io ^ ■ + V^(/ip + /in - 1 + -^)H/' ' ' + +V^(/^P + f^n - 1)H/^ ' 

\/ 2, 3 po 1 pe / — 3 po 1 pe 

1 2m 3 c- rn 3 c- TTi 3 c- 1 4m 3n^ 

i?2- = --(Mp + /in - 1 - — )i//^ - — ifo ^ + —H.^' - + u„ - 1 + —)H,''' (230) 

^ rr'D- , m ■^D- , 3 2m 72 "^No-'sr,3Dr 

iin H Hn H 1 -r 3n- H /in + /in — 1 H -ni • 

Functions C^lsqal's'q' ^nd i/,^^ ^ ^ '^^^^ are given in ref. [p,lCI|| . Rewriting equation for /i+ we find 

/i+ = ^NR + A/i+, (231) 

where the equation 

3 1 

^^NR = (/ip + /in) - -(/ip + /in - 2)^='l?i+ ' 

gives the nonrelativistic formula for the deuteron magnetic moment, while the expression 

Afi+ = R+- ifip + iin)iP3s- + P-^D- + P'Pl + P'P? + P'Pl + P^P?) (232) 
represents the relativistic correction. The final expression for the deuteron magnetic moment is 

/iD = ^^NR + A/i , (233) 
A;U = R+ + Afi^ + /i3_ , 

A/i_ = -(/ip + /in) 



^(P3p. + P3p.) + (Pipe + Pp.) + 2P35- + ^Pa^- 



+ i(P3pe + P3po) + i(Pp. + Pp.) + P3^- - ^P3^- + R,_ + R,_ . (234) 

We introduce the nonrelativistic value for the magnetic moment of the deuteron /iTvj? and group all the 
relativistic corrections into Afi term. 

3.5.3 Isoscalar Pair Current 

We want to understand the physical meaning of the P-waves in the BS amplitude. To this end, we 
consider nonrelativistic reduction of the expressions of the magnetic moment. Starting from Eq. (233) 



we perform fco-integration in functions h'^^^s^l s q ^ QiLSg,iL's'g' _ Ugj^g the Cauchy theorem, we 

take into account only the pole of the positive energy propagator component S^{ko, \k\; s) (see Eq. ([7TD), 
namely ko = ko = M/2 — Ej^ + ie. We disregard the contributions from the singularities of the other 
parts of propagator and vertex functions. 

We now introduce the nonrelativistic analogs of the Bethe-Salpeter vertex functions for ^Si- and 
^D^-waves, corresponding to the positive nucleon energy, 

g3a+{ko,\k\) g3r,+ {ko,\k\) 



with ai = 47rv2M, and the normahzation condition for the nonrelativistic wave functions 

J ed\k\{u\\k\)+w'{\k\)) = ^. (236) 

Performing the integration over ko and introducing nonrelativistic analogs of the functions gsg+ 3^+ , 
we represent the terms appearing in expression 

S+,D+ , S+.P , D+,P /OQVN 
/^D = /iD + /^D + /^D • (237) 

The first term (pd = 2/7r J k'^d\k\w'^ is the probability of D-wave) 



(/ip + /i„) - ^ (/ip + /in - ^) Pd, (238) 



reproduces nonrelativistic value for the magnetic moment in impulse approximation, while terms 

/if'"' = 8M(27r)3 + Di + D2 + D, + D, + D, + D,), (240) 

correspond to matrix elements between ^S^—, ^D^ - and P-states. The functions Si and Di are given 
as 

Si = ^ j \k\'d\k\(^g^peCko,\k\)-g.poCko,\k\)y{\k\), (241) 
5'2 = (/ip + Ain - 1) y \kfd\k\ (^g3p.{ko,\k\) - g3po{ko,\k\)^u{\k\), 
S3 = -^ I k'd\k\(^g^p.{ko,\k\)~grpoCko,\k\)')u'{\k\), 

Do = -ifip + fin-l)^ j \k\''d\k\{^g.p.{koAk\)-g,po{kQ,\k\))w{\k\), (242) 
^1 = ^^ / 1^1 ^1^1 1^1) ~ 5'3p°(^o, |fc|)) w(|fc|), 

^2 = / l^l^^l^l (^3p^^(A;o, |fc|) -^3po(^o, |fc|)j u;(|A;|), 

D3 = ^j k'd\k\ (^gsp.Cko,\k\)-g^po{ko,\k\))w'{\k\), 

^4 = ^1^ y (^5(ipe(A;o, |fc|) - 5'ip°(^o, |fc|)) w(|fc|), 

D5 = ^ J \k\^d\k\ (^gipe{ko,\k\) - gipo{ko,\k\)^w{\k\), 

^^ = ^1 \k\'d\k\ (^g.peCko,\k\)-g.po{ko,\k\))w'{\k\), 

where u'(|fc|) = du{\k\)/d\k\, w'{\k\) = dw{\k\)/d\k\. In the formulas shown above, we perform the 
expansion in k/m to terms of order (|fc|/m)^. It should be noted that the matrix elements in (|233|) 
which correspond to transitions between ^5*^- and ^Z^^^-states, ^S^-, ^D^- and P-states vanish as the 
result of integration with respect to /cq upon performing the expansion in |fe|/m. 

The term /i^ ' coincides with the contribution of nonrelativistic impulse approximation. To con- 



One-Iteration Approximation. The Bethe-Salpeter equation for the deuteron is commonly solved 
iteratively. After a partial-wave decomposition, the BS equation reduces to the set of integral equations 
for the vertex functions gjisgiko, \k\). They are related to the partial amplitudes 4>jLSg{ko, \k\) which 
are used here for simplicity on the right-hand side of the following equation 

2 +00 +00 ^ 

9.ako, |fc|)=E j<j^^^\d\k'\Y^ ^l....(^o, Ife'l) W(^o, Ifc1),(243) 



L'S'q' 



where the index jj, labels the type of exchanged meson, q^nn is the coupling constant and Vj^g^ jj^^ig,^, 
is the matrix element for the interaction kernel between the states JLSg and JL'S'g'. To get fast 
convergence to the solution of this equation one needs a good educated guess for the initial vertex 
function gjisei^o, \k\) and respectively 4>jL'S'e'{Ky \k'\)- The solution of the equivalent nonrelativistic 
problem may be helpful for this task. After several iterations one usually gets correct solution. 

For a zero-order approximation, it is convenient to relate the vertex functions for the ^S^- and 
^Dj'"-states to the corresponding S- and D-wave components and respectively) of the 

nonrelativistic deuteron wave function as 

4+(^0' 1^1) = - ^(1^1)' (244) 

9i%{ko, \k\) = -a, (M - 2^^) w{\k\), (245) 

and to set the vertex functions for all the remaining states to zero: 

9TLS,ik,, \k\) = 0, JLSg Si^'Df. (246) 

Note that we need the vertex functions at a fixed value of the relative energy ko = Jcq = M/2 — 
(see Eqs. ( p41| , p42| )). We substitute relations ( p44H246|) into the right-hand side of Eq. ( p43| ) and perform 
integration with respect to fco- As before, we do this with the aid of Cauchy theorem, taking into account 
only pole of the propagator, which corresponds to positive nucleon energy. This yields 

JLSg=^ P^.}P^, (247) 

where we introduced the functions Vi'tl ,,^^,^^A\k\,\k'\) that can be obtained from the functions 

'^jLSg (35+ 3_D+)(^0' 1^1' ^0' l^'l'' setting /cq = ko and = Ivq and by expanding Ej^ and Ej^' in powers 
of \k\/m to second-order terms. 

We consider now the vr-meson exchange kernel in the ladder approximation and calculate in this 
case the functions gsp^^ipo in first iteration. The matrix elements of the interaction kernel in this case 
are given by (index superscript (/x) is omitted), 

V3pe3s+iko, \k\; kQ, \k'\) = V^peij^+iko, \k\; k^, \k'\) = 0, (248) 
V^poss+ih, \k\] fed, \k') = ^ Q 1^1 (% - ^) Q^i^) + ^ \k'\ Qi{z) - i |fc| {Ef^ + 2m)Qo{z)^ , 

^3p.3B+(^o, \k\; k'„ \k'\) = ^ (^3 \k\ {2Ej^, + m) Q^iz) - m \k'\ Q^{z) - - \k\ {Ey - m)QQ{z 

/2 / 2 1 

l^ipe3s+(/co, k'^, \k'\) = - — |fc| {Ey - m)Q2{z) + m \k'\ Qi{z) - - \k\ {2Ey + m)QQ{z 

V^pe3j,+ {ko, \k\; k'o, \k'\) = ^ Q 1^1 (% + 2m)Q2(^) - m\k\ Qi{z) - ^ \k\ {Ey - m)Qo{z)^ 



where Q£{z) are the Legendre functions of the second kind and z = {k^ + fc'^ + A*^ - {k'^ — ko)'^)/2\k\\k'\ 
with being the vr-meson mass. In the above expressions we have omitted the isospin factor which is 
—3 for the deuteron. The functions Vusg,iLS+{\k\, \ k'\) then become 



V^3pe35^-(|fc|, Ifc'l) = V,p.,^^{\k\, k'\) = 0, (249) 

— 2t?7 

V^3p.35+(|fe|, \k'\) = -j={\k'\ Q,Cz) - \k\ go(5)), 



(250) 



\/2m 

y3po3z,+ (|fc|, \k'\) = Q2{S) - \k\ Q,Cz)), (251) 

Vipe3s+{\k\, \k'\) = -^V^posilkl, \k\), 



(252) 



Kpe3B+(|fc|, \k'\) = V2V,po,j,+ {\k\, \k'\), (253) 

Kp„,35+(|fc|, \k\) = Kp.,3^+(|/c|, \k'\) = 0, (254) 



with S = {k^ + k'^ + ^l)/2\k\\k'\. 



To perform the integration over \k\ in Eq. (|247|) , we introduce functions u{r) and w{r) in the 



configuration space, which are analogs of functions M(|fc|) and w(|fc|) {wq = u, W2 = w, £ = 0,2), 

+ 00 +00 

we{\k\)= j rdrw,{r)n{\k\r), ^ = ^ j k"" d\k\w,{\k\) n{\k\r), (255) 







where j(,{x) are the spherical Bessel functions. Using formulas ( |247|) , ( P49| )- (p5^) and Eq. (P55|) , and 
recalling that the Legendre functions admit the integral representation, 



+ 00 

Q,{i) = 2\k\\k'\ j e-^'-''rdrj,{\k'\r)ji{\k\r), 



we obtain 

g,p.{h,\k\) = Q, (256) 



+ 00 

,2 



^3p.(fco,|fc|)= -(-3)-|^^ j rfr^(l + /x.r)(^r) + i=^(r))ji(|fc|r), (257) 



+ 00 

g.p.{hAk\)= -(-3)^%^ / dr'-^{l + ^,^r){-^u{r)+w{r))jmr). (258) 



g,po{k^,\k\) = 0, (259) 

where the isospin factor (—3) has been taken into account explicitly. 

Summarizing the above analysis, we can say that, by using the one-iteration approximation and by 
specifying the zeroth-order approximation for the vertex functions with the aid of the relations ( p44| - 
P46| ), we related the P-wave components of the Bethe-Salpeter vertex function to the nonrelativistic 
deuteron wave function. Note that the resulting relations are proportional to g'^j^^/An. 

Pair Current We should note how to relate the P-wave components with the magnetic moments of 
the deuteron. Let us now make use of the expressions obtained above. Substituting relations ( p56| - |259|) 
into Eqs. ( |241[ - p42| ), we perform an expansion in terms of g'^^^^/An. As a result we find: 



A^D =NaA [l^p + M TT^ Da - tt;-^ Da-—Db + -r^ Dc + Dd , 



where Na-^ = ai/8M(27r)^. We have also introduced the integrals 



and functions 



Sa 
Sb 
Da 
Db 
Dc 
Dd 
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dr 



g3po{r) [u{r) - ru{r)], 



- / drgipe{r) [u{r)], 



71 f dv 

— g3po{r) [2w{r) + rw'{r)], 
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TT 

2 

7r 
2 

7T r dr 



drgzpo{r) [w{r% 
drgip_e{r) [w(r)], 

gipe(r) [2w{r) + rw'{r)] 



V3 - 4^ ^ (1 + /^.r) (^0 + -^w{r)), 
m Att r a/2 



^?ip.(r) = v^^%^^(l + ^.r)(-i=«(r)+«;(r)). 
^ m 47r r a/2 



(261) 



(262) 



(263) 
(264) 



It is convenient to rearrange the above expressions in such a way as to isolate explicitly the component 
corresponding to the vr-meson pair current. For this purpose we single out three terms: 



8M(27r)3 

(/ip + /in) 



(-(/ip + /in) + (/ip + /in) ^ Da 



(265) 
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47r 4m 



[I + /i^r) 



—=uw w 

V2 10 



The first term is proportional to [fip + /i„) and coincides with the vr-meson pair current of the deuteron 
magnetic moment. The second and third terms have no analogs in the nonrelativistic expressions. 
The final expression for the deuteron magnetic moment has the form, 



f^D = /i^'^ + /^d"^' + ^/^ 



PC 



D: 



A (2) I 



,(3) 



(268) 



where /i^^^ is the result obtained in the nonrelativistic impulse approximation (see Eq. (|3|)). /ig*^ is 
the contribution of the vr-meson pair current, and A/i£) is the additional relativistic contribution to the 



3.5.4 Quadrupole Moment of the Deuteron 



We derive now the expression for the quadrupole moment of the deuteron. Using Eqs. ( ^111 ), ( |217| ) and 
(|220D we get 

a, a' p,p' a, a' p,p' 

Here we label partial-waves of the deuteron BS amplitude by a, a' with p-spin p, p'. The subscripts C 
and M mean the corresponding contributions of charge and magnetic parts of the 'jNN vertex ( [L95| ). 
The terms with the index LB indicate the contributions due to the factor ([A(£)~^]^ — 1) in the matrix 
element (|217| ), which come from the Lorentz boost. 

The matrix elements of the operator Q in Eq. ( P69| ) schematically represent the combinations, defined 
by Eq.( ^ll| ), of the following matrix elements (P = P^/ = -Pr/' ^ = ^rf, k' = k'^^, Jvi = 0, 1): 



c 



2M J i{2n) 



Tr {xm{P, k')-foXM{P, A;)(P ■ 7/2 - ■ 7 + m)} 



{PM\J^^^\PM) 



C 



V^7 



d^k 



2M J i{27r] 



4 Tr {xm{P, k'hoXAiiP, k){P ■ 'y/2 - k ■ 'j + 7/1)7073} 



{PM\J^^^\PM) 



M 



(269) 
(270) 

(271) 



e X 



d^k 



2MAm J z(27r)4 
{PM\J^''^\PM)^J' = 



Tr {xMiP, k')A{C){^oq ■ 7 - 9 " 77o)A(/:)Xa^(P, A;)(P ■ 7/2 - A: ■ 7 + m)} 



(272) 



e >fy^ f d^k 



■ Tr {xMiP, k')i'joq -1 -q- 11q)xm{P, A;)(P • 7/2 - /c ■ 7 + 771)7073} 



2M Am J i{2TiY 
with X = Kp + x„. 

We manipulate the above expressions to find the quadrupole moment, where the main contribution 
to the quadrupole moment comes from the transition between the positive energy components ^S^ and 
^D^ of the BS amplitude. 



a,a'=S,D 



(+.+) , n(+'+) 



(273) 



The explicit expressions for g^'^^ and Ql^^'^> are given in Appendix A. Performing integration in 
Eq. ([A.l|) with respect to /cq and taking into account only the positive energy nucleon pole in the 
propagators we obtain the following expressions. 



Qd 



d^k 
20 J (27r)3 



A; 



,dui\k\) dwi\k\) 
d\k\ d\k\ 



+ 3\k\w{\k\) 



dui\k\) 
d\k\ 



In the last equation we have introduced the expansion over \k\/m up to terms of the second order and 
used substitutions by following Eq. (|235|) . We then find that Eq. (|274|) coincides with the nonrelativistic 
expression. 

Second term Q^}^'~^^ in Eq. ( P73| ) and matrix element from Lorentz transformation (see (|270| )) have 
purely relativistic nature and are responsible for relativistic contributions to the quadrupole moment of 
the deuteron. For instance, for positive energy states after integration of the Eq. ( |A.3| ) in parts we get 

.2 



e idkoI^d\k\ M 2k^ 2 1 \ 2El-mE^ 
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3.5.5 Numerical Estimates 



In this section we presented analytical calculations of the magnetic and the quadrupole moments of 
the deuteron. We showed that both moments in one-iteration approximation have three terms. The 
first is pure the NRIA contribution, the second is related to pair-mesonic current and the third is 
pure relativistic which has no analogy with the NRIA calculations. To perform numerical the RIA 
calculations of the both moments of the deuteron first of all we must take into account the contribution 
of the P-waves in the deuteron. It will be done in future. Now we estimate the contribution of the RIA 
calculations ^5, 116 | into the magnetic moment of the deuteron. 



There are eight states in the deuteron channel (instead of two in the non-relativistic case), viz. Si, 

(see section |2.6.2|) . The normalization condition for this functions 



3 n+ 3 c- 3 p)- 3 pe 3 po 
J-^l ) ^1 ) ) 1 ) 1 ; 

can be written as: 



3 pe 3 po 
' 



P+ + P- 



P_ 



P35++ + 



(275) 



introducing pseudo-probabilities Pa that are negative for the states ^S^ , ^D^ , ^Pf , ^P°, ^Pq, ^Pq , and 
positive for ^S"^, '^D^ []85| , p6| . The calculation with realistic vertex functions give the following values: 



25+1 TP 


+ 




3Pf +3 P° 


3 pe _i_3 po 




[P%] 


4.8 


-6 ■ 10-4 


-0.88-10-2 


-2.5 ■ 10-2 


m 



It is obvious that the main contribution to the normalization is due to the states with positive energies, 
and the contribution of the P-states is larger than that of the negative energies states by at least one 
order of magnitude. 

Now we calculate the relativistic corrections /x^*" and AfiD to fin ( [2681 ). Thus | |116| , |117| |, 



= 2.792847337; 



0.0022; 



-1.9130427; 



exp 



NRIA 



0.857406; 
-0.00058; 



,,NRIA 
H'D 



0.852458703; 



= 0.8538396861. (276) 



We have shown that the expression for the magnetic moment in the Bethe-Salpeter approach can be 
written in a form closer to non-relativistic calculations. The additional terms in equation ( p68| ) can be 
considered as relativistic corrections to the non-relativistic formula. 

The non-relativistic value reflects only the Z)-state probability. Whereas in the relativistic corrections 
P-states play the dominant role and improve the agreement with the experimental data. 

The magnitude of the corrections can be compared with the contributions of mesonic exchange 
currents to the magnetic moment as extracted from ref. [|1^]. The main contribution is due to the pair 
term, which leads to A/i//i7v_R = 0.21 — 0.22% for different forms of the Bonn potential. The same size 
of this correction as compared to ( p76|) may be considered as an indication that both corrections are of 
the same physical origin. 



4 Deuteron Electrodisintegration 

We examine here in particular the nonrelativistic reduction of the Bethe-Salpeter approach. We do not 
attempt to solve this problem in a general way, but merely consider a particular process, namely the 
electro-disintegration of the deuteron, in order to investigate some relativistic corrections. We consider 
the threshold region only, where a single transition amplitude to the ^Sq final state dominates (see, 
e.g., ref. [|113|| ), which offers some technical simplifications. The extension to other partial waves is 
straightforward but tedious and beyond the scope of the present discussion. 

In fact, the deuteron disintegration reaction near threshold energy region ||114|| invariably attracts 



attention from both the theoretical and the experimental sides. It has been an excellent process to 
examine non-nucleonic degrees of freedom and relativistic effects. It is well known that lA alone, which 
considers nucleons as non-relativistic objects, fails to describe the double differential cross section at 



present in the results of calculations [|115|| . The minimum can be filled if meson-exchange currents are 
calculated. The contributions of tt-, p-currents, and AA-configurations further provide a satisfactory 
agreement with data ||118| , |119|| . The question of a consistent inclusion of all relativistic corrections (at 
least for the pion sector) has recently been addressed in ref. ||120|| and supports the above statement. 
Nevertheless, some conceptional problems of the theory still remain open questions. Among them is the 
problem of gauge invariance and the choice of the nucleon form factor to be used in the exchange- current 
contributions 



121 - 124 



By now it became apparent that, some relativistic methods reproduce such corrections as meson- 
exchange currents in the impulse approximation. It has been shown that in the framework of the 
light-cone approach, the so-called extra components /s of the deuteron wave function, and g2 of the 
^S'o-wave function introduced in ref. |79] give an expression in the nonrelativistic limit, that equals to 
the contribution of the pair current [^. In this context, the first iteration assumption that will be 
explained below for solving the dynamical equation is substantial. Another important result is the 
calculation of the static electromagnetic properties within the Bethe-Salpeter approach. We find, that 
the contribution of the P-states to the magnetic moment of the deuteron is numerically close to the 
contribution of mesonic currents and agrees in sign [16, 110]. It has already been noticed earlier, by 



considering covariant reductions of the BS equation to the three dimensional equation, that negative 
energy components in the wave function are responsible for pair-current-type contributions, see, e.g., 
ref. [11251 



4.1 Relativistic Kinematics 

The amplitude of the electrodisintegration of the deuteron, Mfi, in the one-photon-exchange approxi- 
mation has the form (Fig. [T^), 

Mf, = te^ u{l\ <) ru{l, Se) ^ {np\J^\DM), (277) 

where u{l, Sg) is a spinor of the electron with 4-momentum / and electron spin 4-vector Se, q = I — I' 
is 4-momentum transfer, the matrix element of hadronic current {np\Jfj,\DAi) represents a transition 
from the deuteron state \DAi) with 4-momentum P and total spin projection Ai to the final state of 
a rap-pair with 4-momentum P' = P + q, and is the electromagnetic current operator. 

^ /' 

Figure 12: Deuteron electro- 

disintegration reaction, eD e'{np), in 
the one-photon approximation. 




Considering the final np-pair being near the threshold, we can assume ^S'o-state for the np-pair as 
a basic state, and because of the Lorentz covariance, transformation properties under parity and time 
reversal as well as current conservation, the general structure of the transition matrix element for the 
1"*" — > O"*" transition current can be written as 

{npCSo)\MDM)=ze,^p,^JUqf'P^V{s,q^), (278) 

where is the deuteron polarization 4-vector, V{s, q^) is a Lorentz invariant function. Due to the 
asymmetric tensor, e^Q,/3^, the matrix element is gauge invariant, 

A T in A//\ — n 



Using Eq. (|277|) we obtain the differential cross-section in terms of tlie leptonic and liadronic tensors, 



Li^i, and W^'^, respectively. 



«2 

g4 |Z| 



(279) 



with a = e'^/An. The leptonic tensor, L^,^, is 

Lf,u = 2{lf,ll + I'^Q + q^g^, + 2imee^,^pq'' s^, (280) 
where Se is the electron spin 4-vector. The hadronic tensor, W^y, has the form 



(27r)3 2Ep^i2nY 



(281) 



iy'^^ = (np(iSo)|JlZ^A^)(DA<|jt-Mi5o))|^ / 6iP + q-pp-p^)^ 

where Pp (p„) is the 3-momentum of the proton (neutron), and Ep = \/p\, + vn? {Ep^ = a/p1"T"?7?) 
is the energy of the proton (neutron). With the help of Eq. ( |278| ) the hadronic tensor can be written as 



W'' = R G"'' p^p G*'''^ Vis,q 
where i? is a kinematic factor (s = (P + q)"^), 

87r2 2M ' 
and Pai3 is the density matrix of the deuteron 



\k* 



Pa/3 - ^\-9aP + -^^) + l^^^a^l^ 



D 



Pr.P 



(282) 



(283) 



(284) 



,2 PxiPx2 



D ' 



where (W^a)^/? = i^a/S'yxP'^ / M , sd is the spin 4-vector and is the alignment tensor of the deuteron. 
Using the explicit form of the deuteron density matrix ( P84| ) , the hadronic tensor becomes (the electron 
mass is neglected) 



W^i"^ = [qA^^M^ - (P^f) + (^M?- + %P>^){Pq) - Pf^Puq' - q,quM'] \ v{s, q^] 

^RM (sDq) te^^ap q" P" \v{s,q^)^ , 



R 



■^[^tiXiaP^XiU^S + ^fj.XiafS^XiU'ysjP^P'^q'^q^ + ^(-S'AiAa + ^2 



(285) 
(286) 
(287) 



X [g,,{q'M-' - [Pqf ) + {P^q, + q^P,){Pq) - P^P.q' - q^q.M'] p^^"^ V{s, q 



The superscripts (m, v, t) denote unpolarized, vector-polarized and tensor-polarized parts, respectively. 
For the case of unpolarized electrons and deuterons the differential cross section can be written as 



dE'dVt' 



e / unpol 



/da\ M\k*\m ^.^ ^ ^, oOe^ o^e 

(^)m^2^^ 1'^' + - 2^-^'™ i' i 



where 



cos^ 



^ da\ 

Kdin/M "4^2 sin^% 



(288) 



(289) 



is the Mott cross-section. We use the following normalization conditions, 

< DP'M'\DPM >= 2Ed {2-Kf 5mm'5{P' - P) 
< Np\ n'lNp, a >= 2Er. (27r)=^ - p). 



4.2 Asymmetries for the Polarized Deuteron 



With the general form of the hadronic tensor W'^'^ in Eqs. ( |285D -( PS7D we can calculate various polarized 
deuteron asymmetries, 

^ da{],D)-da{lD) 

da{\,D) + da{lDy ^''^> 

where da is the differential cross section, f (|) is the helicity Ae = 1) of the incident electron 

and D the polarization state of the vector-polarized deuteron. We take the momentum of the incident 
electron directed along the Z-axis. The 4-vectors / and /' then take the form, 

/ = (Ee,0,0,Ee), /'= (E;-E;sin^e,O,E^COS0e), (291) 

where 6e is the electron scattering angle, and I' is in the XZ-plane defined by the incident and the 
outgoing electrons. 

We consider, first the case of the vector polarized deuterons. If the deuteron is polarized parallel to 
the Z-axis, then vector polarization asymmetry is 

, ^ 3 {E, + E',){E,-E',co^e:) 

" 2 (Ee + E^)2-2EeE;COs2^e/2 ' ^ ' 

where k is the degree of the polarization of the deuteron. For the asymmetry, the dependence on 
y(s, q^) disappears. For the case of the backward scattering {9^ = 180°) Eq. ( p92|) gives 



3 

All = 2^' ^'^^^'^ 
If the deuteron is polarized parallel to X-axis, the vector asymmetry is 

A -3 , {E^ + E'^)E'^sine^ 

^ 2 {E, + E^)2 - 2E,E'^cos^ej2- ^ ' 

Generalizing formulas for the asymmetries to the arbitrary direction of the polarization {d,ip), we 
write 

,3 {E, + E'^) (E; sin 9^ sin ^ cos + {E^ - E'^ cos Oe) cos 
'^^ " 2 {E, + E':)^-2E,E',cos'^ej2 ' ^^^^^ 

We consider then the case of the tensor polarization of the initial deuteron. If the initial deuteron 
is only aligned due to the Pdzz component, then 

da{pDzz) = da[l + A^^pozz], (296) 
4^2 ^E'J + AE^E'^ - AE^E'^ cos 6^ + 3El ^ cos 29^ 



A. 



mE, + E'^y~2E,E'^cos^ej2) 



where A^z is the tensor analyzing power. For the backward scattering {6^ = 180° the analyzing power 
is Azz = 1. 

4.3 Relativistic Impulse Approximation 

We evaluate the electromagnetic current matrix element {np{^ So)\Jfj,\DAi) in the relativistic impulse 



approximation (see diagram in Fig. jTS]). Within the framework of the BS approach, using the Mandel- 
stam procedure of constructing of the electromagnetic current operator p3|, it can be written as 
(compare with Eq. ( |194| )) 



(np( 



'So)\J!!'^\DM)=z ! rf^fcTrly°0(P^fc')rW(g)y,rP,fc)(P-7/2-fc -7 + ^)1, (297) 



where x^'^i.P'^ is the BS amphtude of the ^5*0 state of the np-system, Xjm^-^^ ^) is the BS amphtude 
of the deuteron, and k' = k + q/2. The vertex of the 'yNN interaction is chosen to be on mass shell, 



The isovector form factors of the nucleon Ff^\q'^) = {F^2{(f') ~ -^i 2^(Q'^))/2, appear due to summation 
of the two nucleons. They are normalized as F^\{i) = 1/2 and F2^-'(0) = (xp — x„)/2. 



Figure 13: Deuteron electro- 

disintegration in the relativistic impulse 
approximation. 




P'=P + q P 



In order to extract V{s, q^) from the general expansion ( p78| ) for the 1"^ — > 0'*' transition, we first take 
the trace. According to the Lorentz structure of Eq. ( p^97| ), seven integrals have to be evaluated, which 
are expressed through Lorentz scalar quantities abbreviated by (...) (the index M. for the polarization 
vector ^ is suppressed for simplicity). 



X7 



ie,e.p, Cq^P^ j d'k (...), 
^e^.ap, CP" j d^k (...) k\ 
j d^k (...) k\ 
eV^^ j d'k (...) k'^ks, 
iea^-ysCq^P'' j d^k (...) k'k„ 

teap-ysCq^'P'Pf. j d^k (...) k\ 

teap-yseq^P^q>.J d'k (...) k'. 



(298) 



From inspection of Eq. ( p98|) it becomes clear that beside Ti, which needs no further consideration, two 
generic types of integrals have to be evaluated, namely. 



d^k (...) k^ = Ciq^ + C2Pa , (299) 
j (fk (...) k^kp = DiM^go^p + D2qaqp + D^{q^Pp + P„g^) + D^P^Pp 



/3 ' 



where we have already indicated that after integration the expressions should depend on the external 
4-momenta only, i.e. the transferred momentum q and the total deuteron momentum P. Due to the 
antisymmetric tensor, the number of terms is reduced substantially, and all the terms with Di, i > 1 
vanish. So, beside Xi, it is necessary to evaluate the following integrals: 



= ^e„,.,^PV f d'k(...)c,(k,q,P), 



Xs = te^af3-yCq^P"' J d^k{...)c2{k,q,P), 
X4 = te^^p,Cq''P^ j d''k{...)d{k,q,P), 

= le^p^^Cq^P^ M"" J d^k (...) dik,q,P), (300) 



where the functions q = Cj(fc, q, P) and d = d{k, q, P) are given by 

M^{kq) - {Pk){Pq) 



Ci{k,q,P) 



M'^q^ - {Pqf 



_ (Pfc)^g^ + M^(fcg)^ + {Pqfe - M'q'k' - 2{Pq){Pk){kq) 
''^^'^'^^ " 2M\M^^-{Pqf) • ^^^^^ 



The integration and the comparison with the structure of the transition matrix element given in Eq.( p78|) 
finally leads to the expression for V{s,q^), 

V{s,q^) = V,{s,q')F^^\q') + V,{s,q')F^''\q'). (303) 
The expressions for Vi^2{.s,q^) are 

Vi{s,q') = Id'k{vl'\s,q^k) + V^'\s,q^k) + Vj'\s,q^k) + V}'\s,q',k)), (304) 

V2{s,q') = I d'k{vl'\s,q\k) + Vi'\s,q',k) + V<'\s,q',k) + vi'\s,q',k)). 

(1 2) 

The explicit expressions for the functions ' are given in Appendix B. 
4.4 Isovector Pair Current 

We calculate now the function Vi^2{s,q'^) (Eq. ( |304| )) from the underlying dynamics in the Bethe- 
Salpeter approach. We have two representations of the BS amplitude for the concrete calculations: 1) 
the covariant representation and 2) the partial-wave representation. Implementation of the covariant 
representation for the initial and the final states allow us to get rid of the functions Vi^2{s,q'^) from 
the proper current matrix elements. However, in order to study the relation to the nonrelativistic 
expressions, it is more suitable to work out the formulas for Vi^2{s,q'^) in terms of a partial-wave 
representation of the vertex functions. Hence, the expressions for Vi^2{s,q'^) can be written as 

Vi^2{s,q'^)=i d'^k gjLSsiKAf^'l) (^JLSeJ'L'S'g' 9J'L'S's'iko,\k\), (305) 

JLSg,J'L'S'g' 

where the partial vertex functions gjsLg{kQ, \k'\) are connected to the partial amplitudes (fco, \k\) 
via simple relations ([70|) . These vertex functions represent the states ^Sq ...^Pq of the np-pair, and the 
states ^S^...^P° of the deuteron (see paragraph p.6.2| ). In general, the lengthy functions Ojslq^j'l's'q' 
depend on the Lorentz scalars and the explicit expressions are omitted here. We will specify them 
below after having introduced appropriate approximations. The expressions for /cq and are given in 
a formally covariant way (compare with Eq. (^)) 



1/2 

(306) 



We note that through the use of the Bethe-Salpeter vertex functions the denominators of O appear- 

( M 10 4-1-^4-1^., -\-n^\ ( A/T 10 -\-l-^-r 'FP, -r I . fZ lO -\- -\- 'fP. . -\- 



and (a/s/2 ± /cq =f Ej^^ =f ie), that stem from the nucleon propagators (see Eq. ([Til)). We evaluate the 
integrals in the laboratory frame (deuteron at rest). At the threshold, because of the small deuteron 
binding energy, it is possible to utilize the static approximation that preserves the analytic structure of 
Eq. ( |305| ) through the following equations, 

k'o = ko, \k'\ = \k + ^\ = + ^ + \k\\q\x^ . (307) 

Here and in the following we use a = a/|a|, 

\q\ = {{M^ + s-q'^y -AM^sy/y{2M), x = k-q, 
q = (tu,q), w = (s-M^ -g2)/(2M). (308) 

We illustrate the approximation by looking closer at the expressions involving the nucleon positive 
energy states. The full denominator of the integrand in Eq. ( |305| ) leads to a complicated pole structure 
and reads explicitly 



It can be simplified by using Eq. (|307|) : 



(^y^ + ko-Ef^+ze^ (^y^-ko-Ej^+ze^ (^ + k,-E^+it 



ie . (309) 



Thus, the static approximation means in particular that uo = Q (no retardation) and the Lorentz boost 
transformation of the np-pair vertex function is neglected. We can go beyond the static approximation 
by expanding the full expression in terms of uo/M and \q\'^ / s that leads to additive corrections. 

The integration over fco can now be performed using the Cauchy theorem, namely by choosing a 
proper integration contour and specifying the corresponding poles, e.g., closing the upper half plane 
leads to poles for k^ at k'^ = \/s/2 — E^i. The vertex functions are then evaluated at {k'^.k). Since 
in the reaction under the consideration s ~ Am? ^ M^, we expand the vertex functions near ko = 
M/2 — Ef^ for gi and Gi, respectively, that allows us to derive analytical expressions in the one- 
iteration approximation. The analogous procedure holds for other partial wave vertex functions. With 
this choice, one of the nucleons in the deuteron is taken on shell. 

We now perform the ko integration for the interaction in Eq. ( P05| ). The angular integration is 



simplified by taking q along the Z-axis, and by replacing x in Eq. (^081) as follows: 



k"- 






|g|V4 




k\ 







2\k\' 

X = ^—^^1^, dx = d\k\', (310) 

Finally, to go to the nonrelativistic limit we introduce a |fc|/m-expansion excluding terms of the order 
Oik^/m?) (i.e. Ej^ = m + 0{k'^/m^) = Ej^'). The resulting structure functions Vi,2(s,g^) are then 
given by 

+00 |fc|+|qr|/2 - ^ - 

^ 1 f f 9^s+{KAk\) /^3s+(A;o, |fc|) 1 c/3i3+(A;o, l^r 







\\k\-\q\/2\ 

2 y/^-2Eu' + ie 



k 

/ ^3g+(fco, |fc|) |g| +2|fc|x 1 g-iD+{ko,\k\) \q\p^{^x) + 2\k\x \\ 



g3pe{ko, \k\) -g3po{ko, \k\)^ - — (^-^3pe(fc[„ |fc'|) +^3po(fc^, |fe'|)^ 



\k\+\q\/2 

^ 1 f \i \ f \r \u\, \'i 9is+{koAk\) /g3s+{ko,\k\) 1 ^3^+(A;o, |fc|) 



||fc|-|qr|/2| 
yS |q|2 j35^(^Ufc1) / ,A X 

3v^ 1-12 



\q\ 
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''\M-2Ef^ \q\\k'\ ^ M-2Ef^ \q\\k'\ J j ' ^^^^^ 

where = po = y/s/2 — Ej^^, and ko = M/2 — Ej^, and ^2(2;) = (Sx^ — 1)/2 is the Legendre polynomial. 
The functions (735-, S's^^-, fl'ipf, 5'ip°, ^'15- disappear in the above expressions after /cq integration and 
because of the \k\/m expansion. Within this approximation, we are left with (+) to (+) and (+) to 
(e, o) transitions only. All other matrix elements, such as (— ) to all and (e, o) to (e, o) vanish. We now 
examine the expressions for Vi^2(s, q"^) more closely. To recover the nonrelativistic result, we neglect the 
vertex functions gape, gspo, gape, gapo that correspond to the negative p-spin components (i.e. do not 
exist in the nonrelativistic scheme). If we replace the functions g3g+ and g-Api+ by the nonrelativistic S 
and D wave functions of the deuteron using Eqs. ( P35|) and gig+ by the ^Sq continuum wave function 
of the np-pair in the following way 

gig+{k'f^,\k\') ^ I 1 

-a2Uo{\k'\), a2 = ^= — , (313) 



and insert them for the respective vertex functions in Eqs. (|311 , 312 ) we obtain 

+00 |fe| + |q|/2 

^^°n^,?') = ^^G'i^^('?') / \k\d\k\ J md\k'\uo{\k'\){um-^H\k\)P2{x)) , 

||fch|qr|/2| 

(314) 

where we have introduced the mag netic isovector form factor G^^^ = F^^^ + F^. In configuration 
space, the respective integral is found using the following transformations for the scattering state (for 
the deuteron states see Eq. ( P55| )) 



hoo +00 



Uoi\k\) = ^ I rdruo{r)U\k\r), ^= J k' d\k\ Uo{\k\) U\k\r) . 


The resulting expression is 







+00 

V^'^(^'l') = ^^GZ\q') I druo{r)(u{r)U\q\r/2)-^w{r)U\q\r/2)). (315) 







This result reflects the so-called nonrelativistic impulse approximation and represents the lowest order 
nonrelativistic expansion of the transition form factors given in Eqs. ( pil| , pi2| ). 



One-Iteration Approximation. Since the one-iteration procedure for the deuteron channel was 
discussed in detail in section |3.5.3| , we consider here the ^S'o-channel of the np pair. 

The inhomogeneous Bethe-Salpeter equation for the amplitudes in the ^S'o-channel reads 

2 +00 +00 ^ 

9f,NN -i f f 1 1^1 



^ —00 

,jL'5'o'(^o, \k\; /cq, \k'\)(f)jL's'e'{K^ Wl), (316) 



where /x is the exchange meson and 4>^j\sg denotes the plane-wave function, 

,(0) 1 1 



0y/5.(^o, \k\) = ^ ^jLs.^^st ^(ko) m\ - (317) 



47r ^0 k 

and k* is the on-energy-shell momentum given by = ^J s/A — m?. 

For the subsequent discussion it is more convenient to spht the system of equations as follows: 

2 +00 +00 ^ 



ifci) = *Sj(*«.ifci)+ Elf 5 / < I ii^w'"'' 

^ —OO 

2 +00 +00 



X E ^izkjL'5V'^^°' 1^1' \m,,,Ak'o. W\l (319) 

As in the case of the deuteron channel, we consider the one-iteration approximation. To this end, 
we chose a similar expression for the zero approximation which is helpful in finding relation with the 
nonrelativistic solution (compare with ( |244[ - p46| ) ) , 

0L5,(fco, |fc|) = 0, LSg St . (321) 

Here Uo(|fc|) is the nonrelativistic continuum wave function in ^5*0 channel given by 

1 mt{\k\Ak*\\Et.A 
um) = ^ m - \k*\) + feJ2_ fe2 ^^ ^ (322) 

and t{\k\, |fc*|; Ej^*) is the nonrelativistic half-off-shell t-matrix for the ^5*0 channel normalized through 
the condition 

t{Ei^.) = t{\k*\,\k*\;Ef^*) = -- sin5oe*'°, (323) 

where 5q is the phase shift, and Ej^* = k*'^/m. 

Analogously to the deuteron case, we finally arrive at the one-iteration solution for the amplitudes 
in the ^5*0 channel, 

gsp.ikoAk'\) = -(+1)^4^ [ rfr^(l + /i.r)«o(r)ji(|fc'|r), (324) 
" irm 47r y r 



gspoCko,\k'\) = 0, (325) 



where (+1) is the isospin factor. 

We have shown in this section that proper choice of zero approximation wave function (i.e. the 
nonrelativistic one) allows one to obtain additional partial amplitudes through the Bethe-Salpeter 
equation after one iteration only. They are connected to the interaction kernel and in electromagnetic 



Pair Current. We are now in the position to turn to the first order corrections to V^^\s,q'^) given 
in Eq. ( |315|) . To this end, we expand the expressions given in Eqs. ( |311| , pi2| ) into a power series of 



sInn/^'^ to extract the pionic contribution only. Also we consider the P-states contribution only, i.e. 
components with one negative p-spin. The resulting expression for the transition form factor V{s, q^) 
will be denoted by V^^\s, q^). Substituting Eqs. Ol^) as well as Eqs. (|324| , |325|) into Eqs. ( ^ ) and 
and using the replacements of Eq. (|235|) and Eq. (PTl) we obtain 



+00 _ +00 |fc|+|qf|/2 

^"'<^''''> = 2^^^%" / + / 1*1 ''1*1 / I^VI*'! (326) 

||fc|-|qr|/2| 

x|(-3) Uo{\k'\)(u{r) + ^w{r)) jmr) 

UQ{r) ii{\k'\r) 



The ^'-integration can be solved analytically which yields 

+00 

y''\^y) = 5^%^ / *^(l + A.r) (327) 







Uo{r) {u{r) + -^w{r)) ji{\q\r/2) 



F^'il') - ~, FPif)) Mr) («(r) + ^»(r)) Ji(|9|r/2) 



hoo 



H{q'f-f^l rfr^(l + /i.r)«o(r)(t.(r) + l=w{r))j,{\q\r/2). 




m'^M\q 

Here we have introduced the function 

H{,') = Fri^')-l^^,FPiq^). (328) 

This first order contribution in s^^tvat supplements the lowest order relativistic expansion as given in 
Eq. ( |315| ). We then arrive at the following expression for the transition form factor: 

+00 

Vis,q') = ^l^|G^)(g2) I druoir) (u{r) U\q\r / 2) - -^w (r) U\q\r / 2) 



+00 

+^^^'^ m^| ^47 /^''^^^ + f^nr)uo{r) (n(r) + ^w{r))j,{\q\r/2)^ . 


The lowest order in the ttNN coupling constant fi'^jvAf leads to an additional contribution after iter- 
ating the P-wave channel once. Comparing this result to the one achieved within the nonrelativistic 
scheme that introduces meson-exchange currents, we find that the first term coincides analytically with 
the nonrelativistic impulse approximation contribution, and the second one with the vr-pair-current 



Nonrelativistic impulse approximation and pair current contribution 

section for the D ^ ^Sq - transition has the form (nonrelativistic case) 



The differential cross 



16 1^ l'^ • 2 ^ //; I r/\2 o/r' 2 ^M/lc II mMag n 7~,\|2 

-a — — - — sm - {{I + I ) — 211 cos -)|( Sq \\ ^ \\ D)\ , 



where q 



dVLduj 3 q2 ^2 — 2 
(cu, q) is the momentum transfer, t 



(329) 



-q^. The momentum A;* is related to the relative energy 



Ej^* of the np system as given before Ej^* = k*'^ /m, and the relation between kinematical quantities is 
given by |g| = ^((2m + Ej^^y - + tf/AM^ + t and /' = {-u + ^u'^ + t/ sin^ ^/2)/2, u = E^- E'^. 

In the general case, the current matrix element is a sum of the nonrelativistic impulse approximation 
contribution and the contributions from the meson-exchange current and the retardation-currents. Our 
concern is only with vr- meson pair-current part, and hence we find for (^5*0 || T^^^ \\ D) the following 
expression: 

(1^0 II Tf II D) = {^S,\\T^t^\\D) + {'S,\\T^:!\\D) , (330) 



where we introduced T^^^, which reflects the impulse approximation operator and T™^^, which is the 
TT-meson pair (contact) operator. 



Results. We have shown that the nonrelativistic reduction of the Bethe-Salpeter approach utilizing 
the one-iteration approximation leads to the results that exhibit the same analytical structure as the 
nonrelativistic result plus pair-current corrections. Some details differ as will be discussed below. One 
may now use the "exact" nonrelativistic wave functions to calculate the different contributions to the 
cross sections. This is done here for an illustration. 

The formula for the cross section is given in Eq. ( |329| ). The dominant Ml transition matrix element 
(i.e. D npi^So)) of the multipole decomposition given there directly corresponds to the nonrelativistic 
reduction of Eq. ( |278| ) given above. The contribution of the nonrelativistic impulse approximation given, 
e.g. in ref. ||118|| coincides with the formula given in Eq. (|315|) if ai and 02 are chosen as in Eq. (|235|) 
and ( |313| ). The analytical structure of the nonrelativistic pair-current contribution equals to that of the 
P-state contribution derived from the Bethe-Salpeter approach given in Eq. (|327|) . We note that the 
nonrelativistic pionic pair-current contribution given in ref. | ]118| | depends on the nucleon form factor 
Pi^\(f)- Sometimes, the electric nucleon Sachs form factor G^{q^) = F^^\q'^) + '§^F2^\(f') is used 
(see ref. ||119|| ). The Bethe-Salpeter approach yields a different dependence on the nucleon form factor 
that is also consistent with current conservation and is given by the function H{q^) defined in Eq. ( p28| ). 
As an illustration of the different behavior we display the form factors F^\q^) (solid line), G^^\q'^) 
(short-dashed line), and the function H{q^) (long-dashed line) in Fig. [14. 
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Figure 14: Nucleon electromagnetic form fac- 
tors which enter the calculations for the nonrel- 
ativistic pair current contribution as discussed 
in the text. Sohd line displays F^^\q'^), dotted 
line — G^\q'^), and dashed line — H{q^). 



As a parametrization for the nucleon form factors we use the one given in ref. ||127|| . The form factor 



pair current is expected to be larger than the one using the other two form factors. The function H{q^) 
being in between the two others, the respective contribution of the P-states in Bethe-Salpeter approach 
is different from the nonrelativistic calculations which normally use either F[^\q^) or G''^\q^). 

To investigate the influence of the nucleon form factors more closely we calculate the impulse ap- 
proximation and pionic pair-current contributions to the differential cross section. The calculation is 



performed with the Paris A^A^ potential |128| at Enp = 1.5 MeV and 9 = 155 



It is well known that some uncertainty is related to the strong nucleon form factors. Without 
dwelling too much on that point we would like to compare three different form factors to see how this 
uncertainty propagates. The introduction of the strong nucleon form factors changes the expressions for 
the pair-current contribution |7^. Three sets of strong nucleon form factors (for vrA^A^- vertex) have been 
employed in calculations. The results are displayed in Figs. 15 a-c, respectively, a monopole vertex |^| 
with a cut-off mass of A = 1.25 GeV (set a) and A = 0.85 GeV (set b). In addition, a vertex inspired 
by a QCD analysis has been used with two parameters chosen to be Ai = 0.99 GeV and A2 = 2.58 GeV 
PI (set c). 




Figure 15: Impulse approximation (ia, solid 
lines) and vr-meson pair-current contributions 
(pc) to the differential cross section: (a)-(c) 
correspond to different sets of strong nucleon 
form factors as explained in the text. Dashed 
line: ia+pc with F^^\q'^), dash-dotted line: 
ia-|-pc with H{q^), long-dashed line: ia-|-pc 



with G|^^(g^). (d) Calculations using F^'''{q'^) 
only but different strong nucleon form factor 
sets as given in the text. Dashed line: ia-|-pc 
with set a, dash-dotted line: ia+pc with set b, 
and long-dashed line: ia-|-pc with set c. 
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It is seen from Fig. that there is a strong dependence of the differential cross section on the 
nucleon electromagnetic form factors. It is also seen that the minimum at t = 12 fm~^ in the impulse 
approximation contribution is shifted by pionic pair current to the region t > 16 fm~^ using G'e(5'^), 
to t > 18 fm~^ using the function H{q^) and to t > 22 fm~^ using F^\q^) as a form factor. The 
largest shift in the cross section is obtained where F^\q^) is used in the calculations. Both the size of 
the shift and the behavior of the cross section considerably depend on the set of parameters as well as 
the type of the irNN-vertex used. This is illustrated in Fig. p!5| d for calculations using Fi^\q^) alone, 
but for different parameterizations of the strong vertex. Thus we conclude that to compare the RIA 
calculations with the experimental data one must take into account the contribution of the P waves 
for the deuteron and include the contribution of the final state interaction within the BS approach. It 
should be done in future, in the meantime we shoewd here the important points of this calculations. 



5 Deep Inelastic Scattering 

Studies of the electromagnetic processes with finite momentum transfer to a bound state within the 
BS formalism allowed to extract important features of the relativistic bound states. It has provided 
the general relations between the relativistic structure of the bound state and its dynamical properties. 
However, such a study triggers many questions. At present, numerical calculations are restricted to the 
region of small relative momenta of nucleons which excludes from the analysis of the BS formalism a 
very important range of high relative momenta in Minkowski space. The role of the non-nucleon degrees 
of freedom should be studied and the BSA should be reformulated in order to take into account these 
states. In this section, we discuss a process which is free from such uncertainties. This is deep inelastic 



the forward amplitude of the Compton scattering, where momentum transfer to the bound state is zero. 
In analysis of this process we can concentrate on a specific feature of the bound state — relative time 
of the constituents. 

The problem of the relative time of the constituents was widely discussed for a long time, and it 
was accepted that this feature is an unphysical property of the relativistic theory and the problem 
can be solved by fixing the relative time in a consistent way. This solution produced large number of 
similar quasipotential approaches. However, the analysis of different methods of fixing the relative time 
performed in recent papers |^ has shown that these different methods are not equivalent. The analysis 
of the DIS off the deuteron and light nuclei has shown that this is a property of a bound state, and 
it allows to solve the long-standing problem of the EMC effect. In this section we will consider it in 
detail. 



5.1 Basic Definitions 

In the process of the deep inelastic scattering (DIS) of leptons from nuclei 



(331) 



a lepton I with momentum k is scattered off a nucleus A with initial four-momentum P transferring 
momentum, q = I — V . Inclusive experiments on DIS record only the final lepton with momentum 
while the state X is unobserved final hadronic states of the reaction. In the lowest order of the electro- 
magnetic constant a = e^/An this process can be depicted schematically by the one-photon exchange 
graph, as shown in Fig. |T6|. 




Figure 16: Deep inelastic scattering off 
a nucleus A in the one-photon-exchange 
approximation. 



X 



In this approximation, the cross section for the reaction (p31|) can be written as a contraction of 
hadronic and leptonic tensors: 



a 



da^—L^^{k, k')W,,{P,q) 



q- 



(332) 



The leptonic tensor describes the hard photon emission by a lepton. Since the lepton is assumed to be 
a point particle, the expression for L^'^ takes the simple form 



(333) 



All the information about the target and its structure is contained in the hadronic tensor, which 
has the form: 

W,AP. Q) = 1 J.|P)(27r) V(P + q - p„). (334) 



This definition allows the hadronic tensor of the nucleus to be related to the amplitude for forward 
Compton scattering by means of the unitarity condition: 



W,,{P,q) = ^lmT,,{P,q). (335) 

In the case of electron (muon) scattering on an unpolarized target, the tensor W^,y can be written most 
generally as 

W'^-'iP, q) = W,{u, q')g^^ + ^^^^^P^P'^ + ^'^^f\ 'q'' + ^^^^^(^V + q'Pl. (336) 

where u = qQis the photon energy, Wi are the target structure functions. Due to the gauge invariance 
condition, 

q,W^''{P,q) = Q, (337) 
the hadronic tensor depends only on two structure functions, 

^ ( .q^,qu\.W2{lW)( \fp 
W^,{P, q) = Wi{v, q ) ^-g^, + + V " ^^'^ J V " 

In the Bjorken limit {—q^ = — > cxd, u — *• oo) the condition for the scale invariance is realized, and it 
is possible to change over to the structure functions independent of g^, 

MWi{v,q^) Fi{x) 

iyW2iiy,q^)^F2ix), (338) 

where Fi and F2 are the scale-invariant structure functions (SF), and x = —q'^/{P ■ g) is a new scale- 
invariant variable, called x- Bjorken variable. Using ( |338| ), Wfj^^ can be written as 

W,u{P, q) = (^-g,u + ^) F,{x) + J- (^P, - ^q^^ (^P, - ^q,^ F^ix). (339) 

The experimental study of the deep-inelastic scattering of muons on the deuteron and the iron 
nucleus has led to the discovery of the EMC effect, which is the manifestation of the nontrivial nucleon 
structure changes in a bound state. Strictly speaking, the EMC effect, which is interpreted by most 
authors as a decrease of the value of the SF of the free nucleon in the iron nucleus in the range 
0.3 < a; < 0.7, most likely reflects the difference in the structure of the deuteron and the helium 
nucleus. In fact, if we restrict ourselves to the range 10~^ < x < 0.7, it is easy to verify that the form 
of the ratio r{x) = F^'^ / F^ is mimicked in heavier nuclei. The universality of the x dependence of the 
modification of the nucleon structure in nuclei with mass A < A was established in Refs. |[130| , |131|| , 



where the world data on DIS of electrons and muons on nuclei were analyzed. These results obviously 
indicate that the saturation of the modification of the structure function F2{x) already occurs in the 
helium nucleus. The evolution of the modifications from A = 4 to A ~ 200 is manifested as an increase 
of the oscillation amplitude a^^^ = 1 — r;^j„ by a factor of ~ 3 and is well described as the effect of 
evolution of the nuclear density function of A. 

The cessation of the modifications of -^2(2^) for A > 4 is demonstrated most clearly by the unchanging 
form of r^(a;), fixed by the location of the three points xi = 0.0615, X2 = 0.287, and 0:3 = 0.84 at which 
r'^(x) = 1 independently of A ||132|| . Thus, the detailed study of the structure function of the light 
nuclei helps to resolve not only the problem of the EMC effect but to understand the nature of the 
short range interactions as well. 



5.2 Basic Approximations 

Lets us consider the standard assumptions used in the analysis of DIS: 



• treatment of the DIS amplitude for a nucleus as an incoherent sum of amplitudes for individual 
constituents; 

• representation of the hadronic tensor of the bound nucleon in terms of scalar functions in the 
same form as the free nucleon. 

The first assumption allows us to solve the corresponding equation for a bound state wave function. 
In case of the Bethe-Salpeter formalism within the meson-nucleon field theory this assumption leads to 
the kernel of the BS equation in the form of one-meson exchange. As it was discussed in the previous 
sections, it results in a certain success in describing the low-energy properties of the deuteron and its 
elastic form factors. However, for a successful description of high energy behavior of bound state we 
need nonperturbative methods for deriving the kernel of the BS equation. As one of the methods, the 
separable form of the kernel can be used, which is beyond this approximation. As we noted above, the 
DIS in the Bjorken limit does not depend on high energy behavior of the BSA, so we may rely on this 
approximation. Below we will use Graz II potential for numerical calculations which also successfully 
describes low-energy properties of the deuteron. 

The second assumption allows to treat the squared amplitude W^^ of DIS on the nucleus as the 
sum of the squared amplitudes for scattering on individual constituents, while the interference terms 
are neglected. As it will be shown below, the justification for this is the suppression of the interference 
terms in as inverse powers of Q^. Therefore, these terms are important at small where they 
can affect the dependence of VT^, while in the Bjorken limit they can be neglected. 

The available experimental data for DIS on nuclei is mainly in the region x > 10~^ and > 1 GeV^, 
and shows that the ratio F2 / F2 is independent of . In the calculations we shall restrict ourselves to 
the Bjorken limit, where the first and second approximations are well justified. 

The third assumption allows the hadronic tensor of a virtual nucleon to be represented in the form 



( 339 ). But this representation is valid when the nontrivial differences between scattering on free and 



bound nucleon are small. There are three such differences which result in the so-called off-shell effects: 

• the impossibility of using the condition of gauge invariance for the bound nucleon in the form ( p37| ); 

• the contribution of antinucleon degrees of freedom; 

• the relative time separating the bound nucleons. 



Since it is impossible to use the condition ( |337|) , which has been formulated systematically only 



for physical particles, the expression for for a bound nucleon turns out to be more complicated 



than ( p39| ). In general, as analysis in the quark-parton model has shown ||133|| , the amplitude for DIS 
on a bound nucleon can be constructed in terms of 14 structure functions, of which only three are 
important in the high limit. From this point of view, the choice of the actual number of the SFs 
parameterizing the Lorentz structure of the hadronic tensor depends strongly on the model assumptions. 
The model independent solution of the problem has been found within the Bethe-Salpeter formalism, 
which gives in the Bjorken limit a rigorous relation between nuclear structure functions and on-shell 
nuclear constituent structure functions -Fi(x), -^2(0;) and their derivatives 0, |S3 . 



The role of the contribution of the antinucleon degrees of freedom in the structure of relativistic 
nucleus is not yet completely clear. Recent studies in the framework of the Bethe-Salpeter formal- 
ism discussed in the previous sections have shown that in the electron elastic scattering and electro- 
disintegration of the deuteron the effects of the antinucleon degrees of freedom can be connected with 
mesonic-pair currents. They are also important for describing deuteron static properties (see section 
3751) . However, as it will be shown below, their contribution to in the Bjorken limit is negligible. 

The third off-shell effect was dropped out from all semi-relativistic approaches. The different ways 
of doing it consistently leads to different quasi-potential approaches. The disregard of the effect was 
justified by the assumption that the relative time of the bound nucleons is an unphysical feature of 
the relativistic bound state. Thus all observables should not depend on this property and different 
quasi-potential approaches should produce physically equivalent results. However, simple analysis of 
the analytic properties of the off-shell nucleon hadronic tensor has shown that the relative time can affect 
observables in the DIS. Another indirect note on such possible effects was made by the analysis 



covariance. Further analysis of the amphtude of the DIS has shown that the relative time effect is 
responsible for the deviations of nuclear to deuteron SF's ratio and can provide universal understanding 
of the EMC effect for all nuclei. In the succeeding sections we will consider in details these results. 

5.3 DIS on the Deuteron 



5.3.1 BSA for Compton Scattering on the Deuteron 

Due to the unitarity relation ( p35| ), the calculation of the hadronic part of the amphtude reduces to the 
calculation of the amplitude for forward Compton scattering on the deuteron, which by the definition 
is the expectation value of the T-product of nucleon electromagnetic currents in deuteron states: 

T^u{P,q)=^ j c?W^^(D|T(J^(x) J,(0))|D). (340) 

Using (0), this definition can be rewritten in terms of the solutions of the BS equation for the deuteron 
r^(P, k) and two- nucleon Green's functions Gq^j^: 

(^r°(P, A:05(2)(P, A:i)G6,.(g, P, ^i, k,)S^,)iP, k,)T^iP, k,). (341) 

According to (p!8|), the function Gq^^ is related to the exact two-nucleon Green's function with an 
insertion describing the Compton scattering of virtual photons on a system of two interacting nucleons: 

Ge,Aq,P,k,k') = J ^^0-^Gl\P,kM)G,,M.PMM)Gl\PM,k'), (342) 

where 

G6^.(g, P, k', k)=i j d^xdSdWd^Yd^Y'e-"'y+'^'y'e-"'''e-'''^^-^'^ (343) 

X < o|T(^(F + |)^(r - 1) j,(x) j.(o)^(r' + ^)^(r - |^))|o > . 

If a specific form for G^ is assumed, the function GQ^,y can be obtained explicitly. To find the amplitude 
for Compton scattering on the deuteron in general, it is sufficient to determine the relation between 
Gg/ii/ and the expansion of Gq^i^ in terms of the functions G4. 

Expressing the functions 6*4 using (|3^, we obtain the following expansion of 6*4 in terms of G4: 



G^{P; k, k') = 5(2) (P, k) \^{27iy 6{k - k') + (344) 
+ / ^■..^G4(P;A:,A:O5(2)(P,A;O...G4(P;A;„,A:0%(P,A:0y 



Furthermore, expanding Ge^iy and substituting this expression into ( p42| ), we obtain a series whose n-th 
term has the form 

G,^;J{q,P,k,k')= (345) 
E / ft dV ^"'^^' ^' k,)G,%'\q, P, k„ k,)Gt'\P, k,, k'). 
Choosing the term of zeroth order in 6*4, we immediately obtain the corresponding contribution to 

G,^'l{q,P,k,k')= (346) 
= f^,n^-^(P. k) In'P^'' In. P k)(27r)'^r)^(k - k'\ + C^"" in. P kM2TT\^l)^ik - k' - n\] Sr..-^ (P. k'\ . 



Thus, in the zeroth order in G4 the function Gq^^ contains both the one-nucleon contribution a) 



gPIm^ P, k) = G4^. f + ^) ® ^ - ^) + ^4^- (q,^-k^®s(^^ + k), (347) 
and the contribution corresponding to scattering on various nucleons b): 

Gf^l^iq, P, k) = G,^ (^g, ^ + k^® G3. [q, ^ - ^) + ^3;. ^ - ^) ® G,, (^^ ^ + ^) ' ^^^^^ 

The Green's functions G/^f^^ and Gs^ respectively describe the Compton and elastic scattering of a 
virtual photon on a virtual nucleon. 

The first order contribution to G^^y depends on G^^^l and Gq^I, and this leads to the expression 
G,%,P,k,k') = S^2)-\P,k)G^^\^{q,P,k,k')S^2)-\P,k') - (349) 
{Si2f\P,k)G^\P,k,k")Ge^'!iq,P,k",k') +G,^;;l{q,P,k,k")+G'i\P,k",k')S^2)-\P,k'] 

where the function Gg^^^^, is expressed in terms of the Green's functions and the zero-order term 
of the function Gq^i^: 

G,^^^{q,P,k,k') = 
d^k" d^k'" 



(27r)4 (2pi)4 
d^k" 



5(2)(P, k)G,^ {q, k, k" + q) G,{P, k" + g, k'" + q)G,, {q, k'" + g, A:') 5(2)(P, k') + (350) 



(27r)^ 



Gi^^P, fc, fc")G6;2(g, P, fc", A;')^(2)(P, k') + 5(2)(P, fc)G6i°j(g, P, fc, fc")Gl'^(^, A;", A;')} 



According to Eq. ([TsD , the function Gs^ is determined by the Green's function Gs^ describing the 
absorption of a virtual photon by a system of two virtual nucleons. 

Following this procedure, we can obtain Gq^^, in any order in G4. However, the general structure of 
Eq. ( |341[ ) is such that all the higher contributions reduce to the leading term already studied. This can 
be easily checked by using (|38|) to go to higher order of G4 in (|341|) . 

Substituting the expressions obtained for G^^^ into ( ^41| ) and taking into account the definition 
(|342| ) , we obtain the amplitude for Compton scattering on the deuteron in general form: 

J4 



Tl{P.<l) = I (^rD(P,A:)G6!f)(g,P,fc)r°(P,A:) + 



d^k rl'^y rl^h" rl^h'" 

T^{P,k)S^'\P,k)G,,{q,PXk')G,{P,k',k")x 



(27r)4 (27r)4 (27r)4 (27r)4 

y.GUq.PXX')s^'\py)T''{P.k 



In Fig. we show schematically the various contributions to the amplitude of forward Compton 
scattering on the deuteron. Transpositions of virtual nucleon lines are implied for all diagrams. The 
explicit form of the expressions represented by these graphs is given by the terms (a) and (b) in ( P46| ) . 
Here the heavy and light lines denote the nucleon propagators with high and low momenta, respectively. 
The graph a) represents the relativistic impulse approximation in which only scattering off single nucleon 
is taken into account. It corresponds to contribution of the first term in Eq. (|346|) . The diagram b) 



represents contribution of interference terms in the impulse approximation (the second term in Eq. ( p46| ) . 
The contribution of the terms contain the BS vertex functions with high momenta to imaginary part of 
the Compton amplitude is suppressed as (1/Q^)', / > 2. The diagrams c) and d) represent contribution 
of interaction corrections to Gq^^, (see Eq. (|349|) . These terms contain the contributions of two or more 
nucleon propagators with high momenta and, therefore, are suppressed as (1/Q^)^ 

Thus the only independent term comes from the relativistic impulse approximation, while irre- 
ducible interaction corrections to the imaginary part of are suppressed by powers of [|135|| . 




Figure 17: Diagrams of forward Compton scattering on the deuteron 



5.3.2 Hadronic Tensor of the Deuteron 

Let us consider now approximation for the hadronic tensor of the deuteron which neglects terms of the 



order Substituting (p46|) into (|341|) and discarding the l/Q"^ terms, we can write the amphtude 

for unpolarised scattering on the deuteron as: 



^r^(p, k)s - (5 + (q,^ + k^s(^^ + k)] T^iP, k). 

The nucleon propagator S can be expanded in terms of Dirac spinors: 

m m''(p)-u*'(p) m f*(p)-u*'(p) 



(351) 



E{po-E + i5) E{po + E- 16) 



(352) 



Here E = y p^ + + {p + fn)G2{p) is the nucleon energy, which becomes the nucleon energy on the 
mass shell {E = a/p^ + m^) for 



2 

m . 



The Green's function G4^,^ is directly related to the amplitude of Compton scattering on the nucleon: 

k,q] =2m^u' I ? + k] G^,,,, ( Q,^ + k]u' I ? + k 



(- 



2m 



2 
P 



2+^1 G^^, {q,- + k]v' 



P 



P 

2" 



k 



(353) 



where N(N) denotes a bound nucleon (antinucleon). Using the representation ( p52| ) and Eq. (353) and 
taking into account the azimuthal symmetry of the Bethe-Salpeter vertex function, we rewrite the 
Compton scattering amplitude in terms of the nucleon and antinucleon amplitudes: 



^E^r'(^ + ^'^)/^'^''(^'^) = Tl{^ + k,q)f{P,k), 



(354) 



s,S 



This leads to the analog of the convolution formula for the Compton amplitude 



(0tn(| + ,,,)Ap,,) + 



{2tt 



The total averaged nucleon Compton scattering amplitude 6*4^,^ (g, y + ^) has singularity associated 

with the continuum in the intermediate state. Therefore, at large Q^, tJ^u^'^ + k,q) can be related 
to the hadronic tensor of the nucleon by the unitarity condition: 

Q)= f S^WI (^ + Kq\ f(P, k)+ f (^ + k,q) f(P, k), (355) 



where distribution functions have the form: 



f{P,k) 



im 



1 



1 



+ 



im 



-ko -{E-^) + i6 -ko + {E + ^)-i6 



+ 



-ko + iE + ^)-i6 -ko -iE-^)+i6 



The functions $ are related to the BS vertex functions as follows: 



^l4Mn,k) = C(Mz,, A;) ^<(k)llKk) 5^ n^(-k)li;(-k)rD(M^,fc), 

s s 

<^l_{MD,k) = -r°^(M^,fc)^<(k)wKk)^t;|(k)t;;(k)r°(Mz,,A;), 

s s 

^l4Mn,k) = -r°^(Mz,,A:)5^<(-k)l;K-k)^..^(-k)^;(-k)rD(M^,A;), 

s s 

^l4MD,k) = r°^(Mz,,A;)^<(-k)uK-k)5^^|(k)^;(k)r°(M,,,A:). 



(356) 



Thus, we have obtained an expression relating the hadronic tensor of the relativistic deuteron to the 
hadronic tensors of the off-shell nucleon and antinucleon bound in the nucleus. 



Structure Function of the Deuteron In order to calculate the deuteron structure function 

F^(a;), it is necessary to express the corresponding hadronic tensor in terms of the scalar structure 



function. This procedure can be performed by using the representation ( p39| ), which is valid only for 
free particles. This makes it inapplicable for a nucleon bound in the deuteron. 

The problem can be overcome by integrating ( P55| ) with respect to ko, taking into account the 
analytic properties of the integrand. The integrand contains singularities of the propagator, and the BS 
vertex functions. The nucleon propagators contain nucleon and antinucleon poles and cuts connected 
with the self energy 6*2 (p)- The latter contributes only at large nucleon energy and can be neglected 
in our basic approximation. The singularities of the Bethe-Salpeter vertex functions can be fixed by 
means of the relation between these vertices and the two-nucleon Green's function: 



r° (P, k)Tl{P, k') = \im(P' - Ml)S7^AP, k)G,,ps-y{P, k, k')S7MP, k'). 



■ (57 



(2)' 



'(2)' 



Thus the BS vertex functions for the deuteron has the same singularities in the relative momentum as the 
two-nucleon Green's function. Since the singularities in this function closest in energy are determined 
by the cut beginning at k"^ = m^, they can be neglected in integrating ( |355| ) with respect to ko, following 
our assumptions. This makes it possible to approximate the integral with respect to ko by the residues 
at the nucleon and antinucleon poles of the corresponding propagators. As a result, we obtain the 
following expression for the hadronic tensor: 



d^k 



m 



{2n)^2E^{MD - 2Ey 



<^l4Mn,k)W'i{k,q) + 



d 



{Mp - 2 E) 
Ml 



2 r 



+ {Mo - 2E)— q)^l^{MD, k)) ^^^^^ + 

+ ^1_{Md, mau{k, q) + ^U{Md, k)Wl{k, q) + 



(357) 



dko 

Ml d 



KT^^iwl{k,qW_{Mn,k)]_ 



Ml ^2 



^^i(Mn,k)WZ(k,q) 



I. 



Thus we have expression relating the hadronic tensor of the deuteron to the hadronic tensor of the 
on-shell nucleons and their derivatives near the mass shell. Now we can use Eq. ( p39| ) and obtain 
by means of a projection operator: 



Wf{q,h)=PrW;^^{k,-q,q',k^). 



In the Bjorken limit we can use the metric tensor g^^i, as this operator: 

hm 9'''W^}^\P,q) = --F^^^\x). 

This operator is independent of the relative momentum, and so the derivative of the hadronic tensor 
has the following form: 

W^ih ■ q, g^ k^) = g^^W^ih, q). (358) 

Here the first term reflects the modification of the structure of the bound nucleon. The second term 
reflects the variations of the nucleon hadronic tensor with nucleon energy and its value is proportional 
to (Md — 2E)/Md. This allows us to neglect the dependence of on kf: 



4- lim g''''wl{P,q)Lf,N 



1 ^ / N ^ d ^ , . 
-F,{x) - -—F,{x) 
x^ X ax 



Neglecting terms of order (M^) — 2EY we can write the deuteron structure function in the following 
form: 



Xd = / 7 r^T T7 TTT < Ff(xN) + — $ Md, A; V 



(360) 



Md-2E rfF2^(xN)^2,,^ M , ,E-ks,^^ d 



XN y "^ $1Md, k) + F^^(xn)— rr^(MD - 2^)— <I>^(Md, A;) 



MrmD dxti Md dko 



ko=E-Mo/2 



The normalization condition for the Bethe-Salpeter vertex function (|4^) gives two normalization 
conditions for the distribution function in this expression. They are the momentum sum rule: 



and the baryon sum rule 



/ f\ ^ P"~^$^(Md, k) + (Md - 2i?) A$2(^ ^2. 

J (27r)^2E3(MD-2E)2 \ Md ^ D' ^ ^ d j k d, 

That means that in frame of the Bethe-Salpeter formalism the baryon and momentum sum rules are 
different form of the same normalization condition (|4^). What gives solution of the longstanding problem 
of the simultaneous satisfaction of these sum rules in the framework of the phenomenological models of 



The Nonrelativistic Limit. Let us expand the energy of the bound nucleon in (|360D in powers of 
k^/m^. This leads to 



(2^ 



m 



-T 



m 



e dF^(x^) 



(362) 



where T = 2E — 2m is the nucleon kinetic energy and e = — 2m is the binding energy. 

We introduce the analog of the nonrelativistic wave function \E'^(k) (see also Eq.( p44|) ), related to 
$2(Md,A;) as 

{"^ (^^^'^)}.o=i.-Mo/2• 



^2(k) 



4E3Md(Mb -2E)2 
The normalization condition for \l/^(k) has the form: 

d^k 



{2ny 



^^(k) = 1. 



We compare (|362|) with the calculations in the nonrelativistic limit of the meson-nucleon field theory 
with synchronous nucleons: 



F^{xu) 



(2vr) 



m 



m 



-xb- 



dXB 



(363) 



Here ^^(k) is the solution of the Shrodinger equation ||136 . 

Equations ( |362| ) and ( ^631 ) obviously have the same structure. In the nonrelativistic calculation the 



term containing the derivative of the nucleon structure function arose as a result of the inclusion of the 
meson corrections associated with the nucleon potential, while the analogous term in ( p62| ) is ensured 
by the relative time of a bound nucleon. It is this contribution which makes the deuteron to nucleon 
structure functions ratio to differ from unity in the range 0.3 < x < 0.6. 

From the obtained results the x-rescailing model can be derived. Introducing the following variables, 



e=<T > -e = 2E - Mb and y 



xd Md {E - ks 



m 



we rewrite the expression ( p62| ) in the form: 



F^ixB) 



de j dy 





F^ 



y m 



dF^ 



y "I 



m 

E 



d{xY)/y) 
h 



X 



(364) 



m 



6{e - (Md - 2E)) 



The variable y has the meaning of the fraction of the deuteron 4-momentum projection onto direction of 
the photon momentum transfer carried by the struck nucleon. It is normalized on the nucleon mass in 
our notation. The variable e characterizes nucleon separation energy. If we assume that the mean value 
of this variable is small, then the expression in the square brackets can be considered as an expansion 
with respect to the parameter e. Thus we can convert it to the following expression 



xdMd 



+ 



dF^ 

e xb 2 



y m J 



y m J m y"^ d^x-o/y) 



xb Md 
y — - m 



To put this another way, the assumption is correct only when the states with high separation energy 
e give negligibly small contribution to the integral in Eq. ( P64| ) . Then the expression ( |362| ) takes the 
following form: 



m 



M, 



-X 



D 



oo M-o/m 



de 



y — e/m 



(365) 



Here x = x-oM-Q/m is the Bjorken x normalized to the nucleon mass, f^^^{y, e) is the deuteron spectral 
functions for a bound proton: 



These formulas precisely reproduce expressions obtained in the x-rescailing model ||139|| . It is clear that 



this expression is incorrect if the nucleon separation energy becomes large, which seriously constrains 
this model to small values of the parameter e. This fact can be important in calculations of the 
structure functions of heavy nuclei at large values of x and can explain why the x-rescailing model fails 
to reproduce experimental data at large x. In next section we will consider extension of our formalism 
for nuclei heavier then deuterium. 

Thus, the Bethe-Salpeter formalism allows the deuteron structure function to be expressed in terms 
of the structure functions of the bound proton and neutron. The antinucleon contributions are sup- 
pressed as the square of the mass defect. The inclusion of the dependence on the relative time in the 
amplitudes for DIS on bound nucleons leads to a modification of the nucleon structure reminiscent of the 
EMC effect in heavy nuclei. This allows us to conjecture that the nature of the effect can be attributed 
to the evolution of the bound nucleons relative time effect from nuclei with A = 2 to the nuclei with 
values of A at which the saturation of binding effects sets in. In the next section we will consider in 
detail the application of the Bethe-Salpeter formalism for light nuclei and provide an analysis of the 
relativistic effects in the evolution with A. 



5.4 Structure Functions of Light Nuclei and the EMC effect 

Here we shall study the derivation of the relative changes of the structure function, F^, in relation to 
the structure function of the isoscalar nucleon, = ^[F2{x) + where p and n denote the free 

proton and free neutron obtained from the recent world data fit. On the other hand, comparison with 
the experimental data can be made only for ratios of the structure functions of a nucleus A and the 



deuteron — A/D. This is why we also calculate the ratios A/D using the results of the section ^ 
In this way we will analyze also the evolution of the effect of the relative time in the bound nucleons, 
discussed in section |5.3|, from the deuteron to light nuclei and its manifestation in the EMC effect. 



5.4.1 Amplitude of the DIS for light nuclei 

We consider the generalization of the formalism developed in the preceding section for the analysis of 
DIS off n-nucleon bound system, n = 2-=-4. According to Eq. (|20D the nuclear Compton amplitude can 
be written in the form: 

T;^uiP, q) = j dlCdJC'T^iP, /C)5(„)(P, /C)G2(„+i)^^(g; P, /C, /C')%)(P, /C')r^(P, /C'), (366) 

where /C denotes a set of momenta which describes the relative motion of nucleons, JC = ki, . . . , 

dK. = d^ki/{2'7i)^ . . . d^kn-i/{27i)*, and P is the total momentum of the nucleus. The function r'^(P, /C) 

is the BS vertex function in momentum space: 

71 

5(„)(P, /C)rf (P, /C) = / d'x, . . . d^a;„e"^-i''"'$„,p(xi . . . (367) 



Here the Green's function G2(n+i)^^ represents Compton scattering of a virtual photon on a system of 
virtual nucleons. As in the case of the deuteron (see section |5.3|) , the only independent term comes 



from the relativistic impulse approximation, while irreducible interaction corrections to the imaginary 
part of are suppressed by powers of l/Q"^ ||135|| . This justifies the consideration of the zeroth order 



term of G 



G2(„+i)^,(g; P, /C) = J] G4^.(g; P, h) ® S^^'_,ik,, . . . h_^, . . . fc„_i)(5(/C - /C') + 0{1/Q^). (368) 



Then can be rewritten in terms of the off-mass-shell nucleon Compton amplitude Tj^^[ki, q): 

T^P.q) = I t//C5^r^^,(A;„gMk,)5(„)(P,fcO«(k.)r^(P,/C)5(„)(P,/C)r^(P,/C) (369) 

T^ih, q) = 7l(kOG4^,(g; P, h)u{\^i) (370) 

Following the procedure described in the previous section, we can make integration over fcjQ this 
expression, relating with on-shell nucleon Compton amplitude. 



Tl{p.q)=i I ciW^^(N|T(J,(x)J,(0))|N). 



However, this can be realized only after the singularities in nucleon propagators and the BS vertex 
functions are taken into account |p.35|| . Unlike the deuteron case, where singularities in the BS vertex 



function can be neglected, in case of A > 2 there are additional singularities in related with nucleon- 
nucleon bound states. These are poles in the range of low relative momenta, and we have to take 
them into account in the integration over fcjo- To this end, the "bare" BS vertex function can be 
introduced, which is regular with respect to the relative nucleon momenta ||53|| : 

r^(p, ^) = - j dJCg^niP, ^, /c')%)(p, }c')g{p, /C), (371) 

where g2n denotes the regular part of n-nucleon Green's function at P^ M\: 

n-l 

5(2n(P,/C,/C') = ^ G2ni{P,h---krn]k[---^'rn)®G2(n-rn){P,k^+l---kn,k'^+l---K)- (372) 
m=l,{l...n} 

The sum implies all the possible transpositions of bound particles. The function g2n, however, contains 
singularities of m-nucleon Green's functions {m < n). For example, in the case of ^He the function qq 
depends on the exact two- nucleon propagator G4, which contains the deuteron pole and the nucleon- 
nucleon continuous spectrum g^. 

c(^^ + kk kA - r°(2P/3 + fc,fci)r"(2P/3 + /c,fc;) /2P \ 

G,^- + k,k„k,j- (2P/3 + kf - Ml + ^4^^ + fc,fc„fc,J . (373) 

For ^He one has, additionally, the ^He and poles. Thus, unlike the deuteron case, the nuclear 
amplitude of DIS is determined not by the nucleon amplitude alone but also by amplitudes of all 
possible bound fragments of the nucleus. 

This can be demonstrated by the example of the He^ nucleus. Substituting expression ( P71|) into 
Eq. (|369|) and using the relation ( p35| ) we obtain the hadronic tensor of ^He in the form: 

r"(2P/3-Fi^,fc)rD(2P/3 + ir,fcO _ 3H 



^^^^ Gi{2P/3 + K, k, ki)S{P/3 + K/2 + k^) ® S{P/2> + K/2 - k^) 



3 Up 

G4(2P/3 + K, k, fci)] ® ^(P/3 - K)g «'=(P, K, ^) (p/3 _ _ • (374) 



Integrating over the zeroth component of the relative momentum of the fragments we obtain the hadronic 
tensor for ^He expressed in terms of physical amplitudes of the fragments and its derivatives over ko at 



The scalar structure functions can be extracted from the hadronic tensors with the help of the 
projection operators (p59|). Introducing now Bjorken variables for a nucleus = Q^/{2Pa ■ q) and for 



a nucleon = Q /(2P^ ■ g), where N is either p or n, we find F2 for He and H in the form: 
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2 [^'■^n) = ^2 la;3Hej|p^n 
and for "^He in the form: 
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where = —Ma + -En + Ea-i can be interpreted as the removal energy of the corresponding nu- 
clear fragment. The three-dimensional momentum distributions $^(k) are defined via the bare Bethe- 
Salpeter vertex functions. For example for ^He one has: 



<^>3He(k) 



mM, 
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4EpEDM3He(MD - Ep - Ed) 



d% d% -3He,„ , , ^o /^2P , , ^ 



(377) 



2P 



k,k,]T'' (^ + k,k'^S2(^ + k,k[ 



J2<i^H{k)]g'''%p,k,k[) 



ko=kop 



where /cop = M3h/3 — E^. Since at present there are no realistic solutions of the BS equation for a 
bound system of three or more nucleons, we have to use the phenomenological momentum distributions 
for numerical evaluations. 

The momentum distribution (377) describes the motion of a nuclear constituent (N, D, ...) in the 
field of the off-mass-shell spectator system. It is directly related with the nuclear momentum distribution 
measured in the e-A scattering when only a struck nuclear constituent is detected. It is reasonable, 
thus, to assume that the momentum distributions in Eqs. ( p75| ) and ( p76| ) can be related with those 
extracted from the experimental data. In the calculations we make use of the distributions available 
from 1 113711 and | ]138| |. The contribution arising from the continuous spectra {ppn for ^He and ppnn for 



^He) is small in the considered kinematic range and does not change comparison of the final result with 
the data. This justifies some simplifications which results in a rather transparent form of Eqs. ( 375 ) 



and ( P76|) . The contributions neglected in the derivations have been consistently taken into account in 



the normalization of the momentum distributions ^sjjg and $4^3. 

It should be stressed that both the modification of and its evolution from A 



1 to 4 found in 

the approach |135|] are the consequences of the relativistic nature of the nuclear structure. In 

the analytical calculations, it is essential to use the fact that the nucleons in the nucleus are separated 
by the relative time r^. It is this feature which is the cause of the nuclear effect in F^, which appears 
as a result of the dependence of the hadronic tensor of the bound nucleon on r^. 



5.5 Results 

We emphasize that both the modification of the F2 and its evolution from A = 1 to 4 obtained in the 
framework of our method result from the relativistic consideration of the nuclear structure. In the deriva- 
tions we essentially exploited the fact that the nucleons behave in a nucleus as asynchronous objects. 



of the bound nucleon hadronic tensor on Xj. The developed approach has a twofold merit. First, we 
can naturally reproduce the results of nonrelativistic models (e.g. ||139|| ) which offer the parametrization 
of the relativistic binding effects. Second, the outcome of the present study is particularly easy to 
understand when compared with the results of the x-rescaling model ||139| . 

As an example let us consider the structure function of ^He. If we regard the integrand in Eq. (|375|) 
as the terms of the expansion in the binding energy and take into account the fact that terms above 
the first order in this quantity are negligible, we can add higher-order terms in such a way that the 
resulting series can be represented by the expression. 



m 



-X 



3 He 



X 



y - e/Mo 



(378) 



Here e = Ap^*^ has the meaning of a nucleon (deuteron) separation energy, x = xsn^Msii^/m is the 
Bjorken x normalized to the nucleon mass, and f^^^^^^^'^iy, e) are the ^He spectral functions for a 
bound proton (deuteron): 
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Indeed, from the comparison of Eqs. (|359|) and ( |375|) with Eq. (|378|) , we find that the relative time 



dependence in the off-shell nucleon Compton amplitude results in the rescaling of the nucleon Bjorken x. 
We notice also that due to the relation between the nucleon mass and the four-dimensional radius of its 
localization region, ~ the dependence on x; has to lead to the increase of the localization region 

of the nucleon. In a way this result resembles the model considerations of the effect of the increase of 
the deconfinement radius or nucleon swelling ||140 |. 



The binding effects are expressed in Eq. ( [j75| ) and Eq. ( |376[ ) as the first-order derivatives of the 
structure functions of nuclear fragments. Thus the input structure functions, F2^^\x), are responsible 
here not only for the internal nucleon structure but also for the dynamics of the two-nucleon interactions 
as well. Similarly, F^ is responsible for the structure of the two-nucleon bound state and for the 
dynamics of three-nucleon interactions. The derivative of F^ is expressed in terms of the first- and 
second-order derivatives of F2 with the corresponding coefficients as shown in Eq. ( p62| ). Since the 
off-shell deformation of the bound deuteron structure is determined by the second derivative of F2 , this 
very term accounts for the three-nucleon dynamics. However, the second derivative of F2 contributes 
to ^2^"^ with a very small coefficient, A^^'^Ap , and the three-nucleon dynamics can thus be neglected 
in the consideration of the binding effects in DIS. 

The nucleon structure functions are introduced by parameterizations based on the measurements 
of the proton and the deuteron structure functions in DIS experiments. We have used the most recent 



parametrization of -Ff (x, Q"^) found in [ 141 | and fixed the value of to 10 GeV^. The structure function 
F^(x) is evaluated from -F|'(a;) and from the ratio F^(a;)/F|'(a;) determined in ||142|| . We have verified 
that the uncertainties in F2^^\x) are suppressed in the obtained ratio r^{x) and, thus, can be neglected 
in the considered kinematic range. On the other hand we have checked that an unrealistic input -F|^(a;) 
would have completely destroyed the evolution of the modifications we find in the lightest nuclei. 

The results of the numerical calculations are presented in Fig. |l8|(a). Here, we show how the free 
nucleon structure function -F^(a;) {A = 1) evolves to the deuteron {A = 2) and helium {A = 3 and 4) 
structure functions. The evolution, which starts from F^{x), is shown in Fig. |l^(b). In contrast to the 
modifications observed for nuclei with masses A > 4, the pattern of the oscillation of r^{x) changes in 
the range of A < 4, causing the coordinate of the cross-over point, 0:3, to move toward a larger value 
of X. 

The modifications with respect to -F^(a;) (Fig. 0(a)) are not only of academic interest. We use them 
to demonstrate that the change of the nucleon structure function in the deuteron cannot be regarded as 
negligible, and, therefore, the relation F,^{x)/F^{x) ~ F,^{x)/F2{x) cannot be justified. The position 
of X3 is displaced by 0.08 when F2{x) is replaced by F2{x) for the case A = 3 (Fig. ^(b)). The 
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Figure 18: (a) The ratio of structure functions F2 / ■ (b) The ratio of the structure functions / ■ 
The dashed curve in the left panel shows the result of the calculation for the deuteron [A = 2) made in 
section |5.3| . The results for A = 3,4 are shown by the solid lines. For A = 3 the isoscalar combination 
+ F;^]/2 is used. 



measurements of the ratios F,f{x)/F^{x) [ |132| |. According to |132| , X3 = 0.84 ± 0.01 independently of 



A ii A > 4. This accuracy allows one to reliably discriminate the effect of modification of the deuteron 
structure from that of the structure of the free nucleon. 

It is remarkable that the value of (1 — 0:3), which is found for F2{x)/F2{x) to be ~ 0.32, decreases 
for the ratios F^^^ (x) / F^ (x) and F2^'^ (x) / F^ (x) to ~ 0.16 and ~ 0.08 respectively. Further evolution 
of the modifications of F2{x) beyond A = 4 is forbidden by Pauli exclusion principle. As it follows 
from the pattern displayed in Fig. 0(a) and from the relation between the cross-over points X3, the 
modifications of the nucleon structure resemble a saturation-like process which is fully consistent with 
the rapid saturation of the nuclear binding forces. This phenomenon allows us to introduce a class of 
x-dependent modifications caused by the binding effects. Within this class there are no mechanisms 
which could lead to further changes in the pattern of r^(x) formed at the first stage of the evolution, 
A < 4. The evolution of the modifications to heavier nuclei, where the EMC effect was discovered, has 
to proceed independently of x and should be viewed as the second stage ||131| ]. The two-stage concept 
of the evolution of the free nucleon structure in nuclear environment is crucial for understanding of the 
long-standing problem of the EMC effect. 

As long as the experimental data for A = 2 and 3 are not available, our predictions can be only 
confronted with the results on F2^'^{x)/ F^{x) reported in Refs. ||143| , p,44|| . This comparison is shown in 



Fig. 0. The position of the cross-over point, obtained from our calculations is X3 = 0.913, which is in 
good agreement with the extrapolated data. It is of course of high importance to improve the accuracy 
of the data. 

On the other hand, we note particularly a good agreement between the corresponding point for A 



= 3, X3 = 0.845, and the average 0:3- value for nuclei in the range A = 9 ^ 197 found in ref. ||132|| . Such 
an agreement naturally follows from the two-stage concept of the F2{x) evolution which is x and A 
dependent for A < 4 and A dependent only for higher masses. The remarkable feature of our result is 
that the x dependent pattern of the EMC effect found experimentally in Fe develops already at A = 3 
and therefore can be viewed as a particular case of the class of the modifications. This is demonstrated 



by comparison of our calculation with data for A = 197 presented in Fig. Thus the EMC effect in 
heavy nuclei can be considered as the reflection of the relative time effect of bound nucleons scaled by 
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Figure 19: Results of the calculations of F2^'^(x)/F2°(x) performed in Refs 0, ^ (solid line). The 
experimental values are shown by the dark ||143|| and light ||144|| circles. 




Figure 20: Evolution of the nucleon structure in nuclei from A = 3to A = 197, obtained by multiplying 
of F^^^ (x) / (x) by a scale parameter. The A dependence of the parameter is determined by the 
Wood-Saxon potential and is shown by different degrees of shading. The results of measurement on a 
gold target were obtained in ref. 



143 



A fundamental relation follows from the obtained results. Since binding corrections for ^He and '^H 
have the same form (see Eq. ( |375| )) we can write 



1 





The result represents the Gottfried sum /, which has often been studied experimentally from the 
combination of -Fj'l^) -^2^(^) (cf- Ref. ||145|| ). Such a combination is equal to / to within a correction 



proportional to — 0)- Indeed, as follows from Eq. (R6 
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Apparently, such tests cannot be performed rigorously because is unknown at x = 0. On the other 

hand, if the difference of the '^He and binding energy is considered negligibly small, an experiment, 
which used these targets, would be able to measure the nucleon isospin asymmetry independently of 
the model uncertainties in the binding corrections. 



5.6 Conclusions 

We have shown in this section that the BS formalism can be applied to solve the topical problems of 
DIS of leptons off the lightest nuclei. In particular, the method for calculations of the evolution of the 
nucleon structure function as a function of A has been developed. The method allows us to express the 
structure functions F,f{x) of the bound nucleon in terms of the structure functions of nuclear fragments 
and three-dimensional momentum distributions. 

The pattern of the evolution of the modification of the nucleon structure function in the lightest 
nuclei, D, ^H, '^He and ^He, is consistent with the saturation property of the short range nuclear binding 
forces. The evolution is totally different from that observed previously for heavy nuclei, in which only 
the amplitude of deviations of F.^/F^ from unity increased with A. The quantitative predictions for 
^He and ^He nuclei, which have to be verified in future experiments at HERA or TJNAF, imply that 
the EMC effect in heavy nuclei can be naturally understood as distortions of the parton distributions 
in ^He or which are modified by the nuclear density effects. 



6 Summary 

In the present review we have considered the formulation of the Bethe-Salpeter equation. It is realized for 
the two-nucleon system by using the multipole expansion with the spinor structure of the two nucleons. 
The separable ansatz for the interaction kernel for each partial wave has provided a manageable system of 
the linear homogeneous equations for the BS amplitude. We have demonstrated then the construction of 
the separable interaction by taking only one term in the Yamaguchi form. Even with the two parameters 
for the ^5*0 and ^5*1 channels, we have found good reproduction of the phase shifts up to about 100 MeV 
in addition to the deuteron binding energy. This part has demonstrated the details of the formulation 
of the BS equation with the separable interaction. 

We have switched then to the case with the use of the covariant revision of the Graz II separable 
potential with the summation of several separable functions. The calculated results have been compared, 
first, to the static deuteron properties. The comparison shows very good agreement. We have applied 
then the BS amplitude for the calculation of the nucleon form factors that determine eD elastic scattering 
cross sections. Comparison of the obtained results with the experimental data is good in general. But 
there exist definitely some defects in the comparison with the data as the charge form factor and the 



The comparison with the elastic scattering has brought us to discuss the ingredient of the BS 
formahsm by taking the simple cases as the deuteron magnetic moment, deuteron quadruple moment 
and further the electro-disintegration of the deuteron. In this discussion, we have taken all the possible 
channels in the BS amplitude. We have found that the relativistic covariant description automatically 
include the meson exchange currents, in particular, the pair current through P-wave (negative energy 
state) component in the BS amplitude. We find it necessary to extend the partial waves to include the 
P-wave in order to construct the relativistic deuteron state. The P-wave components were not included 
in considered separable interaction kernels. 

Reactions of the elastic lepton scattering off the deuteron and the deuteron electro-disintegration 
served as a testing ground for the method under investigation and helped to outline both strong and 
weak points of the approach. The analysis has proved the technique to be very promising, even if we 
find a few evident discrepancies with data at this stage of development. Several items can be suggested 
for the program of further theoretical studies. 1) Construction of the separable potential of higher rank 
in order to reach better understanding of the properties of the deuteron, phases of NN- scattering, and 
to study hadron-dcutcron processes (for example, the reaction p + D ^ p + X); 2) Research into the 
relativistic two-body currents; 3) Studies of the off-shell effects in lepton-deuteron scattering. 

One specific feature of the BS formalism deserves a special comment. The BS amplitude depends 
on the zeroth component of the relative coordinate (relative time) of the bound nucleons, which is 
reflected in the dynamical observables of the n-nucleon bound state. In momentum space, this leads 
to the dependence on the zeroth component of the nuclcon relative momentum (relative energy). The 
dependence is manifested as observable effects in DIS of leptons off the lightest nuclei. 

The extension of the approach to the process of deep inelastic lepton scattering off the deuteron has 
been realized in a model-independent way. This quality is of particular importance for the consideration 
of the relatively small effects of the modification of the nucleon structure function P2^(^) by the NN 
binding forces. Based on the the model-free technique, the method for calculations of the evolution of 
the nucleon structure in the lightest nuclei as a function of A has been developed. 

We have found that the effects from asynchronous nucleons which naturally follow from the relativis- 
tic treatment of the two-nucleon bound state are decisive in obtaining differences between the structure 
functions of bound and free nucleons. The characteristic modification of the nucleon structure functions 
found for A = 2 serves as a priming for the modifications in the three- and four-nucleon systems and 
plays, therefore, a fundamental role in the evolution of the bound nucleon structure. The EMC effect, 
which was essentially the observation that partonic structures oi A = 2 and A = 56 nuclei were different, 
can be now regarded as a particular case of the whole class of modifications of the free nucleon structure 
in nuclear environment. 

When translated to a nonrelativistic language, the event of asynchronous nucleons can be associated 
with the increase of the localization region for the bound nucleon and is observed as the modification of 
F^ix). The developed approach does not require consideration of the three- nucleon forces to describe 
correctly the data available for the ratio F2^^/F^. The two-nucleon interactions can be, therefore, 
considered as the dominant mechanism for the evaluation of the nuclear binding effects in the kinematic 
range 0.3 < x < 0.9. 

The results obtained for the lightest nuclei, D, ^H, ^He and ^He, demonstrate that the BS formalism 
can have interesting practical applications for a certain class of electromagnetic reactions, namely DIS of 
leptons in asymptotic regime. We see future prospects of the covariant BS formalism in the confrontation 
of forthcoming high precision data with quantitative theoretical predictions. 

We have reviewed the covariant description of few-nucleon systems in the framework of the Bethe- 
Salpeter approach. This trial is still at primitive stage as compared to the widely used quasi-potential 
description of the BS equation. We thought it important though to write up the present status of this 
new development in the review form. Hence, the purpose of this article is not to tell all the consequences 
of the present approach but to show what is included in the covariant description of the BS approach 
and what has to be done further in this framework. In this sense, we would like to mention also the 
connection of the covariant BS formalism to the one of the light front formalism discussed in the article. 
Even at this simple stage it is clear that covariant four- dimensional approach give qualitatively new 
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Appendix A 

In this appendix we give the matrix elements used for calculation of the deuteron quadrupole moment. 
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Here we omit arguments {ko, \ k\) in the BS amplitudes ^35+ and 4>3jj+- 



Appendix B 



In this appendix we give the functions V^^'"^^ which defines the Lorenz invariant part V{s, q^) of the 
matrix element for the deuteron electrodisintegration. This part is defined by the Eqs.( |303|) and ( P04| ). 
The function used in the Eq.( p04|) are follows: 
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with the notations 

uji = 1/m, LJ2 = M^/m'^, U3 = l/m^, LJ4 = u)^ = M"^ /m^. 
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